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PREFACE 


This study analyzes the design of an experiment to 
measure communication characteristics of wideband satellite- to- 
ground links. Of special concern are the effects of rainstorms 
and atmospheric turbulance on path -attenuation and phase fluctua- 
tion. The proposed experiments are to aid substantially in the 
accurate determination of the tolerance of typical communication 
systems to the various interference producing mechanisms. 

The study estimates the expected range of channel parameters 
such as delay-Doppler scattering functions, range and frequency 
resolution requirements and dynamic range. Key hardware require- 
ments are considered and data collection and analysis methods 
are investigated. 

The study concludes with the recommendation of a modification 
of multi-tone techniques used in the ATS-5 and ATS-F experiments. 
Hardware analysis indicates advantages to the in-phase/quadrature 
approach of measurement implementation. 

Recommendations include development of key components such 
as space-qualified, stable reference oscillators, IMPATT oscilla- 
tors capable of operation above 20 GHz, and an efficient frequency 
synthesizer for the prober. 

The authors of this report are Howard Salwen who was project 
engineer and Dr. Peter Alexander who was concerned with the analy- 
tic portions of the study. Contributing to this investigation was 
Dr. P.A. Bello. The Goddard Space Flight Center technical monitor 
was Leonard Kleinberg. Dr. Julian J. Bussgang and John N. Pierce 
edited the final report. 
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SECTION 1 
INTRODUCTION 


1.1 Millimeter Wave Communications Systems 

It is rapidly becoming evident that the expanding spectrum 
requirements of various radio communication services will soon 
stimulate the development of millimeter wave communication 
links. The millimeter wave portion of the spectrum, particu- 
larly the band from 30 GHz to 90 GHz, offers the use of poten- 
tially wide bandwidths and high gain antennas, and appears 
eminently suitable for high data rate satellite -to -ground links. 
This is especially true when one considers the performance 
gains and communication efficiency which may be realized by 
the use of time division multiplexing (TDM) . Typically this 
approach results in very wide signal bandwidths thus necessi- 
tating a reexamination of the link characteristics, to deter- 
mine whether or not such bandwidths can be supported. 

Visions of the ultimate role of millimeter wave communi- 
cations systems have been outlined in the literature for both 
synchronous satellite -to- ground and terrestrial links [l. 1,1.2], 
In both situations the overriding concern is that of attenua- 
tion due to intervening rain storms. For line of sight ter- 
restial links this necessitates the use of reduced repeater 
spacings down to several miles to ensure adequate signal 
strengths. As a consequence wideband digital modulation with 
repeater regeneration becomes desirable if not necessary, to 
alleviate effects of interference and repeater noise. 
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Satellite systems must also be designed to take into 
account high attenuation levels due to rain in the atmos- 
phere. In most instances the beam center-lines subtend an 
angle of 10° or more with the horizon, so that the path 
traverses the atmosphere for only a few kilometers. Gener- 
ally speaking this results in less severe scattering and at- 
tenuation effects than would be encountered over terrestial 
links, but nevertheless a substantial deterioration over 
free space propagation conditions takes place. This report 
is concerned exclusively with the modeling and analysis 
of satellite-to-ground data links. 

There are several new problems encountered with the use 
of both higher carrier frequencies and wider bandwidths. 

For example, the move to higher frequency bands involves sub- 
stantial increases in path attenuation and antenna beam dif- 
fraction effects, while the frequency dependence of the medium 
and multipath scattering result in an effective band limiting 
of the channel. The basic objective of this study is to 
investigate the link parameters which determine communication 
performance, and to outline various measurement techniques 
which might be used to evaluate such parameters in a milli- 
meter wave experiment. Typical parameter values and their 
range of variation, are to be estimated from the experiment 
and related to the local climatological conditions where pos- 
sible. The link characteristics of interest, which must be 
examined to allow prediction of communication link perfor- 
mance, can be grouped in the following way. 
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(a) Signal Fades 


These will occur on a short time scale (seconds) due to 
multipath scattering, and on a larger time scale due to gross 
variations in the intervening medium. The depth and frequency 
of fades must be determined so that an adequate fade margin 
can be specified in the link design. 

(b) Pulse Distortion 

Because the signal transmitted between satellite and ground 
station traverses various layers which have frequency depen- 
dent refractive index properties, wide bandwidth signals are 
subject to an effective filtering and distortion. This will 
be particularly evident when a short pulse is transmitted. When 
the pulse length is reduced sufficiently, the original trans- 
mitted pulse becomes smeared in time. This effect severely 
limits the transmission of wideband data primarily because of 
inter symbol interference. 

(c) Multipath Spread 

When the medium has random structure, (e.g., atmospheric 
turbulence and rainstorms) the transmitted energy is scattered, 
thus attenuating the incident field. More importantly, the 
receiving system can then acquire signal components from other 
than the direct path. These scattered signals arrive with a 
variance of diverse time delays and contribute to further smearing 
of the signal as in topic (b) . The principal difference here 
is that the scattered components are random in their behavior, 
hence requiring a characterization based on statistical 
properties. 
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(d) Doppler Spread 


The scattering of energy from inhomogeneities along its 
path gives rise to a time varying distortion of the transmitted 
signal. The scattering medium is continually in motion so 
that individual scattering elements impose Doppler frequency 
shifts on the received components. The resultant spectral 
properties of the signal are thus determined by bulb and tur- 
bulent motion of the medium. Gross spectral parameters such 
as RMS Doppler spread of a transmitted sine wave are impor- 
tant to the communication engineer, since certain communica- 
tion techniques are more susceptible to Doppler spreading 
than others. 

It is apparent therefore that before the millimeter wave 
region can be seriously considered for high rate data links, 
a carefully designed characterization experiment must be im- 
plemented. Although such an experiment is closely coupled to 
studies of the physics of the atmosphere, it is essential 
that certain parameters which are important from the point of view 
of the communication problem, be measured. 

A second area of research which will ultimately be of 
great importance is that of co-channel interference. By 
taking advantage of PCM modulation techniques, which are in- 
herently less sensitive to interfering signals then conven- 
tional analog methods [1.3], and by utilizing the narrow 
antenna beamwidths in the millimeter wave region, it appears 
that dense communication networks will be feasible. Such 
networks, which will involve considerable overlapping of the 
frequency band assignments, are likely to be interference 
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rather than noise limited. There are various ways that energy 
from co-channel transmitters can be coupled into a given link, 
and the prediction of the resulting interference levels is 
obviously of an important problem. Some attention is given in 
Appendices C and D to interference which arises due to coupling 
of antenna beams. 

1.2 Phenomena to be Observed 

The transmission medium between satellite and receiver 
will be subject to marked variations due to changing weather 
conditions. Two of the most important atmospheric effects 
which distort the signal are associated with rainstorms and 
atmospheric turbulence. The former will be the one 
determining the worst case system performance but the latter 
is important because it will be the dominant propagation mech- 
anism most of the time. 

The channel characteristics are determined by the nature 
of the propagation medium and geometrical factors such as 
antenna pointing angles and beamwidths. 

The received field for the millimeter wave channel can 
conveniently be thought of in terms of the sum of two com- 
ponents; a scattered incoherent signal, and a coherent but 
time invariant component. This can be clearly demonstrated 
for scattering by raindrops, since the exact solution for the 
field scattered by dielectric spheres can be stated explicitly 
[1.4]. Similarly, many theoretical analyses have been pre- 
sented on the scattering effects of atmospheric inhomogeneities 
[1.5,1.6]. In both cases we are primarily concerned with a 
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statistical description of the field which results when a 
plane wave is incident on the medium. The total channel ef- 
fect can be conveniently summarized by the use of a time vary- 
ing transfer function H(f,t) such that, 

H(f , t) = T (f)[l + T ( f , t ) ] 

C K 

i.e., the response at time t to a sinusoid of frequency f. 

The function T c (f) represents attenuation and phase 
shift characteristics of the medium as it affects the coher- 
ent transmitted component. It appears as a function of fre- 
quency because the clear atmosphere and rain storms both be- 
have as dielectric material with a refractive index which is 
frequency dependent over very wide frequency bands. This 
gives rise to distortion of the broadband signals [1.7, 1.8]. 
Note that pure constant attenuation or linear phase shift 
characteristics with frequency do not introduce any distor- 
tion effects; the received signal is a delayed, attenuated 
version of that transmitted. 

The time varying filter T (f , t ) represents the scattered 
signal components. Its behavior as a function of t is deter- 
mined by movement of the scattering elements with time (e.g., 
downward or turbulent motion of raindrops). The spectral 
width of T is usually referred to as the Doppler spread. The 
frequency dependence of T is related to the geometrical pro- 
perties of the scatterer distribution, which also determine 
the multipath spread of the channel [1.9]. Note that T R (f,t) 
is a random function of both f and t. From a two dimensional 



correlation analysis, both multipath spread and Doppler spread 
parameters can be obtained [l.lO]. 

Provided sufficient information can be extracted about 
H(f,t), the communication link performance for any arbitrary 
signaling scheme can be computed. In general H(f,t) is a 
complex process, i.e., it contains both phase and amplitude 
information; therefore channel measurement techniques based 
on incoherent processing are often inappropriate. 

1.3 State of Knowledge Summary 

The current body of knowledge pertaining to millimeter 
wave propagation can be separated into three distinct areas: 
theoretical, related terrestrial experiments , and satellite- 
to-ground experiments. 

The theoretical area has been the most active one . The basic 
theory of propagation through a random medium has long been a 
subject of interest, particularly for dielectric sphere scat- 
tering [1.4] and propagation through a turbulent atmosphere 
[1.5,1.6]. Application of these stochastic wave propagation 
theories to the specifics of millimeter wave links has also 
been pursued [1.7,1.8,1.11-1.15] with major emphasis placed 
on relating attenuation characteristics to meteorological con- 
ditions. This is true also for the second category which 
covers experimental programs [1.16,1.17]. However several 
experiments have been developed with the objective of obtain- 
ing data which would be useful from a communication theory 
point of view [1.19-1.21]. 
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The direct measurement of satellite to ground links has 
received only secondary attention. Most of the effort has 
been associated with the ATS-5 and ATS-F millimeter wave ex- 
periments [1.22-1.24]. 

At this date the total available data relevant to mil- 
limeter wave communication system design is much too 
sparse. Furthermore the link characteristics above 30 GHz can 
only be obtained by extrapolating this sketchy data, clearly 
pointing to the need for the implementation of a pilot exper- 
iment which can be easily adapted to a range of frequency 
bands. In designing an experiment it is useful to know at 
least the broad range of parameters to be encountered, so the 
following numerical values which have been extracted from the 
literature, are listed to help establish design guidelines. 

Attenuat ion: 

Atmosphere (Zenith) - 0.3 dB (total) 

Rain Storm (1 inch/hr) - 8 dB/km 

Multipath Spread < 1 nanosec. 

Bandwidth of T (f) >1 GHz 

c 

Doppler Spread < 1 Hz 

Scattered Signal/Coherent Signal - -40 dB (6"/hr rain, 

30 GHz) 

The parameters of course vary with antenna beamwidths and 
pointing angles ro that the values given above must be used 
with caution. They were obtained for narrow beam receiving 
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antenna (0.2°) pointed at the zenith. The multipath and 
Doppler spread parameters increase with increasing antenna 
beamwidth. 

1 . 4 Experimental Techniqu e 

In Section 1.2 an equivalent filter model for the link 
was presented. The design of the experiment is based on this 
model and its relationship to the basic physical phenomena. 

The filter characteristics which must be determined are phase 
and amplitude for T c (f) and two dimensional correlation pro- 
perties for T (f,t). 

R 

The principal difficulty involved with such measurements 
is the very large bandwidth required for the probing signal. 
There are two different philosophies that one could adopt - 
either measure channel characteristics with a signal of band- 
width comparable to that available for present day data modu- 
lation techniques, or alternatively probe with signal whose 
bandwidth is well in excess of the expected channel bandwidth 
even though the full band may not be utilized in contemporary 
communication systems. The former simply gives a verification 
that the channel properties will be suitable for immediate 
needs, while the latter will demonstrate constraints for future 
systems, and perhaps provide useful scientific data about pro- 
pagation through the atmosphere. 
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The choice of probing signals lies mainly between a fre- 
quency domain approach, i.e., the use of several tones centered 
about some nominal carrier frequency (multitone signal) , and a 
time domain approach, i.e., the use of a short pulse or pseudo 
random noise phase modulating sequence [1.25, 1.261. The ob- 
jective of covering at least 1 GHz bandwidth is difficult to 
achieve using the time domain methods. On the other hand, the 
multitone spectrum is easily generated and can be processed as 
a set of narrowband signals at the receiver. However, some 
problems do arise in connection with phase stability of the 
transmitter and movement of the satellite as will be discussed 
subsequently. 

Consider a transmitted waveform consisting of a set of tones 
Cs n (t)); 

s (t) = a sin (u) t + 8 ) n=0,+l, . ... ,+N (1.1) 

n n n n 

This type of signal structure is employed in the ATS- 5 and ATS-F 
propagation experiment. 

After phase and amplitude distortion by the channel, the 
received signals essentially represent the transfer function H(f,t) 
sampled at frequencies { w ^/ 2 ti } i.e.. 


r 

n 


(t) = |H ( 




sm 


( 


uj t + 0 

n n 




( 1 . 2 ) 


Thus the required properties of T (f) and T R (f,t) can be ex- 
tracted from the phase and amplitude behavior of the received 
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set of signals {r n (t)]. The received spectrum is a set 
of Doppler broadened spectral lines as shown in Fig. 1.1. 

Unfortunately the phase processes {8^} do not remain 
constant due to transmitter oscillator instabilities and satel- 
lite motion. It can easily be demonstrated that oscillator 
instabilities give rise to phase errors which increase lin- 
early across the band, and hence are indistinguishable from a 
variable propagation delay. The satellite also suffers from 
systematic and random perturbations in its position thereby 
introducing an additional variable delay component. These two 
effects can be lumped together as in Fig. 1.2. Thus the fil- 
ter T (f ) , representing coherent effects, is related to path 
delay perturbations in the line of sight propagation. 

It has been shown that the usual carrier phase locking 
procedure is inadequate for this kind of wideband measurement. 
That is to say, the correction for phase variations at center 
frequency leaves significant residual phase errors in the other 
received tones, and the magnitude of these errors is large 
enough to obscure the true phase properties of H(f,t). It is 
desirable to separate out the phase instability effects from 
channel induced perturbations so that the data is applicable 
to other satellite configurations. 

The time varying delay t introduces a phase shift versus 
frequency expressed as; 
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Transmitted Tone Output Spectrum 

Spectrum i 
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Figure 1.1 Frequency Domain Channel Response to a Multi-Tone Spectrum 
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Incremental Phase Correction Procedure 









0(f) = 2nf t 


= 2TT(f Q + A f ) T 

= e o( 1 + if) <1 ' 3) 

where 0^ is the phase at the carrier (center) frequency. Thus 
knowledge of 0 Q allows the receiver to make the required phase 
correction as depicted in Fig. 1.2. However the transfer func- 
tions T and T both distort the phase characteristics further, 

R C 

so that 0 q is never observed directly. This aspect of the mea- 
surement has been studied in detail in Chapter 3. In particular 
the following points are relevant; 

(a) When the scattered signal level is sufficiently low 
it can be shown that a suitable estimate of 0^ can be obtained. 
The detailed performance and precise region of validity must be 
studied further for various link configurations. 

(b) The phase instabilities and channel phase character- 
istics have different spectral properties in general. Our 
experience in the design of phase locked loop receivers leads 
us to believe that the two effects can be separated by suitable 
choice of filters within the loop structure. 

Returning to Eq. (1.2) it can be seen that a phase cor- 
rection scheme as discussed above, allows us to remove the phase 
parameters (.0 ) by suitable receiver design (see Fig. 1.2) « 

Hence phase and ; mj_ litude characteristics of H(f,t) are obtained 
directly. These quantities provide a complete description of 
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the instantaneous transfer function H(f,t) sampled in frequency. 

If the samples are taken close enough together, the complete 
function can be reconstructed. Since H(f,t) is a random pro- 
cess, its statistical properties are of interest rather than 
its instantaneous values. This aspect of the measurement is 
discussed more fully in Chapter 2, while Chapter 6 contains 
information on the required form of processing for the tone 
responses . 

1.5 Organization of the Report 

As just noted Section 2 contains a discussion of the channel 
characteristics and measurement philosophy. The discussion 
emphasizes the communication aspects of the channel's charac- 
teristics rather than the physical mechanisms that give rise 
to those characteristics. The discussion of Section 2 leads to 
the expression of channel characteristics in terms of parameters 
that can be measured or derived from an experimental configuration. 

Chapter 3 discusses one candidate for a channel measurement 
configuration, namely a multi-tone channel prober. This technique 
is currently employed in ATS- 5 and ATS-F K-Band experimental equip- 
ment. The discussion of Section 3 reviews the characteristics of 
such channel probing equipment and shows that the capabilities of 
such equipment can be augmented by the use of an incremental phase 
correction technique. The phase correction technique allows for 
measurement of the channel parameters in a manner which is more 
easily related to the channel characterization parameters in 
Section 2. 

Section 4 discusses an alternate probing technique, namely 
the pulsed probing technique* Several variations of the pulsed 
probing technique are described and discussed. Although at first 
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glance the pulsed probing techniques appear to offer advantages 
over the multi-tone techniques , it is found that component limita- 
tions make this technique less desirable. The pulsed techniques 
are desirable because the data derived through their use is in a 
form more suitable for deducing communication properties of the ■ 
channel. 

Section 5 presents detailed block diagrams for a multi-tone 
probing transmitter and receiver system. The system is designed 
to capitalize upon advances in the state— of— the— art of component 
technology which have occurred s ince the initial design of the 
ATS experimental package. The design furthermore# incorporates 
some of the data processing features derived in this report — 
namely# incremental phase correction and some minor modifications 
in the data detection circuitry. Assuming that it is undesirable 
at this time to make any modification to the spacecraft hardware 
for the ATS-F flight experiment, a modified receiver configuration 
is suggested which is compatible with the ATS— F K— Band multi-tone 
transmission, but which provides incremental phase correction at 
the receiver. Thus, if such an approach were taken, additional 
data could be extracted from the existing K-Band downlink signal. 

Section 6 contains a discussion of data processing requirements 
for a K-Band satellite-to-ground link channel measurement experi- 
ment. A data processing technique which deals with in-phase and 
quadrature components of the received tones rather than amplitude 
and phase components is presented and analyzed in detail. After 
describing the data processing algorithms to be employed, a sta- 
tistical analysis of expected performance using the in-phase and 
quadrature algorithms is carried out. Finally, the measured 
parameters are related to the simplified channel model which was 
derived in Section 2. Section 6 goes on to consider meteorological 
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data collection and processing. Rainfall rate, atmospheric 
turbulence and refractive index profiles are considered. 

Section 7 presents a brief summary of the results of 
this study effort. In summary, the report concludes that 
the multi-tone channel probing method remains the most desir- 
able candidate technique at this time in spite of the great 
advances in the state-of-the-art since the ATS experimental 
configuration was first conceived. The report suggests minor 
modifications in data extraction and processing techniques 
which will provide additional data about the channel. Block 
diagrams of a new transmitter and receiver are presented 
which incorporate these minor design advantages and which 
capitalize on the advances in the state-of-the-art in com- 
ponent technology which have occurred in recent years. The 
discussion of Section 7 includes some recommended development 
areas which will ensure success of future programs in this 
area. 
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SECTION 2 


CHANNEL CHARACTERIZATION AND MEASUREMENT PHILOSOPHY 

2 . 1 General Channel Characteristics 

The statistical cha racter ist ics of the millimeter wave 
channel are determined by the nature of the propagation medium 
and geometrical parameters such as antenna beamwidths . In the 
sections that follow two specific media will be considered in 
detail, i.e., rain and turbulent atmosphere. These will il- 
lustrate the general nature of the channel characterization 
task and with slight modification the analysis can be used 
directly for other meteorological conditions. 

2.1.1 Description of the Model 

The medium, which gives rise to scattering and absorption 
of millimeter waves, can not be characterized too accurately 
because of the natural inhomogeneous structure of rain, fog 
and other atmospheric phenomena. Instead an idealized model 
must be employed to gain insight into the general behavior. 
Figures 2. 1 and 2.2 indicate distinctions between the real and tne 
ideal satellite-to-ground communication links. The following 
assumptions are also necessary to simplify the description. 

(a) The distance R Q is much greater than the altitude 
of the scattering medium. (Typically the latter is 0 - 15 Km 
for atmospheric effects.) Moreover, the receiving beam is 
sufficiently narrow and the transmitting beam on the satellite 
so broad, that the field incident at the top of the layer can 
be approximated by a plane wave over the region which will ul- 
timately be of interest (several receiving beamwidths). 



Zenith 

Path 
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Fig. 2.1. Typical Path Geometry. 






(b) Oblique pointing angles may be of interest. These 
will be accounted for by increasing the effective thicKness of 
the layer using an appropriate factor, i.e., 


T 

sin 8 


( 2 . 1 ) 


where T^ is the physical thicKness of the layer and 8 is the 
elevation angle. Hence we need only consider the zenith 
pointing angle in this discussion. 

2 . 1.2 Discussion of Terminal Characteristics 

The received field for the millimeter wave channel can 
conveniently be thought of in terms of the sum of two compon- 
ents; a scattered incoherent signal, and a coherent time 
invariant component .This is clearly demonstrated in Appendix A 
for scattering by raindrops, where the exact solution for the 
field scattered by dielectric spheres can be stated explicitly. 
This is true also for scattering by atmospheric inhomogeneities, 
where the coherent field is simply the time invariant statistical 
average field. In both cases we are primarily concerned with a 
statistical description of the field which results when a plane 
wave is incident on the medium. In particular typical parameters 
describing the antenna terminal voltage are obtained as a function 
of various atmospheric conditions and link geometries. 


Of additional interest are the effects of multipath spread- 
ing, i.e., the time smearing of a short transmitted pulse, and 
temporal fading characteristics. Both of these are determined 
by scatterer properties and system geometry, and some details 
are available in the Appendices. The total channel effect is 
conveniently summarized by the use of a time varying impulse 
response h(t, t) , the response at time t to an impulse t seconds 
earlier. This function will always provide a valid description 
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of the channel since it is founded only on the assumption of 
linearity between input and output; this assumption is in turn 
assured by the linearity of Maxwell's equations which describe 
the received field. 


The channel can be described equally well by the time 
varying system function H(f,t) which is related to h(t,T) by 
Fourier transformation on t, i.e., 


H( f , t ) 



h( t , t ) e 


- j2rrf t 


dT 


( 2 . 2 ) 


H(f,t) then, is the complex time varying response to a tone 
of frequency f. Note that variation of H(f,t) with respect to 
t will be quite slow (tens of Hz or less) compared with 
the data rate, so that H is easily visualized in terms of 
conventional system functions, except that it is randomly 
varying in time. 

The measurement of channel properties can be approached 
at several different levels. We could attempt to measure the 
complete function H(f,t) over particular time and frequency 
bands. However such Knowledge would generally provide an 
overly detailed description of the channel: To provide data 

useful in predicting communication performance we need broad 
trends rather than detailed sample function behavior. Generally 
a sufficient description of the channel is provided from the 
average properties of H(f,t). When the channel statistics are 
Gaussian it is only necessary to measure certain underlying 
correlat ion functions, but even if the channel is non-Gaussian, 
correlation functions may still be appropriate for most pur- 
poses. Note that a given correlation function describing sta- 
tistical properties of H(f,t) will only be valid over time 
periods for which the bulK medium properties are fixed. 
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2 . 2 An Equivalent Filter Channel Model 

The effect of the atmosphere and discrete scattering 
particles (such as raindrops) on the transmitted signal, can 
be viewed in terms of equivalent filters. In particular we 
need to consider a model with two different Kinds of filters, 
one corresponding to the effect of the medium on the coherent 
signal, and the other corresponding to the scattered signal. 
Consider the former which will be denoted by T c (f). As the 
number and shape of constituent scattering elements within 
the coherent scattering zone changes, so does the equivalent 

refractive index of the medium, but as we have seen in Appen- 
dix A the effect is small. If the medium scattering pro- 

perties were invariant with the frequency of the incident 
field, we could assume that the equivalent refractive index was 
constant with frequency. Hence T^Cf) would have a constant 
amplitude and linear phase shift across the band. A refined 
model would allow variable attenuation and quadratic or higher 
order phase shifts across the band. These characteristics/ 
which have been investigated for rain in a theoretical study 
by Crane [ 2-1], can be expected to give rise to pulse distortion 
if a wide enough bandwidth is used. Note that if T c (f) was 
expanded as a function of frequency using a Taylor's series 
for phase and amplitude, the sample calculations of Appendix B 
suggest that only the first few terms would be required. 

With a single scattering model the energy received from 
a collection of scatterers is computed using only the original 
field incident on each scatterer (i.e., the field which would 
exist in the absence of other scatterers). One can visualize 
that this incident field, before reaching the scatterer, 
initially passes through a region which imposes its attenua- 
tion and phase characteristics in the same way as previously 
discussed for the coherent field. Similarly the field originating 
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from a single scatterer must traverse the medium between the 
scatterer and the antenna. Thus we can use a random time vary- 
ing filter T' (f,t) to represent the total influence of the 
incoherent scattered signal as in Fig. 2.3(a), or alternatively 
separate out the coherent distortion effects which result from 
propagation through the medium, as in Fig. 2.3(b). A channel 
model which is frequently of use in describing communications 
links when scattering is involved, is the Wide Sense Stationary 
Uncorrelated Scattering (WSSUS) model of Bello [2—2]. As indi- 
cated in Appendix A it is not necessarily valid to assume that 
the scattered component represented by T (f,t), fits into this 
category. The use of the WSSUS model requires system functions 
which are statistically stationary in both the f and t variables, 
and it is therefore clear that the product T R (f,t) T c (f) does 
not fit this category. 

In the material that follows we shall generally refer to 
the model given in Fig. 2.3(b) without assuming T to be a WSSUS 
filter. The overall input-output relationships can be described 
in terms of a composite transfer function H(f,t); 

H ( f , t) = T c (f)[l + T r ( f , t ) ] (2.3) 

Measurements of phase and amplitude or in-phase and quadrature 
components can be made using either tones or wideband signals 
and these provide data for the subsequent compilation of statis- 
tical properties. Because the transfer function T (f,t) has been 
defined as a zero mean process (see Appendix A for details) , it 
is possible to extract information about T c (f) by directly averag- 
ing the in-phase and quadrature components of a given tone re- 
sponse. Second order properties of H(f,t) are also of interest, 
but are more difficult to obtain, so that this subject will be 
left for a more complete development in later sections. 
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Fig. 2.3(a). Equivalent Filter Model. 
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2 . 3 Statistical Properties of the Receive d Signal 

2.3.1 Definition of Parameters 

The time varying impulse response h(t,T) is derived from 
the super-position of many scattered signals for each t and t. 
The scattering elements may be phase coherent, but there is 
randomness associated with their scattering properties, e.g., 
because of variations in shape and diameter. This leads to 
the assumption of a Gaussian statistical description for the 
impulse response h(t, t)., its transform on t, H(f,t), or for 
that matter the response to any deterministic waveform. The 
response to a tone of frequency f can be expressed as; 

r(f,t) - A(f,t) cos ( 2rrf t +Y(f,t)) 

= X(f,t) cos 2TTft - Y(f,t) sin 2nft (2.4) 

where 

X(f,t) = x(f,t) + X(f) 

Y(f,t) = y(f , t) + Y(f) { 2 . 5 ) 

Thus the quadrature components X(f,t) and Y(f,t) have been de- 
fined as the sum of a frequency dependent mean and a zero mean 
random term which is consistent with previous definitions in 

terras of specular and scattered signals. 

The time varying transfer function H(f,t) can be related 

to r (f, t) using, 

r (f , t) = Re (H(f, t)e j2TTft t (2.6) 

T (f) , T (f,t) and hence H(f,t) have been implicitly defined 
c R 

_oo < f < oo (although generally only H(f,t) f > 0 will be 

required since H(f,t) = H* (f, t) ), which implies that, 
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X( f , t ) = Re { H( f , t ) 3 


Y(f,t) = Im { H( f , t ) } (2.7) 

Thus 

H(f,t) = X(f,t) + jY(f,t) (2.8) 

In much of the material which follows, the properties of H(f,t) 
will be examined at a particular value of t. However for the 
present it is convenient to define the following statistical 
properties in terms of the complete time-frequency function 
r(f,t). The in-phase and quadrature components, which are 
assumed to be stationary in time, have mean and correlation 
properties denoted; 


E (X( f , t ) ] = x:(f ) 

(2.9) 

E {Y( f , t ) ] = Y ( f ) 

(2.10) 

E{x(f 1 ,t 1 )x(f 2 ,t 2 )} = R x (fi> f 2 ’ W 

(2.11) 

E{y(f 1 ,t 1 )y(f 2 ,t 2 )] = R y (f 1 .f 2 ' t l“ t 2 ) 

(2.12) 

ElxUj.tpyU-j.t.,)} = R xy (f l’ f 2’ t l' t 2 ) 

(2.13) 

Referring to the original definitions (Eqs. (2.4) 

and (2.5)) it 


is apparent that one can define a pair of covariance functions 

for the complex process H(f,t) and express them in terms of 

the quadrature component correlations above. For example define, 
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K H (£ 1 , f 2 ,t 1 - t 2 ) I EtWfj.tj) - - H*(f 2 ,t 2 )]J 

(2.14) 

% 

H denotes mean value. Using Eq.(2.8) and subtracting out the 
mean values we find; 




= EtWf^t^ - jy(f 1 ,t 1 )][x(f 2 ,t 2 ) + jy(f 2 ,t 2 )]) 
= V f l’ f 2’ t l- t 2 ) + R y (f l’ f 2’ t l" t 2 ) 

•+ Jt R xy (f 1 ff 2’ t 1 " t 2 > “ R xy (f 2’ f l’ t 2- t l )] 

(2.15) 


Similarly without the conjugate in Eq. (2.14) we get; 

/ 

V f l’ f 2’ t l" t 2 ) = ECtmfi,^) - - H(f 2 ,t 2 )]} 

= ECtxtfj^j) - jy(f 1 ,t 1 ) ][x(f 2 ,t 2 ) - jy(f 2 ,t 2 )]3 

= R x ^ f i > f 2’ t l“ t 2* “ R y^ f l ,f 2 ,t l" t 2^ 

- j^ R xy^ f l’ f 2 ,t l _t 2^ + R xy^ f 2’ f l ,t 2 _t l^ 

(2.16) 


Thus if the quadrature and in phase covariances (Eqs.(2.11) - 

(2.13)) are Known, both K R and K R can be reconstructed. Conversely, 

these two functions alone provide sufficient information to 

allow the extraction of R R and R uniquely. 

x y 
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In the next section we examine more closely the various 
correlation functions which have been defined and relate their 
properties to certain physical properties of the millimeter 
wave channel. 

2.3.2 Correlation Function Properties 

It is clear that estimates of the functions R^^, ant ^ R X y 
along with knowledge of T^, allow one to predict the performance 
of any potential’ communication scheme. For the present we shall 
be concerned only with measurement of the above functions, and 
in order to capitalize on any simplifying channel properties, 
we need to study them further. 

Consider first the question of stationarity in time. It 
can be shown that the response r(f,t) at a single frequency, 
which will be denoted by r(t) has a covariance function, 


• (t l’ t 2 ) = \ Re C K H (f i' f 2' t r t 2 )e:i2TTf(tl " t2) + K H (£ l ,f 2' t r t 2 )e52TTf(tl 2> 1 

(2.17) 


K 


Here we see that the bandpass process r(t) will be a 
function of .(t^+t^ and hence nonstationary unless K H is zero 
despite the fact that the in-phase and quadrature components 
were themselves assumed stationary. Returning to Eq. (2.16) 
which defines K u it is seen that stationarity of r(t) (for 
fixed f) requires; 

R x (f,f,t 1 -t 2 ) = R y (f,f,t 1 -t 2 ) (2.1i 


and 


R (f.f.t.-t,) = 

xy ’ ’ 1 2 


-R (f,f,t_-t.) 
xy ’ ’ 2 1 


(2.19) 


T> 
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These constraints have frequently been observed and satisfied 
in certain other problems concerned with bandpass stochastic 
processes. In particular they hold when r(t) is bandpass noise 
or when we can associate with r(t) a random initial phase com- 
ponent y which is uniformly distributed over [-tt,tt]. This 
o 

latter approach must be applied with care however. If the 

parameter is in fact constant over a given measurement interval 

(0,T), then subsequent processing of the data for that interval 

will not effectively average over y . Generally speaking it can 

be shown [2-3] that a large number of channels, e.g., tropospheric 

and ionospheric scatter channels, do give rise to a signal r (t) 

such that the above conditions are satisfied. These channels can 

be characterized in terms of K only, since K = 0. 

H ri 

A further simplification can often be made if the spectrum 
of the process r (t) is symmetrical about the carrier frequency. 

It can then be demonstrated that; [2-4] 

R xy< f ' t ' t r t 2> “ ° 

so that r(t) is completely characterized by a single correlation 

function R . Suitable extensions can be made to the above conclu- 
x 

sions for a pair of responses r(f^, t^) and r(f^, ^’ 2 )' ^ Ut t ^ le 
necessary symmetry properties (i.e. symmetrical delay power spec- 
trum) are rarely satisfied in practice. 

For the millimeter wave link it is expected that the chan- 
nel characteristics will preclude the assumptions of true station- 
arity for the response to a single tone, because the requirements of 
Eqs . (2.18) and (2.19) will not hold. The in-phase and quadrature 
components will _n general be strongly correlated and may have 
different variances. Evidence to support these statements is 
contained in Appendix A for rain scattering while the effects of 
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atmospheric turbulence have been covered in detail by Bello 
Ref. [2-5]. Broadly speaking we can expect the in-phase and 
quadrature components to be correlated and of unequal variance 
whenever the majority of the scattered energy originates from 
the first few Fresnel zones. As a consequence of this coupling, 
the resultant tone responses cannot be classified simply as 
having Ricean amplitude distributions; the behavior of the chan- 
nel is instead determined by the complete set of functions given 
in Eqs . (2.9) to (2.13) (see Ref. [2-6] for an evaluation of the 

resultant amplitude distribution) . 

These channel characteristics, which arise because of the 
extremely narrow beamwidths envisioned for millimeter wave links, 
obviously make the required measurement program more complex, 
and underline the need for coherent (in-phase and quadrature) 
receiver processing. The subject of channel measurement is 
broached in Section 2.4, while Sections 3 and 4 give a more de- 
tailed description of two different schemes. 

In Section 2.3.3 we present a discussion concerning the 
channel behavior for situations where the random signal component 
is much smaller in strength than the coherent signal. 

2.3.3 Channel Properties for Low Scattered Signal Levels 

The existence of a low level scattered signal allows the 
use of certain approximations in the determination of a suitable 
measurement strategy. The composite time varying transfer func- 
tion is. 


H (f , t) = (1 + T (f , t) ) T (f ) (2.20) 

c 

where T R is the transfer function corresponding to the scattered 
signal component, and T^ which includes the gross path loss 



characteristics, is representative of the average medium behavior. 
When the ratio of unscattered power to scattered power is high 
(say greater than 20 dB) we have, 


Tr ! < < i 


( 2 . 21 ) 


and from the phasor diagram of Fig. 2.4 we see that. 


-ir l T c T R |sin 11 
jrf(f,t) Tan ^ j T . + | T T 

o' c R 


“! 
- I 


cos r)J 


= Tan 


-1 


r 1 r' 

L i.+ it r 


sin T) ^ 

7 ! 

I COS T| J 


( 2 . 22 ) 


Using Eq. (2.21) this expression can be approximated by; 

! T r 1 sin n 
~ 1 + T R cos ti 

~ ]T | sin rj (2.23) 

The phase angle v has been left as a parameter to allow the dis- 
cussion to proceed in terms of an arbitrary reference. In 
practice, the receiver phase reference may be locked to the 
received carrier using a phase locked loop. In addition the 
observable amplitude |H(f,t ) \ can be written as; 

|h (f , t) | = |tJ + 1 T c T r ! cos T| 


= |t !(i+ It | cos n) 

C K 


(2.24) 
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Fig. 2.4. Vector Sum of Random and Coherent 
Signals 
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Thus it is seen from Eqs . (2.23) and (2.24) that the observed 

phase and amplitude fluctuations of the tone response corre- 
spond directly to the components of the scattered signal in phase 
and in quadrature to the time invariant component T . Phase mea- 
surements will of necessity be relative to a reference, so that 
the actual observed phase will be v + ^ rather than $ alone. 
However v can easily be removed after averaging the phase pro- 
cesses . 

The discussion above is useful in providing insight as well 
as specific approaches to measurement of the channel parameters. 
For example the phase fluctuation properties are of interest in 
developing designs for pilot tone phase locked loops. Eq. (2.23) 
suggests that the phase of a narrowband signal can be described 
as a zero mean Gaussian random process, since it represents the 
quadrature component of a random vector. 

The approximations used here are of further interest in 
establishing the utility of certain measurement procedures as 
will be seen for example in Section 3.2.2. 
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2.4 Channel Measurement Philosophy 


The composite channel filter H(f,t) can be examined 
using either a time domain or frequency domain approach. In 
particular the time varying impulse response h(T,t), which is 
the Fourier transform of H(f,t), can be examined directly 
using periodic signals having narrow autocorrelation proper- 
ties. The period must be chosen smaller than the inverse 
Doppler rate which determines the time variation of h(T,t), 
and yet large enough to exceed the duration of the impulse 
response. To satisfy both of these constraints we require 
that the delay spread - Doppler spread product be less than 
unity (see Ref. [2.8]). This is easily satisfied for the 
millimeter wave channel. 

There are three prime candidate techniques which might 
be suitable for millimeter wave channel measurement. (see 
Refs. [2.9, 2.10]). 

(a) Sampling of H(f,t) using tones spaced in frequency. 

(b) Probing of the channel in the time domain using a 
periodic pseudo random phase modulated carrier. 

(c) Time domain probing using a sequence of short pulses. 

Note that in all cases the required bandwidth of the signal 
must be of the same order or larger than the effective band- 
width of H ( f , t ) . 

In the material which follows we shall assume that the 
complete impulse response g(t,§) or its transform H(f,t) is 
available. The former can be most easily obtained using coherent 
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processing of a sequence of pulse responses (see Chapter 4) 
while H ( f , t ) is available in sampled form from a multitone 
probing signal (see Chapter 3). 

2.4.1 Extraction of Coherent Filter Parameters 

The term T (f,t) contained in the expression for the 
R 

channel transfer function has been defined as a zero mean 
process. The obvious way of establishing estimates of T^(f) 
then, is to average the quantity H(f,t) over a long time 
interval (see Section 6 for a discussion of averaging times). 
Then we can proceed on the assumption that, 

<H(f,t)> ~ T (f) 

— c 

where <> denotes time average. In this section we discuss 

the characteristics of T (f) and their relationship to the 

c 

physical properties of the link. 

For a total path length x through the medium which has a 

frequency dependent refractive index n(w), the phase angle, ^ 

of T (f) is found from 
c 

$ = -kx (2.25) 

where 

k = “ n (uj) 
c 

n = refractive index 
c = velocity of light in free space 

For this discussion we will assume a lossless medium. 
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Expanding about the center frequency 


n(u>) ~ a Q + a 1 (o)-uu 0 ) + 


(2.26) 


Thus 


2 

jH*(uu) = - ^ x[a Q + a x (uj-U)q ) + a 2 (u)-uJ 0 ) ] 


(2.27) 


Upper and lower sidetones at therefore have phase 

parameters ; 


■ 


-<V"m )x . 2 


[a. + a oi + aw ] 
0 1 m 2 m 


(2.28) 


and 


-(<u -u>)x ? 

f& = [a - a.u) + a_w ] 

^2 c 0 lm 2m 


(2.29) 


Note that terms linear in frequency give rise to delay effects 
only. At the center frequency, uu=uj q , we have 


. 'V 

= — r a o 


(2.30) 


Then 


*1 = - ftWW + “’m'Wo 1 + Vm 1 <2 - 31) 


^2 - **0 - c C -% (a 0 +a lV + 'V'Wo’ - a 2 U, m ] (2 ' 32) 
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Thus it can be seen that the currently proposed ATS-F pro- 
cessor will give a phase measurement; 


~ 2 ^o 


-2x 

c 


2 

uu 

m 



+ Vo' 


(2.33) 


i.e., in the general case only terms in the phase expansion 
which are even functions with respect to the carrier are 
retained. To relate the parameters a Q , a^, a^ and the correspond 
ing frequency dependent attenuation characteristic to spreading 
of a pulse envelope, one must represent the pulse by its spectrum 
multiply by the equivalent channel filter (using the phase 
characteristic of eq. (2.27)) and transform back to the time 
domain. The details will not be carried out here, but we can 
note the following points; 

(a) The parameter a Q is related to delay effects only 
and is not generally of interest in characterizing 
the channel. 

(b) Parameters a^^ and a 2 are directly related to pulse 
distortion effects and they can be used to predict 
bandwidth limitations due to dispersion [2.1]. Thus 
the quantity provides useful experimental data 
for this channel characteristic. 

Although the role of amplitude variations 1 ( f ) | has not 
been fully discussed, it is clear that amplitude variations 
across the band also give rise to pulse distortions. However 
the required measurement 'technique is straightforward and needs 
no elaboration. 
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Note that the coherent channel function T^Cf) is time 
invariant provided that the bulk medium properties are fixed. 

Hence knowledge of T (f) allows equalization of the channel to 
counteract intersymbol interference effects which are likely 
to occur. 

2.4.2 Measurement of Scattered Signal Properties 

This section describes the rationale behind a measurement 
technique for the scattered signal properties. The technique 
is presented in detail in Section 6 below. The method relies 
on detection and processing of in-phase and quadrature components 
of the complex envelope of the received tones, rather than 
their amplitudes and phases. 

First, let us consider the statistical characterization 
of the amplitude and phase processes A(f,t) and \(f,t) (see 
Eq. (2-4)) in terms of quadrature component properties. Clearly 
A and y will not be Gaussian processes, so that a complete 
description would entail substantially more than the specification 
of correlation functions. However some useful relationships 
between in-phase-quadrature process correlation functions and 
phase— amplitude correlation functions can be established.* The 
results available are applicable only to wide sense stationary 
processes (i.e., satisfying Eqs. (2 . 21) and (2. 22) ) , and often only 
for symmetrical spectra (i.e., satisfying Eq. (2.23)). The con- 
clusions may be summarized as follows; 

*See for example the treatment of this subject given in Daven- 
port and Root, p. 158 or in Middleton, p. 396. These authors 
are primarily concerned with the statistics of a deterministic 
signal plus narrow band Gaussian noise after passage through 
various nonlinear devices such as detectors. The narrowband 
noise can be shown to have equal variance quadrature components 
and be wide sense stationary in character. 
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(a) When & =0 and the channel function H(f,t) is sta- 

ll 

tionary in time and frequency, we have a complete characteriza- 
tion in terms of 


K H (f l" f 2 


t l" t 2 ) 


- 2[ V f r f 2' 


W 


xy 1 2 


W 1 


(from Eq. (2 . 18) ) . 


(2.34) 


(b) The amplitude process A(f,t) has a correlation function 

R which is a complicated function of the magnitude squared of 
A 

Kjj , i • e . , 

|K I 2 = 4 (R 2 + R 2 ) (2.35) 

1 H‘ x xy 


(c) The phase process Q(f,t) similarly has a correla- 
tion function R which is a complicated function of both phase 

and amplitude of K . (Note the simplifications which arise 

H 

when the spectrum is symmetrical and R^ = 0.) 


(d) It is difficult to use the relationships provided 
in (b) and (c) to estimate R and R from R and R . For the 
millimeter wave channel it would be necessary to generalize 
these results to include three parameters R^, and R^ so 
that the two functions R and R would obviously not provide 
sufficient information even if closed form relationships could 
be found. 


The real motivation for pursuing this particular topic is 
related to the following question. Suppose that correlation 
functions for the ohase and amplitude were available.. Could 
these be used to determine the complete correlation properties 
of the process r(f,t) over a given frequency band? Furthermore 
it is important to know how the phase and amplitude correlation 
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functions relate to the communication performance, particularly 
whether or not they provide enough information to determine the 
bit error rate. 

Of course it is not necessary to restrict attention to 
the phase process correlation function. Any suitable statistic 
such as moments of the observed phase quantity, could be esti- 
mated. For example if the measurement was 


© = 9(f 1 ,t 1 ) - 0(f 2 ,t 2 ) 


(2.36) 


/v 2 

then E(9 } can be related (with some difficulty) to the para- 
meters R , R and R . Note however that a phase measurement 
x y xy 

of the form. 


9* = 9(f 1 ,t) + 9 (f 2 , t) - 29 (f Q , t) 


A 


e. 


+ e. 


29 


o 


is even more difficult to handle because it must be related 

to R , R and R at three different frequency pairs, 
x y xy 

In summary we can say that the measurement of channel 
parameters using directly observed in-phase and quadrature 
components appears to eliminate many of the difficulties 
inherent in the task of interpreting phase and amplitude 
measurements. We should note that the implementation of a 
phase difference measurement of the form 

9 = 9(f ,t x ) - 9(f 2 ,t 2 ) 
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requires the use of delay lines for on-line processing, or 
alternatively ©(f^t^ and 0(f 2 ,t 2 ) must be measured individually 
with respect to a stable reference and recorded for off-line 
processing. 

In a recent paper, Mondre [2.7], established some useful 
bounds on the behavior of a complex covariance function in terms 
of the covariance function for the transfer function amplitude 
|H(f,t) | . Although a Ricean channel model was assumed, the 
analysis yields useful insight. Briefly, the result was 
obtained that the covariance function for the amplitude process 
was always bounded by the magnitude of the complete complex 
correlation function. Furthermore as the ratio of coherent 
non- fading power to average random power increases, the bound 
tends to become tighter, and in fact leads to true equality in 
the case of symmetrical spectral properties (i.e., = 0; 

see Section 2.3.2). This latter result is easily explained 
in terms of the material given in Section 2.3.3 where the 
case of low power scattered signals was discussed. When the 
cross — covar iance function R^^ i- s zero we can express the random 
signal in terms of a pair of independent equal variance in-phase 
and quadrature components for any arbitrary reference. Thus 
when the coherent signal is relatively large, amplitude measure- 
ments yield an approximation to the sum of T^ and the co-phasal 
T component. Hence the amplitude covariance is simply the 
correlation function for the in-phase random component which 
provides a complete characterization for R = 0. When this 
latter condition is not satisfied, as will be the case for 
reasons given previously, it is obvious that amplitude informa- 
tion will not be sufficient even for high coherent/random signal 
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power ratio, although the bound is still valid. This is even 
more true for a non-Ricean channel which will frequently occur 
in millimeter wave systems, (see Section 2.3.2 and Appendix A). 
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Section 3 


CHANNEL MEASUREMENT USING COHERENT TONES 
The basic idea behind this scheme is the sampling of 
H ( f , t ) using a set of spaced coherent tones. If the spacing 
is sufficiently small, the complete function H(f,t) can be 
reconstructed and processed to obtain its correlation proper 
ties. Generally speaking it is only necessary to correlate 
individual tone responses thus obtaining an estimate of the 
sampled correlation function. It will be seen that a major 
difficulty in this method is in maintaining phase stability 
for the tones. In fact the true channel response will be 

masked by, 

(a) Instabilities in the transmitter oscillator. 

(b) Movement of the spacecraft giving rise to Doppler 
shifts in the received signal. 

(c) Instabilities in the receiver oscillator output. 

One could argue that (a) and (b) will be part of any 
ultimate communications system anyway and should be measured 
as part of the channel. However the choice of oscillator for 
an initial experiment may be somewhat restricted and not en- 
tirely representative so that it would appear better to try 
and separate the effects of oscillator and propagation medium. 

Alternatively it can be shown that for certain configura 
tions both (a) and (b) give rise to a phase error which is 
linear with frequency. Careful processing at the receiver will 
eliminate this variable linear term, but also unfortunately 
introduce certain errors in the estimation of the medium in- 
duced linear phase component. A linear phase shift with fre- 
quency corresponds to an instantaneous delay in the time domaxn. 
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and constitutes quite an important aspect of the millimeter 
wave channel. It would be desirable to know the typical or 
RMS delay fluctuation and its rate of variation. 

Thirdly, we can consider the use of averaging times that 
are short compared with the oscillator instabilities and satel 
lite induced phase drifts [ 3 . 1 1 . For the millimeter wave chan 
nel, the probing bandwidths of interest do not allow this kind 
of approach, as will be seen subsequently. 

The individual channel probing schemes will now be ex- 
amined in more detail beginning with the approach that looks 
the most promising to date, which is based on the multitone 
signal. 

3.1 Form of the Transmitted Signal 

The basic transmitted waveform consists of several co- 
herently phased tones. These will be denoted by, 

S (t) = a sin (oj t + 8 ) n=0 , +1 , . . . , +N (3.1) 

n n n n 

After phase and amplitude distortion by the channel, the re- 
ceived tones are processed so as to yield estimates of the 
desired channel properties. 

The tones must be generated from a single low frequency 
reference oscillator to ensure coherency. The reference is 
multiplied up to the desired RF band using a synthesizer. 

In addition the reference can be used to phase modulate a 
carrier so that a constant amplitude signal is made available 
to the final stages of the transmitter. 
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The following analysis pertains to the latter technique 
which is to be used to generate the ATS— F satellite millimeter 
wave probing signal. 

Consider for example the 20 GHz signal. Starting with 
a reference oscillator of frequency f ^5 MHz, a carrier is 
obtained by multiplying up to 10 GHz i.e., using a multiplica- 
tion of K = 2000 we have, 

. c (t) = sin (Kid t + K0 ) (3.2) 

s s 

where id and 0 are the frequency and phase of the reference 
s s 

and, 

Kuo 

— — = 10 GHz. 

2 TT 

The carrier is phase modulated in such a way that 65% 
of the output is in the multitone spectrum. The phase modulator 
is driven by the output of a multitone processor which is in 
turn driven by a harmonic of the oscillator. The multitone 
processor consists of an amplifier and matching network and. 
provides a sinusoidal drive to the modulator, i.e., 

m(t) = p.sin (kiD t + k0 ), k = 36 (3.3) 

1 s s 

Hence the transmitter output after x2 multiplication is; 

s ( t ) = sin ( 2K (iDgt + 0 g ) + 0 sin (ku) g t + k3 s )) 

0 = 20 1 (3.4) 

Figure 3.1 shows the basic transmitter structure. 
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Fig. 3.1 Transmitter 
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In terms of a carrier frequency uu and modulation frequency 


u>; 

m 


s(t) = sin (a' Q t + 0 Q + 8 sin ( uo^t + 9^) ) 
which results in sidetone expansion of the form; 

i(uu t + 0 +8, sin ( u) t + 0 )) 

s ( t ) = Im{e u u i 


( 3 . 5 ) 


where 


Thus 


j ( uu t + 0 ) jniu t 

= Im{e 0 0 E a n e m } 

n 


a = 


1 _ 


j ( 8 sin ( uu t + 9 ) 
J m m 

- mu it) 
m 

2tt . 


e 


1 

J n % 

p 11 + j ( P sin x - 

nx) dx 

2tt 

© 

-TT + 0 

m 


jn0 m 




dt 


J (B) 

n 


s(t) = F J n ( P ) sin (u) Q t + 9 0 + n(iu m t + 0 m ) ) 
n = -co 


( 3 . 6 ) 
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Thus if we represent the composite signal s(t) as the sum of 
tones, 

s (t) = E a s (t) = E a sin(uu t + 8 ) (3.7) 

n n n n n 

we see that. 


a = J (p) (3.8) 

n n 

(o = + nuj = ( 2K + kn)io (3.9) 

n 0 m s 

0 = (2K + kn) 6 (3.10) 

n s 

The index (3^ can be chosen so that the desired coefficients 
{a n } are similar in magnitude and negligible outside the de- 
sired band. Differential amplitude effects can easily be 
accounted for provided that the system is correctly calibrated 

We now proceed with a discussion of the tone instabili- 
ties induced by oscillator variability and satellite motion. 

3.1.1 Phase Instabilities 

The following areas are of importance in this discussion 

(a) The effect of satellite oscillator instabilities. 

(b) Phase shifts due to satellite movement. 

(c) Receiver local oscillator instabilities and their 
relationship to derived references. 

It will be demonstrated that (a) and (b) induce a phase devi- 
ation which is linear across the band, while the influence of 
the receiver local oscillator depends on the precise nature 
of the processing employed. 
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Let us assume that the satellite reference phase is 


u) t + 0 (t) 
s s 


where 9 (t) is a random time varying term which accounts for 
s 

the oscillator drift. From Eq. 3.10 it is seen that the in- 
dividual tones are affected also. 


0 (t) = ( 2K + kn) 0 (t) 
n s 


(3.11) 


Alternatively in terms of the carrier phase 0 , 


- l 1 + S 


(■ 


0 (t) 

o 


(3.12) 


Thus for a phase reference (or equivalently frequency) pertur- 
bation, we experience a corresponding linear phase perturbation 
across the band. 

Similarly movement of the satellite causes phase errors. 
If the total transmitter-receiver path length is x, the phase 
shift at the carrier corresponding to this free space path is. 


x 


U) 


t|( = kx 

o 


k = 


0 


Hence for a variable path length, 


0) 

* n U) = f x(t) 


-£ + £) 


qjx(t) 


■ C 1 + i) *c (t > 


(3.13) 
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i.e. a linear phase variation with frequency. Because the 

form of i (t) is identical with the transmitted phase proces- 
n 

ses 9 n (t) °f Eq. (3.12), we will not distinguish between the 
two effects in future, and simply lump f variations in with 
those caused by oscillator instability. Equations (3.12) and 
(3.13) indicate that provided the carrier phase 0 q ( t) is 
tracked, the incremental phase effects at the sidetones can 
be accounted for. In other words the carrier phase contains 
complete phase information for the total spectrum. An incre- 
mental phase correction procedure of this kind will be dis- 
cussed in later sections. For the moment we merely add the 
warning that transmitted signals must also pass through a dis- 
persive medium which does not have linear phase frequency 
properties, and hence any correction scheme based on received 
carrier phase will in fact introduce phase errors. 

The discussion of item (c) , instabilities in the receiver 
reference oscillator) will be postponed until the next section 
which deals with the basic receiver structure. 

3.2 Basic Receiver Structure 

The total bandwidth of the transmitted spectrum must be 
greater than 1 GHz, thus dictating the use of individual tone 
processors i.e. each of the tones must be supplied with a 
reference such that after mixing, individual tones are con- 
tained in separate but similar channels at a second IF of 
about 10 MHz. At this frequency they can easily be manipu- 
lated and processed to yield the desired phase and amplitude 
information. Figure 3.2 shows the general receiver structure 
of interest. 
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Oscillator 


Fig. 3. 2. General Receiver Structure 
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The signals which would exist at the receiver for a 
free space path can be written as, 

s (t) = sin (ou t + 0 ) n=0, ±1***±4 (3.14) 

n n n 

Note that unity amplitude has been assumed to simplify the 
notation. These tones are received subject to phase and am- 
plitude distortions due to randomness of the propagation 
medium; 

r (t) = a sin (id t + 0 + ) (3.15) 

n n n n n 

where <6 and a are the channel induced phase and amplitude 
n n 

fluctuations (relative to the free space response) . 

The received signals are first mixed down to a conveni- 
ent IF band. 

The resulting IF signals can be represented as, 

r (t) = a sin (u> t + B + + 0 ) (3.16) 

in n In n n I 

where uu ^ Q is the center IF frequency and 0^ is the phase of 

the IF reference. Typically the IF frequency will need to be 

a few GHz to accommodate the total spectrum bandwidth. For 

the proposed ATS-F system it is 1.8 GHz. The radian frequency of 

the n-th tone at the first IF is in. . 

In 

To reduce each tone to the second (10 MHz) IF, two pos- 
sibilities exist. 

(a) Use i icividual references offset by m „ below each 

IF 

tone frequency. (See Fig .3 .3a) . here denotes 

the second IF frequency. 
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Fig. 3.3a. Phase Difference References 
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(b) Use a reference between the upper and lower side- 
tones of interest. (See Fig. 3.3b) with subsequent 
translation to tu . 

Figures 3.3a and b show these two alternatives in terms of 
a pair of tones n = ±1. 

For the processing outlined in (a) the signals repre- 
sented by Eq. (3.16) become ; 

r n (t) = a n sin (u^t + ^ \ + e IF > (3.17) 

0 IF is the phase introduced by the references. At this point 
it should be noted that the references will be synthesized 
from a stable reference oscillator with frequency and phase 
denoted by u)* and 0* respectively. Hence 0 will depend on 
n, although the precise nature of this dependence will vary 
with the method used to derive the references. For example 
the references could be obtained directly from the receiver 
oscillator, 

= (2K+kn)ro* 

e R„ = (2K+kn) K - 

where to. ~ uu and N is chosen so that the desired second IF 

* s 

results. 

IF 

Alternatively using a phase locked loop to track the 
carrier at the 1st IF (m ) » we could derive references of 
the form; 

W Rn = W I0 + kn,JJ * " 

V = e io + kn9 * - N0 * 

Again N is chosen so that a second IF at uu results. 
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For this latter situation it can be shown that the re- 
sulting phase difference between tone n and tone — n can be 
expressed as; 



8 —0 + <6 —d> — 2nk9 * 

n -n n -n 


= f6 (t) (t) + 2nk ( 0 (t) -9*(t)) 

n -n s 


(3.18) 


The difference 0 -0 was obtained using Eq. (3.11) . Note 

n -n 

that; 

(i) The usefulness of such a measurement over a total 
interval (0,T), depends on the stability of both 
transmitter and receiver oscillators. 

(ii) 8 and 0. will always be unknown parameters; at 
s * 

best it can be hoped that they will be constant 
over the measurement interval, thus introducing 
only a fixed unknown path delay. 

(iii) Similar results are obtained for phase differences 
between say tone n and the carrier with an obvious 
alteration of the factor 2nk to nk. 

(iv) The phase 8^ introduced by the first mixer refer- 
ence is removed by the phase differencing process, 
since it is not a function of n. 
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It is immediately apparent from Eq. (3.18), that the diffi- 
culties suggested in (i) will be alleviated by locking 9* to 

9 . It is not possible to do this exactly since neither 0 
s s 

nor 8 can ever be observed directly, but estimates of 0 ob- 
0 s 

tained by tracking the center frequency phase could, in prin- 
ciple, be obtained. This is the substance of material presented 
in Sections 3.2.1 and 3.2.2 where we describe an effective phase 
tracking scheme, referred, to as the incremental phase correc- 
tion approach. 

Let us now return to item (b) , which suggests an alter- 
native processing scheme using a central reference. 

Assume that a reference locked to the carrier at IF, is 
available to the second mixer. i.e. the reference is repre- 
sented by, 

sin ( UJ io t + 6 0 + ^0 + V (3.19) 

After mixing with tone n we choose the frequency difference 
signal. 


r = a sin i (uu 
n n L 10 


- “in' 1 ±(9 n- 9 0 + W] (3 - 20) 


The plus sign applies to tones above the IF center frequency 
(«I 0 # while the negative sign applies to those below. Subse- 
quent translation down to uu introduces an arbitrary phase 
constant which need not be considered, provided we are only 
interested in the phase relationship between r^ and r 
Their phase difference will be; 
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= 0 + 0 + (6 + ^ - 2 (9 + rf n ) 

n n -n n -n 0 0 

From Eq. (3.12) it is seen that 


0+0 
n -n 



Hence 

A^ n = + *L n " 2 ^ 0 (3 * 21) 

Thus it appears that is not influenced by either satellite 

or receiver oscillator instabilities and it removes the effect 
of total path delay. Unfortunately it also eliminates any 
linear phase vs frequency trends in the channel phase param- 
eter f6. In effect this measurement removes the instantaneous 
linear phase variation across the band. More generally it can 
be shown to remove all terms with odd exponents in a Taylor's 
series expansion of the phase process. To conclude we note 
that there is no difficulty in making amplitude measurements 
once the signals are available at the second IF. 

In phase and quadrature components can easily be extracted 
from each of the IF tones. The resulting in-phase and quadra- 
ture processes will be directly related to r^ (t) of Eq. (3.17) 
mixed with the tracked carrier at the second IF. The demodu- 
lator will have a reference signal; 


sin (<W-j.pt + 


e n + & 


0 


+ !0 + 


V 


(3.22) 


and after mixing tone n will yield in-phase and quadrature 
components X^lt) and Y (t) respectively; 
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X = a cos 

(6 - 

9 „ 

+ jrf 


- nk6*) 

n n 

n 

0 

n 

0 


Y = a sin 

(0 - 

0 A 

+ 


- nk0 ) 

n n 

n 

0 

n 

0 

* 


or using Eq. (3.11) we get. 


X (t) = a (t) cos [Vjt) -^ n (t) +nk(0 (t) -9*(t)) 
n n L n u ^ -i 


Y (t) = a (t) sin 
n n 


^ (t) -*L(t) + nk ( 9 ( t ) - 0* 
n u s 


(tnl 

—j 

(3.23) 


Unless the reference oscillators are highly stable, the in- 

phase and quadrature components will be inadvertently phase 

modulated. For example consider an oscillator stability of 1 
g 

part in 10 averaged over 10 msec. The expected phase error 

over this interval for n=4, k=36 and u; =5 MHz can be computed 

s 

from. 


e = stability x knu^T 

-9 6-3 

=10 7200 x 10 10 

-3 

= 7.2 x 10 cycles. 

While this is not excessive for the short time interval and 
high oscillator stability specified, it indicates that serious 
difficulties will arise for realistic averaging times on the 
order of 10-100 seconds. Furthermore it demonstrates that an 
oscillator of high quality is mandatory. The consequences 
of processing over time intervals using a poor quality oscil- 
lator may be seen by representing the introduced phase error 
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in terms of a random variable uniformly distributed over 2tt. 
Averaging in an ensemble sense or time averaging as in the 
actual measurement, the result will approach zero. 


3.2.1 Carrier Tracking and Incremental Phase Correction Schemes 

In this section we explore whether a receiver structure is 
capable of correcting for the effects of satellite mostion and 
system oscillator instabilities. Eq. (3.18) and (3.23) indicate 
explicitly the consequences of transmitter and receiver oscillator 
phase fluctuations from the point of view of both direct phase 
comparison measurements and quadrature and in-phase measure- 
ments. The satellite motion introduces phase variations simi- 
lar in effect to transmitter instabilities so that they can be 

included in 0 (t) . The basic idea of this section is to extend 
s 

the concept of the phase lock loop in its role of tracking the 
carrier phase. In the absence of any scattering effects the 
channel can be represented as a time invariant filter T c (f), 

(See Section 2.2) coupled with a variable delay . Figures 3.4 and 3.5 
illustrate this simplified channel model. The delay accounts 
for all linear phase vs frequency characteristics introduced 
by oscillator and satellite variations. Since the dispersive 
medium is physically separated from these effects, it can be 
dealt with independently in the system model. 

The phase locked loop essentially establishes the phase 
at the carrier frequency. An incremental phase correction based 
on Eq. (3.12) can then be implemented, i.e. using the relation- 
ship; 

e n (t) = ( X + t) 0 O (t > 
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If 0 (t) is available at the reference synthesizer, correc- 
0 

tions can be made. Figure 3.6 shows the structure of a phase 
locked loop which would be compatible with the currently pro- 
posed ATS-F receiver. 

It is important to note that the first mixer reference 
must be generated from the reference u)*. The use of an inde- 
pendent oscillator at the first mixer introduces phase fluctua- 
tions, so that the incremental correction scheme would respond 
incorrectly. 

While the configuration shown in Fig. 3.6 is desirable 
from a practical point of view, it is instructive to first 
analyze the operation of the incremental correction scheme in 
terms of a simplified structure. Let us assume, as shown in 
Fig. 3.7, that the loop encompasses the first mixer so that 
all of the loop frequency translations can be combined using 
a single mixing operation. Under locked conditions both in- 
puts to the phase detector have the same phase, and at the 
mixer, (tu^ and 0^ relate to the VCO-multiplier output) 


%*■ + 0 o " + V = UJ * t + 9 * 


(3.24) 


By providing the offset between the mixer and the VCO we 

are able to choose m „ ~ (ju . Hence the VCO output has phase, 

vO s 


“vO * 1 + 


vO 


1 

2K 


(u^t + 


6 1 + 


w*t + 


6*) 


_1_ 

2K 




V 


(3.25) 


The latter relationship follows from Eq. (3.24). Then from 
Eqs. (3.9) and (3.10) we see that 
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Fig. 3.7 Simplified Phase Locked Loop Structure 









(3.26) 


uu _t + 

vO 


vO 


0) t + 0 
s s 


and the references can be generated using the relationships. 


(3.27) 

(3.27) are employed, the 
individual tones are mixed down to baseband, which is desirable 
if in-phase and quadrature components are to be measured. Al- 
ternatively the references can be offset by some nominal IF 
frequency uj . 

The simple minded approach depicted by Fig. 3„ 7 can be 
expanded to avoid various impracticalities but the basic pro- 
cedure remains the same. For example Fig. 3.8 indicates the 
way in which the signal can mixed down to an IF for easy pro- 
cessing, without destroying the phase tracking properties of 
the receiver. In addition the phase detection is carried out 
after a further reduction in frequency to 100 kHz. The choice 
of and K is governed only by practical considerations sub- 
ject to the condition that; 


“Rn ’ (2K + nk) “Vo 


0_ = (2K + nk) 0 n 

Rn vO 


When the references specified by Eq. 


K + K = 2K 

J- « 

where 

2K = 2000 in the proposed experiment. 

In general .it would not be desirable to have the first 
IF lower that about 2 GHz because of the wide signal bandwidth, 
hence 
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Fig. 3.8 Phase Locked Loop Design 














K ~ 200 

1 t 

K 2 1800 

The extraction of in-phase and quadrature tone components must 
be carried out at a second IF band to avoid complications as- 
sociated with the 90° phase shifting of the references. Figure 
3.9 illustrates one possible method of implementing this mea- 
surement. 

At this point it would seem that the procedures involved 
with an incremental phase correction, are not too difficult to 
realize. However we have not yet broached the question of the 
utility of such measurements when the channel exhibits arbitrary 
phase characteristics due to the presence of a scattered signal 
component. This topic will be investigated in the next section. 

3.2.2 Analysis of the Incremental Phase Correction Receiver 

For the simplified channel model shown in Fig. 3.5 which 
involved free space propagation, the incremental phase correc- 
tion procedure completely removes time variations in the free 
space path delay. Furthermore when the actual medium can be 
represented by an equivalent (possibly time varying) refrac- 
tive index, the same conclusions follow. In these cases all 
of the tone phases are coherently related to the carrier phase 
8g(t). However when the signal received is the sum of scattered 
plus direct energy, no such relationship exists. Thus a linear 
phase correction scheme incorrectly interprets changes in the 
mid-band phase by assuming that the linear portion of the com- 
posite phase response is in fact proportional to 8^ + 

Fig. 3.10depicts the various phase processes involved, while 
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Fig. 3.11 indicates the complete channel model. (See Section 
2) . Note that the correction scheme does not actually elimi- 
nate the instantaneous linear phase component as one might sup- 
pose. This would only occur if + 0^ was deterministically 
related to the linear phase term, which is not necessarily 

true when the filters T_(f,t) or H(f,t) are involved. 

R 

To best illustrate the performance of the phase correction 
technique, we first examine a continuous model for the phase 
process. The instantaneous received phase expressed as a func- 
tion of both time and frequency is, 

y(t,t) = *$(f,t) + 0 (f , t) (3.28) 

where 9(f,t) represents the time varying phase shift corres- 
ponding to the model shown in Fig. 3.5 (i.e. filter T^ plus 
a variable delay). 0(f,t) will generally deviate from a linear 
function of frequency because of the frequency dependent nature 
of the medium equivalent refractive index. (See Appendix B ) . 
However its time variation originates almost entirely from the 
variable delay. In any case 9(f,t) can be expressed as, 

0 (f , t) = 9 (t) f"l + 7“J + £ ( f ) (3.29) 

L j-qJ 

where 0 Q (t) = phase at carrier 

Af = f - f Q , deviation from carrier 
e = remainder in phase expansion 

^(f,t) on the other hand is the phase perturbation introduced 

by the scattered signal contributions (See Fig. 3,5). It can 

be expressed in terms of the real and imaginary components of 

the filter T as, 

R 
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1 r Im{T R (f,t)} ^ 

= tan ‘ Ll + Re{T R (£.t)U 


(3.30) 


However one should note that 


T r| « l 


(See Appendix B) ; thus we have the approximation, 

Im{T R (f,t) } 

~ 1 + Re [T r ( f , t ) } 


~ Im{T_(f,t)} (3.31) 

R. 

These approximations can be justified because the strength of 

the scattered signal is expected to be at least 20 dB below 

the attenuated direct component assuming narrow beamwidth earth 

station antennas. Consequently we can conclude' that jzf(f,t) can 

- 2 

be quite small, perhaps less than 10 radians on the average 
(for a 20 dB unscattered/scattered power ratio) . 

Consider an expansion of about the carrier frequency; 

rf(f.t) ~ jrf (t) + a x (t)Af + 6(f,t) (3.32) 

« 

where Af = f - f^ 

Then the total observed phase is, 

.y(f,t) = jz$(f,tj + 8(f,t) = ^ Q (t) + a 1 (t)Af 

+ 8 Q (t) [l + jr\ + e(f) + S(f.t) (3.33) 
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The phase locked loop (PLL) tracks the carrier phase 


V° = 8 0 (t) + fl 0 M 

In discussing phase processes such as this, we are implicitly 
dealing with the phase relative to some fictional reference. 

By representing oscillator and satellite instability effects 
in terms of an equivalent delay as in Fig. 3.5, we can assume 
this reference to be some ideal, perfectly stable oscillator 
on board the satellite. Note that y (t) is the true carrier 
phase observed at the receiver; S^Ct) includes the free space 
phase shift corresponding to a large number of wavelengths 
separating transmitter and receiver. The receiver phase locked 
loop supplies an estimate Y Q (t) , or more exactly, a reference 
with phase V Q (t) , i.e.. 


A cos (uu Q (t) + y Q (t) ) 


Y 0 (t) = e 0 (t) + *yt) 


(3.34) 

(3.35) 


Let us now consider the phase corrected demodulator re- 
ferences which, as discussed in Section 3.2.1 are generated from 
the PLL output. Although this procedure is to be implemented 
in the receiver IF sections, the principles are identical re- 
gardless of the frequency band, and for simplicity we carry out 
the details in terms of the RF signals. 

The references are generated with a phase of the form; 



e R (f,t) 


(3.36) 



and this is essentially subtracted from Y(f*t). Ideally 
Y (t) = Y (t) in which case; 

Y(f.t) = Y (f , t) - 0 R (f , t) 

= f6 Q ( t) + a^ (t) Af + e Q (t) (l + I 1 ) + e(f) + 6(f,t) 

- i 1 -ff)( e o (t > + <3 - 37) 

Therefore 

Y (f # t) = (a 1 (t) - -j— ) a f + 6(f,t) + e(f) (3.38) 

In Eq. (3.38) we see that aside from the possibility of 
errors in the estimate of Y Q (t) , there is a time varying phase 
error present. i.e. the desired form of Y(f,t) is, (from Eq. 
(3.32) and (3.29)) 

Y ( f , t) = ^(f,t) - ^ Q (t) + e(f) 

= a x (t)Af + 6 (f , t) + e (f ) (3.39) 

while the measurement phase given by Eq. (3.36) contains an 
additional term, 

^o (t) 

^ (f,t) = — Af (3.40) 

£ o 

The severity of this error depends on the nature of the 
channel characteristics and the information we are trying to 
extract from the received data. In general it can be seen 
that the factor 
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iC 0.036 = 


720 MHz 
20 GHz 


ensures that the error is small whenever 0 Q (t) and the phase 
process of interest 0(f,t), are not strongly correlated. At 
the other extreme it is clear that when T(f,t) corresponds to 
a simple time varying delay, the error component reduces the 
phase 0(f,t) to zero. 


For situations between these two extremes we need to ex- 
amine some specific measurements. Consider the true in-phase 
and quadrature components of the scattered signal relative to 
the so called specular component represented by T c » (See 
Fig. 3.11). These quantities will be referred to as 


X = I t c t r' cos 11 

and 

y = |T C T R | sin T| (3.41) 


respectively. Following the spirit of previous approximations 
we note that, 


| H ( f , t) | ^ | T c (f ) | + x(f,t) (3.42) 

and 

0(f,t) ~ y(f,t) (3.43) 


However the measured quadrature component will be in error by 
an amount. 


y_ <f » t) ~ 0 (f , t) 
e — e 


y(( 0 .t) 
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(3.44) 


\ 



Figure 3,11 Vector Sura of Scattered and Coherent 
Signal Showing Phase Correction Error 
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Since the scattered signal is Gaussian in nature and zero 
mean by definition (See Section 2.2), y(f,t) is characterized 
by the function, 

V f l- V t l ' V = y(£ l- t l ) y(f 2' t 2 ) 

If we were to try and measure this quantity by averaging the 
product of two tone responses, we would obtain an estimate 
which approached, 

R y ( f l, f 2 ; t 2 ) = y(f 1 ,t 1 ) y (f 2 , t 2 ) (3.45) 

where 

y ( f i» t 1 ) = y(f 1 > t x ) - y e (f 1 , t ± ) 

etc. Thus 


ft (f 

y 


i' 



ty(f r t i ) - if y(£ o'V K y(f 2 - t 2 >- if y'W 3 = 


R y (f l' f 2 :t l t 2 > + (f Q ) R y (f O' f O :t l t 2 ) 





)+R 

y 


(f i' f o 



(3.46) 


This relationship clearly indicates the way errors are intro- 
duced in a correlation function measurement. To be more spe- 
cific let 
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and 


R y (f l' f 2 :t l _t 2 ) R y* f l f 2’ t l” t 2 ) 


f i - £ o - af 


= + ^ ' so t ^ iat - ^2 = 


At = 



Then 


R (2Af , At) 

y • 

R (2Af, At) 

y 


(^)r (0, At) - R (Af, At) 

' r n ' y r n y 


R (2Af. At) 

(3.47) 




The fractional error contained in the brackets is insignificant 

if. 


and 


R y (2Af, At) 
R y (0, At) 



R y (2Af, At) 
R (Af, At) 



(3.48) 


(3.49) 


It can be seen that because A f/f is at most 0.04 ( for the ATS-F 
system, these conditions do not pose serious problems. Note that for 
increasing Af we would normally expect Ry (2Af,Af) to approach zero. 

In addition the error parameter Af/f^ becomes larger, so it is 
clear that the fractional error will be most severe in the 
tails of R , i.e. for large Af. Similar conclusions can be 

y 

drawn for estimates of the in-phase and cross correlation func- 
tions R and R , which were discussed in Chapter 2. 
x xy 
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This concludes discussion of the incremental phase correct- 
ing receiver. Subject to the limitations imposed by Eqs. (3.48) 
and (3.49), the method appears feasible. Further study is how- 
ever required to ascertain whether available receiver compo- 
nents will support the ideal assumptions implicitly made in 
previous discussions. In particular it must be verified that 
the receiver does not superimpose any delay distortion on the 
signal . 

From this and previous sections we can intuitively see 
that if the true channel correlation functions are to be mea- 
sured, some form of phase correction to eliminate oscillator 
instabilities is a mandatory requirement. The technique out- 
lined here appears to provide a valid approach only if the 
actual channel- induced phase fluctuations ^(f,t) are quite 
small, and this in turn is due to the presence of a large spe- 
cular signal. If the phase variations due to the scattered 
component were in fact larger, it might not be possible to 
achieve such good separation of the channel from terminal in- 
stabilities. An exception occurs if the spectral characteristics 
of each are such that judicious choice of phase lock loop filter 
parameters allows separation on a frequency domain basis. This 
possibility is discussed next. 
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3.3 Application of Phase Locked Loops for Incremental 

Phase Correction 

It may be possible under certain conditions to separate 
the phase effects due to satellite motion and oscillator drift 
from the phase effects due to scattering in the channel. This 
is possible, for example, when the scattering fluctuations are 
fast and the phase fluctuations due to satellite motion and 
oscillator drift are slow. In this case a phase locked loop 
can be employed to separate the two phase erfects. 

Specifically, assume that the instantaneous received phase 
is defined as shown in Eq. (3.28), repeated here for convenience. 

Y ( f , t) = jrf (f , t) + 0 ( f , t ) (3.50) 

Where 8(f,t) includes the phase shifts due to oscillator in- 
stability and satellite motion and also includes the phase 
characteristic of the time invariant portion of the channel's 
transfer function. f,t) represents the phase perturbation in- 

duced by the scatter signal contributions at the receiver. 

Now, if the bandwidth of 0(f,t) is much less than the band- 
width of jzS(f,t) then a phase locked tracking loop can be used 
to separate the two effects at a particular frequency. 

A simplified block diagram for this purpose is shown in 
Fig. 3.12. The transfer function of the phase locked loop, from the 
input to the output at point A, H(s)» is typically a simple low- 
pass filter of the form given in Eq. (3.51). 
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H (s) 


(3.51) 


k (ts + 1) 

2 

s + kTS + k 

That is the transfer function is unity from DC to roughly k 
rad/s and then the transfer function rolls off at 6 dB per 
octave. The transfer function from the input to the error 
point, point B in Fig. 3.12 is given by l-H(s). Thus from the 
input to the error point, the loop highpass filters the phase 
data . 

In order to determine the applicability of the phase 
locked loop technique, the dynamics of the phase variations 
due to scattering and the phase variations due to oscillator 
stability and satellite motion must be defined. Appendix A 
shows that phase variations due to scattering create a 
doppler spread process whose bandwidth is on the order of 1 Hz, 
maximum. It is not expected that the Doppler spread will ex- 
ceed this value by very much. However it is possible that 
considerably smaller Doppler spreads will be encountered. The 
large Doppler spread cases are the ones which are of more in- 
terest because they involve heavy weather conditions. The dy- 
namics of the phase variations due to oscillator stability and 
satellite motion are not usually defined in terms of a band- 
width. Instead they are usually specified in terms of power 
series expansion coefficients. For example, in the case of an 
oscillator the predominant source of instability may be defined 
as an aging rate. Of course, for satellites the dynamics are 
defined in terms of maximum velocity and maximum acceleration 
along the line of site path. It is reasonable to assume that 
phase variations due to satellite motion will predominate over 
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those due to oscillator instability. Therefore let us investi- 
gate the bandwidths required by a PLL to accomodate a given 
amount of satellite dynamics. 

Three system configurations will be considered. First a 
digitally compensated second order loop will be presented. 

Then a digitally compensated third order loop will be discussed. 
Finally the use of rate aiding with either of these loops will 
be described. 

The loop compensation for a second order loop is 


F 2 (s) 


TS + 1 

S 


(3.52) 


The overall closed loop transfer function with such compensa- 
tion is of the form 


H 2 (s) 


k(TS + 1) 


2 

s 


+ kTS + k 


(3.53) 


Where k is the overall loop gain and t is a coefficient set by 
the specific digital design. For the loop to be critically 
damped, 

kT^ =2 (3.54) 


and the 3 dB bandwidth of the closed loop response of Eq. (3.53) 
is 

BW = ^ Hz (3.55) 

3aB 2n 
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The dynamic capabilities of tracking loops are typically 
expressed in terms of steady state errors due to position, 
velocity, and acceleration inputs. The steady state error of. 
the second order loop just defined is zero for both position 
and velocity inputs. However the steady state error due to 
an acceleration input is not zero. In particular 



oj 

(2l,BW 3dB> 2 


( 3 . 56 ) 


where is the phase error at the phase detector and a) is the 
doppler rate caused by satellite motion (and, perhaps, oscil- 
lator drift) . 

It is clear that the 3 dB bandwidth of the loop cannot be 
made indefinitely small. The lower bound in this case is set 
by the doppler rate caused by satellite motion. In particular 
the phase error at the phase detector must remain small so that 
the loop remains in its linear region. Typically an upper 
bound on the loop steady state phase error is set at 0.2 rad. 
From Eq. (3.56) it is seen that this constraint creates a lower 
bound on the 3 dB loop bandwidth given by 


BW 


3dB 




(3.57) 


In terms of satellite acceleration, a, the lower bound on 3 dB 
bandwidth can be recast as 


BW 


3dB 


> 



a f 
0 . 4nc 


(3.58) 
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where f is the carrier frequency to be tracked (in Hz) , a# is 
the radial acceleration of the satellite as seen at the re- 
ceiver (in ft/sec 2 ) and c is the speed of light(in ft/sec). 

At 20 GHz the maximum tolerable acceleration which per- 
mits a 3 dB bandwidth of 1 Hz using a second order loop is 
0.0628 ft/sec^. 

A third order loop provides a significant increase in dy- 
namic performance for a given bandwidth. Thus higher vehicle 
dynamics can be tolerated. The loop compensation for a third 
order loop is of the form 

F 3 (s) = ( JtiL T- 1 ) 2 (3 - 59) 


and the closed loop response is given by 


H 3 (s) = 


k ( TS + 1) 


+ kT 2 s 2 + 2krs + k 


(3.60) 


The equivalent of critical damping is achieved with the third 
order loop by setting 

kT^ ft* 4 (3.61) 

The loop 3 dB bandwidth under the constraint given in Eq. (3.61) 
is approximately 


BW 


3dB 


3.2 

5t 


(3.62) 
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The third order loop has zero steady state error for 
position, velocity , and acceleration inputs. The first non- 
trivial steady state dynamic error is produced by a rate of 
change of acceleration, J. Specifically the error due to a 
jerk, J, is given by 




8.4 uf 

125 (BW 3db> 3 


(3.63) 


Thus , 


BW 3dB > 


r 8.4 v -j 

-(125) (0.2) 


1/3 



(3.64) 


This implies that the maximum tolerable rate of change of ac- 
celeration at 20 GHz for a 1 Hz 3 dB bandwidth is roughly 0.024 

. 3 

ft/sec . 

Usually rates of change of acceleration on the order of 0.02 
r* _ 

ft/sec are not experienced with satellites in synchronous 

orbit, thus it may be possible using a third order loop to 

track the satellite motion using a loop 3 dB bandwidth on the 

order of 1 Hz. 

It is desired that even lower tracking bandwidths be used. 
This would provide capability for separating out phase fluctua- 
tions due to scatter under conditions when the Doppler spread 
associated with such processes is less than 1 Hz. Under these 
conditions it would be necessary to add rate aiding to the 
phase locked loop. Implementation of rate aiding requires a 
good estimate of the satellite ephemeris. The ephemeris data 
is processed to derive an estimate of the expected Doppler as 
a function of time at the receiver. This estimate is added to 
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the PLL loop compensation network output at the input to the 
loop VCO. Thus the loop need not track the satellite dynamics 
per se. In this case the loop need only track the difference 
between the estimated dynamics and the actual dynamics. In this 
way the use of rate aiding can reduce the dynamics to be tracked 
by at least an order of magnitude. 

In summary the brief discussion above has shown that the 
use of phase locked loop techniques for separating out the phase 
effects due to satellite motion and oscillator instability from 
phase variations due to scattering are marginal. To accommodate 
narrower Doppler spreads relatively sophisticated equipment would 
be required. Detailed evaluation of the technique requires speci- 
fic orbit characteristics. 
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SECTION 4 

CHANNEL MEASUREMENT USING PULSES 

It can be seen that frequency domain channel measurement 
techniques discussed in the previous section, involve certain 
difficulties with respect to oscillator stability. Because 
of the additional receiver complexity required to circumvent 
these problems, it is natural to investigate the potential of 
more direct measurement strategies. ■ In particular we will con 
sider certain alternative approaches which involve the trans- 
mission of pulse trains. 

4.1 Basic Approach 

In the following sections we shall examine the use of 
two different kinds of pulse trains; 

(a) Coherent pulses, i.e., individual members of 
the pulse sequence are phase coherent at least 
over time intervals large compared with the 
repetition period. 

(b) Incoherent pulses, i.e., because of the method 
of generation, adjacent pulses in the sequence 
are unrelated in phase. 

Both of these signals, when processed correctly, sample the 
complete channel impulse response at a given instant of time. 
With (a) , the coherent pulse train, we must take the same 
precautions as with the multitone spectrum to suppress the 
phase instabilities introduced primarily by transmitter 
oscillator effects. 

The basic pulse modulation approach is illustrated in 
Fig. 4.1. It will be assumed for the purposes of analysis 
that the pulse is perfectly rectangular, of length A, 
amplitude A, and is to be repeated every T secs. 
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Before discussing the two methods we first define the 
general time domain measurement problem. Using complex low 
pass notation as in Ref. [2.2] to symbolically represent the 
bandpass processes, we have a channel input— output relation- 
ship of the form, 

r (t) = J g(t,?) z (t-?) d? + n(t) (4.1) 

where g(t,5) is the complex time varying channel impulse re- 
sponse at time t to an impulse % seconds earlier, and z (t) is 
the complex envelope of the transmitted signal. n(t) is addi- 
tive receiver noise. An ideal channel measurement would yield 
in-phase and quadrature components of the channel impulse re- 
sponse as functions of time for each value of delay 5* These 
correspond to the real and imaginary components of g(t,£). 

From such measurements all of the necessary statistical charac- 
teristics can be estimated as discussed in Section 2.2. 

The rate of change of g(t,?) with respect to t, indicates 
the Doppler spreading properties of the channel, while the vari- 
ation with % indicates the delay spreading properties. Generally 
the degree of delay spread will be discussed in terms of param- 
eters such as L , the maximum extent of the delay dispersion, 
max 

or L, the RMS delay spread width. Similarly for Doppler spread 

we shall use B and B to denote maximum and RMS spectral 
max 

widths. (See Ref. [2.2] for a more precise definition of these 
quantities in terms of statistical properties of g(t,?)). It 
should be apparent after the discussion presented in earlier 
sections, that the measurement of g(t,£) is limited in resolu- 
tion by the signal bandwidth. In essence we require that 


where A is the pulse duration, and W is the corresponding 
signal bandwidth. To satisfy this relationship for the 
millimeter wave channel we need bandwidths on the order of 
1 GHz, but because of equipment constraints these conditions 
may be difficult or impossible to satisfy. 

A suitable signal for estimating the impulse response 
g(t,§) is a short pulse of width A, which is narrow compared 
with the delay spread L. Consider the signal. 


z (t) = Ap ( t) 


p (t) = 


0 | t 1 > A/2 

1 |t| < A/2 


< delay spread L (See Fig. 4.2). 


(4.2) 


Provided the variation of g(t,5) is slow in the t variable 
over an interval of length L, the response to a pulse at time 
t Q can be approximated using Eq. (4.1) as: 


r(t) = A J g(t Q ,§) p(t-§-t Q ) d§ + n ( t ) 


We can relate the observed sequence of responses directly to 
the function of interest, i.e., g(t,§). Thus let 


g. (t) = f g (t , §) P (?-t ) d? (4.3) 

t 0 0 

then 

r ( t ) = A g + (4.4) 
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Transmitted 

Pulse 


Fig. 4.2 Pulse Transmission and Response 


4-5 


Eq. (4.3) can be interpreted as a smoothing of the in- 
stantaneous impulse response g(tQ,§) by the transmitted pulse. 

If a sequence of pulses is transmitted with an interpulse 
spacing of T secs, the resulting sequence of responses resembles 
the quantity g(t,c) sampled at t=nT, except for the convolution 
with p(t). (See Fig. 4.3.) Obviously the period T needs to 
be chosen so that, T »L while we also require the pulse sam- 
pling rate to be higher than the Doppler rate; T < 1/B. Vari- 
ous statistical properties of the impulse response g(t,§) can 
therefore be extracted directly from the sequence of pulse re- 
sponses. Note that the coupling between the transmitted signal 
and channel impulse response demonstrated in Eq. (4.3) pre- 
cludes the possibility of obtaining useful channel data when 
the pulse width exceeds the delay spread L. Consequently we 
need, A <<: L. 

The sequence of pulse responses illustrated in Fig. 4.2 
allows direct measurement of the channel delay spread pro- 
perties, i.e., the behavior of g(t,5) with respect to %. 

Doppler information can also be extracted from the variation 
along the pulse sequence or alternatively by transmitting an 
auxiliary C.W. carrier either adjacent to or superimposed on 
the pulse spectrum. 

Most of this discussion has been concerned with the basic 
channel measurement theory involved; we now examine some of 
the difficulties involved in producing the required hardware 
for an incoherent pulse probing experiment. 

The most severe problem facing the system designer is 
that of obtaining the 1 GHz or so video bandwidth. This 
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aspect can be partially sidestepped by the use of a delay 
line processor shown in Fig. 4.4. 

The response r(t) propagates along the delay line after 
translation to an IF band. At the delay line taps, which 
are spaced by A seconds, we observe the processes (r (t-ni) ], 
so that each tap output is a 1 GHz bandwidth process around the 
IF center frequency. Following the detection of these waveforms, 
which can either be coherent or incoherent depending on our 
ability to realize the required equipment, the resultant 
low pass processes are sampled*. The sampling rate for each 
tap must be locked in with the pulse repetition frequency so 
that the same point on the pulse responses is sampled each 
time. i.e., the n tap is sampled for the k time at 

t = k (K q T) + nA 

T = interpulse period 
A = delay line tap spacing 
Kq= arbitrary integer 

There are various ways the sampler can be shared 

amongst the delay line outputs as long as the period T is 

much shorter than the inverse Doppler bandwidth 1/B, so that 

no information is lost. For example each T secs, a different 

th 

tap could be sampled in sequence, resulting in the n tap 


High speed samplers are currently available with rise 
times on the order of 30 picosecs. 
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Pulse Probe Receiver 








being sampled, for the k 


time at an instant. 


th 

t = (kN+n)T + nA 
n = 0, 1, ... N-l 
N = total number of taps. 

Thus at time zero we sample tap zero, time T + A tap 1, 

2T + 2A tap 2, etc. The sampling trigger signals can be 
generated from a clock reference with period T connected to a 
switched delay line of length nA. Note that the original 
sequences of output samples from each tap are of low 
bandwidth (i.e., the Doppler bandwidth), so that conventional 
digital filtering can be implemented, or a D/A converter 
can be used to regenerate analog signals. The IF processes 
have a noise bandwidth on the order of 1 GHz which means 
that the samples will be quite noisy. The post detection 
filtering suggested above helps to alleviate this situation 
as will be demonstrated subsequently. 

4.2 IF Signal to Noise Ratio Estimates 

Before continuing with implementation details for the 
pulse probe technique, it would be prudent to compute the 
SNR to be expected. 

The impulse response g(t,§) used in Eq.(4.1) to relate 
channel input and output quantities can be represented as 
the sum of two functions; 

g(t,§) = g c (§) + g R (t,§) (4.5) 


4-10 


where 9^?) is the time invariant coherent response, and 

g (t,5) is the time varying scatter component.* 

R 

Cbnsider the transmission of a single pulse of ampli- 
tude A at t = 0. From Eqs. (4.3) and (4.4) we have, 

r (t) ~ A Jt9 c (?) + p(t-S)d§ + n ( t) 


A A g (t) + n ( t) 


(4.6) 


First we compute the IF signal to noise ratio choosing 
the ATS-F experiment as an example. At 20 GHz the following 
parameters are typical; 


Transmitting antenna gain 
System and free space loss 
Receiving antenna gain 
Noise power/Hz 


37 db 
212 db 
56 db 
-159 . ldbm. 


Hence relative to a milliwatt transmitted, we have for a 
1 GHz bandwidth, an IF carrier to noise power ratio, 


Y = -212 + 37 + 56 - 90 + 159.1 
= -49.9 db 


Referring to Eq. (4.6) , we see also that a transmitted 
pulse is spread in time by the channel, hence resulting in 
an even lower SNR at the delay line taps. Let us assume 
that 9g(t) is a factor of 10 wider than the pulse p(t) 


In terms of previously defined transfer functions g is 

the transform of T while g is the transform of T T . 

c R c R 
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(resulting in N = 10 active taps on the delay line) . This 
implies a 10 db drop in signal to noise ratio. Conseauently 
a zero db y can only be obtained by transmitting a pulse 
which is approximately 60 db above 1 mW in power. Thus, 

3 

P =10 watts. 

T 

where P^ is the transmitted pulse power level. The energy 
of each pulse is on the order of 


— 6 

< 10 Joules. 

Note also that for a pulse repetition rate in the range of 
6 

10 pps , we require an average transmitted power of approxi- 
mately 



~ 1 watt 

The use of post detection filtering allows substantial 
improvement in the effective signal to noise ratio however, 
and this aspect will be discussed in the following sections. 

4.3 Incoherent Pulse Measurements 

In this section we assume that the pulse waveform p(t) 
has associated with it a phase angle which is random from 
pulse to pulse. This necessitates the use of incoherent 
detection. Here we consider the use of a square law detector. 
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Referring to Eq. (4.6) we can form an estimate of the 
quantity, 

p(t) = e( lgr Q (t) l 2 } (4.7) 

by computing, 

P(t) = |A g Q .(t) + n(t) | 2 - o n 2 (4.8) 

where a 2 = E[ln(t)| 2 } 

n 

and E denotes statistical expectation. 

This function indicates the estimated average power profile 
of the channel impulse response and as shown, is constructed 
from a single pulse response. In practice the estimate is 
of course refined by repetition of the pulse and super- 
position of the responses. For the above one shot estimator 
it can be shown that the estimate has the following properties. 
The average is computed to be; 

e{£ ( t) } = A 2 P(t) + e{ |n (t) | 2 } - Q 2 

= A 2 P(t) . (4.9) 

while the estimate variance is, 

E{(£(t) - A 2 P(t)) 2 } = E{<|Ag Q (t) + n(t)| 2 - a 2 ) 2 } 

4 2 

-Ap'(t) 

= a 2 (A 2 P(t) + o 2 ) (4.10) 

n n 
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This result was obtained using the assumption that the 
composite impulse response g^(t) is approximately time 
invariant, and hence non-random, i.e.. 


g(t, §) ~ g c (§) 


(4.11) 


The more general case which includes the random response 
g {t,s) can be obtained if desired, but leads to more com- 
plicated expressions. Roughly speaking the variance is 
increased over that due to noise alone because we are trying 
the estimate the average strength of a random process. As 
an example consider the situation where 

g(t,5) ~ g R (t, ?) {4 * 


i.e., the scattered signal component dominates the direct 
signal. If in addition we assume that g R (t,§) has the 
properties* , 

E[g R (t,? 1 )g*(t,? 2 )} = 6 (Sj.-^) Q (Sj_) 
E{g R (t,§ 1 )g R (t,§ 2 ) } = 0 


then the channel is defined by the function Q(?) which is 
estimated as in Eg. (4.8), i.e.. 


0(t) = |Ag 0 (t) + n ( t) \ 4 


- a 


n 


(4.14) 


See section 2.3 for an interpretation of these conditions. 
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where from Eq. (4.6) we have, 

g Q (t) ~ J g R (0»§) p(t-§) d£ (4.15) 

Note that, 

~ 2 rp * 

Et|g 0 (t)| } = JJ E Cg R (°'?i ), ? R (o, ? 2 ) }p(t-§ 1 )p(t-? 2 ) 

d? l d? 2 

= J Q(§) |p (t-?) ! 2 d§ 

~,Q (t) J |p(t -§) | 2 d§ 

= Q (t ) A (4.16) 


Eq. (4.16) results from the use of Eq. (4.13) and the 
assumption that Q(§) is smooth compared with the pulse 

A 

envelope. The estimate Q(t) given in Eq. (4.14) has the 
property; 


E{0(t) } = A 2 A Q (t ) 


(4.17) 


and variance, 

Var{6(t)} = E{(<$) 2 } - [E[^}] 2 

= Ett lAgJt) + n(t)| 2 - ] 2 }-A 4 A 2 Q 2 (t ) 

= E { I Ag Q (t ) + n (t) ! 4 }-2CT n 2 (a n 2 +A 2 AQ (t) ) 

4 4 2 2 

+ a n 4 - A A Q (t) 

= [A 2 AQ(t) + cr n 2 ] 2 (4.18) 
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To summarize we note that in both cases, i.e., the pure 

A A 

specular and the pure' scatter channels, the estimates P(t),Q(t) 
(Eqs. (4.8) and (4.14)) converge to within a constant of the 
desired quantity (Eos. (4.9) and (4.17)). Hence they can 
be made unbiased by renormalization. The variance in each 
case is better expressed as a fractional standard deviation; 
i.e. , 


st. dev, of estimate 
mean of estimate 


(4.19) 


[( 


_n 

A 2 P(t) 


) (1 + 


n 


A 2 P(t) 


-)]* 




pure specular (Eqs. (4.9) and (4.10)) 


[1 + 


n 


] 


V. A LQ (t) 

pure scatter (from Eqs, 


(4.17) and (4.18)) 


Note that we are generally interested in situations 
where the IF signal -to -noise ratio is low. The parameter 
y previously defined as the IF signal to noise ratio for 
a CW signal, can be expressed here as; 


Y 


A 2 / a 

n 


2 


It should be noted that the function Q(t) has the dimen- 
sions of power per unit time while P(t) has units of power. 
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Consider the consequences of averaging the estimate by 
filtering the samples as in Fig. 4.5. The effective stan- 
dard deviation index p given by Eq. (4.19) is reduced in a 
manner which depends on the time constant of the filter. 

For a sequence of samples spaced in time by seconds, 
and a filter bandwidth such that. 


w o > 


B 


B = Doppler BW. 


we average over a total of N samples each perturbed with 
independent noise components. N is given by. 


N = 


W T 
0 0 


It can be shown that the effective index p becomes, 

2 

^(W (~7T7) 7 -r « 1 


0 0' 'yP(t) 


A 2 P(t) 


- < 


pure specular 


i + (w t ) h ( : 1 


0 0 'yAQ(t) 


) 


pure scatter 


(4.20) 


Generally speaking the estimates involved will not be of any 
value unless p 1 can be made much less than unity. In the 
case of a pure specular signal this is achieved by making 

w 0 

( — ) very small. Similarly for the pure scatter signal. 

In this latter instance however we obviously require, 


W Q < B 
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Fig. 4.5 Incoherent. Pulse Receiver 
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and the performance is approximately. 


W 


= <— ) 
' B 


o.H 


(Vo 1 


( 


yAQ(t) 


) 


(4.21) 


If we require p '<1/10 and assume a unity SNR, i.e., y=l, 
it can be seen that when P(t) and Q(t) are 10 times wider 
than the pulse p(t) , we require on the order of 


V7 T 
0 0 


S 10 


samples for the pure specular case and. 


with 

-2 

W Q ^ 10 B 

for the pure scatter signal case. 

Figure 4.5 shows the structure of the receiver required. 
Note that the square law detection can be implemented using 
wideband microwave diodes biased correctly, although the square 
law characteristic could be replaced by some alternative non- 
linearity with little change in performance. The high speed 
sampler is required to have a gate time of between 30 and 100 
ps, performance which can be readily achieved using off the 
shelf samplers. Note that the 10 GHz or so bandwidth of this 
device provides the dual function of low pass filtering, thus 
eliminating the higher order harmonics from the detector output. 

4.4 Coherent Pulse Processing 

Now we take the case where phase coherence is retained 
between pulses in the transmitted pulse train. Such a signal 
could be generated by chopping a coherent CW source at the trans- 
mitter. The receiver is now required to extract in-phase and 
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quadrature signal components from the delay line taps, and 
can take the form shown in Fig. 4.6. The IF signal must 
first be mixed with in-phase (I) and quadrature (Q) references. 
The references are obtained from a transmitted pilot tone. 

The IF signal at each tap can be expressed as; 

A(t) cos Out + 8) 

The mixing operation then gives a signal part of the form, 

m^t) = A(t) cos (rut + 9) cos «jt 

m (t) = A(t) cos Got + G) sin tut 

The high speed sampler again acts as a low pass filter (BW ~ 

10 GHz) which results in sampled versions of the processes, 

X(t) = A (t ) cos 0 (t) 

Y (t ) = A (t) sin 9 (t) 

at the sampler output. Along the delay line we therefore 
observe in-phase and quadrature samples of the process r (t) , 
the response to a single pulse. Repetition of the pulse gives 
a sequence of sample pairs at each tap. (See Fig. 4.6.) 

Each sample contains an additive noise component which can be 
suppressed by filtering. If a filter of bandwidth B is em- 
ployed, the channel function g(t,5) will be approximately con- 
stant over the filter time constant, and the RMS noise level 
will be reduced by a factor; 

(T q B) 
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Hence the real and imaginary components of the impulse re- 
sponse will be observable directly with a signal to noise 
ratio obtained using an equivalent noise power 

% = °n ( V> 

g 

Pulse repetition rates on the order of 10 should be easily 
obtainable/ while a conservative estimate for the Doppler 
spread B would be 10 Hz. Hence the equivalent noise power 
is in the vicinity of; 


Considering the expected IF signal to noise ratio computed in 
Section 4.2, it appears that this method of directly obtain- 
ing in-phase and quadrature components of the channel impulse 
response will be feasible. The impulse response estimate is 
then processed further to obtain the complete mean, Doppler and 
delay characteristics. 

4.5 Summary for Pulse Probing Techniques 

The basic difficulty with this and most other millimeter 
wave channel measurements is that tremendous bandwidths are 
required to obtain useful data. The figure of 1 GHz has been 
arbitrarily chosen to be the signal bandwidth objective as a 
result of theoretical predictions and data from related experi- 
ments. Such bandwidths make video signal processing very dif- 
ficult and require state of the art components for generation 
of the -transmitted millimeter wave signal. The main non- ideal- 
ities and problem areas which must be considered in conjunction 
with the previous analysis are; 
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(a) Generation of sub-nanosecond pulses at relatively 
high power levels. 

(b) Arrival time jitter due to transmitter oscillator 
instabilities and motion of the satellite. 

(c) Accurate tracking of the pulse response arrival 
times which is necessary for extraction of average 
properties by overlapping the sequence of responses. 

The first two topics mentioned above will be discussed 
briefly. 

The state of sub-nanosecond switching technology is 
changing very rapidly, being stimulated to some extent by the 
prospect of high data rate communications systems. The funda- 
mental requirement is the availability of a high speed switch 
which must be used in amplitude and phase modulators and also 
in baseband data switching circuits. A wide variety of micro- 
wave semiconductor switching devices can be obtained with sub- 
nanosecond rise times. In particular we can consider the fol- 
lowing performance guidelines; 

Tunnel diodes 20-100 ps 

Common Emitter Transistor ~ 150 ps 

Step recovery diode 50-150 ps 

The primary difficulty is that of obtaining reasonably high 
power levels in conjunction with fast switching. However as we 
have seen in the previous sections, the use of high pulse repeti- 
tion rates allows somewhat lower power levels to be employed. 

Arrival time jitter is a manifestation of the oscillator 
drift and satellite movement (See Section 3.2.2). It is easy 
to see that for multipath spreads on the order of nanoseconds. 
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the accurate estimation and tracking of the pulse response 
sequence either using leading edge or center of gravity tech- 
niques, constitutes a critical receiver operation. Figure 
4.7 shows the consequences of incorrectly lining up a sequence 
of pulse responses used to form the impulse response estimate. 
Clearly the result is an over-estimate of the channel delay 
spread. Note that when the satellite moves only one foot dur- 
ing the averaging interval, a differential delay of 1 nano- 
second is introduced thus producing a relatively large error 
in the measurement of delay spread. 

Similarly the transmitter reference oscillator is not 
perfectly stable. Using the oscillator output 

f (t) *= COS (f») t + 0 ) 

s s 

the pulse train can be triggered by counting down to a repeti- 
tion rate of 1/T. 



The pulse generation reference is determined by the zero cross- 
ings of 


f T (t) = cos - (m s t + 0 s ) 


/ 2nt 0 s (t) \ 
= cos + — r ~) 


Following along the lines of previous discussion it is apparent 
that the properties of timing jitter A defined as. 
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Fig. 4.7 Pulse Response Averaging 
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can be found using the phase stability characteristics of the 
local oscillator. For example an oscillator stability of 1 
part in 10 averaged over 1 second implies 

_ ( 10~ 9 X 1 sec) ’“s = 1 nanosecond 


where A is the change in A (t) due to phase drift (over 1 sec), 
e 

This indicates the need for either a highly stable oscillator or 
some form of arrival time tracking. The one second averaging 
time relates to the time constant of the post detection filter 
which might typically be used to smooth the estimates derived 
in Sections 4.3 and 4.4. 

To conclude it can be seen that although the pulse probing 
technique appears promising, there are many areas which would 
require substantial investigation and development effort. 
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SECTION 5 


EQUIPMENT CONFIGURATION 

This section discusses the equipment configuration for 
measurement propagation characteristics of a 20-30 GHz satel- 
lite to ground link. The equipment includes a satellite-borne 
transmitter and a ground-based receiver and demodulator. The 
equipment described is basically similar to that employed in 
the ATS-5 and ATS-F experiments. The emphasis in the discus- 
sion below is placed on modifications of the current designs. 
These modifications are based on the experience obtained by 
NASA personnel with the current equipment, on improvements in 
the state-of-the-art of solid state components, and, in part, 
based on the results derived in this study. The proposed in- 
cremental phase correcting receiver is configured to be com- 
patible with existing ATS-F satellite 20-30 GHz experiment. 

5 . 1 Transmitter Configuration 

Figure 5.1 shows a suitable transmitter for generation 
of a multi-tone channel probing signal. The design takes full 
advantage of the current state of the art of solid state micro- 
wave signal sources. 

As shown in Fig. 5.1 the transmitter is controlled by two 
feedback loops. The inner loop is a phase locked loop which 
synchronizes the phase of the transmitted carrier component 
with a suitably multiplied reference signal. The second loop 
feeds back a demodulated tone which is compared against a refer- 
ence tone. The di fference between the two is used to modulate 
the microwave signal source. 

For the purposes of illustrating the operation of the 
transmitter, assume that the transmitted output is of the form 
given in Eq. 5.1. 
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Tone Loop . Fig. 5.1 Transmitter Block Diagram 









Where to is the microwave carrier frequency (say 20 GHz) and 
c 

ui) is the tone frequency (say 180 MHz). The deviation, 6, is 
set at roughly, 6 = 3.0, in order to produce four reasonably 
large- amplitude spectral components on either side of the 
carrier. The spectrum of the transmitted signal is approxi- 
mately as shown in Fig. 5.2. It is noted that for the purposes 
of the propagation measurement experiment the power in each of 
the probing tones need not be the same. It is only necessary 
that there be sufficient power in each of the tones so that a 
considerable margin is available above the receiver noise level. 
The transmitted power spectrum with the deviation set at 3.0 
yields three tones with almost equal power and a fourth tone 
which is not too far down. 

Part of the transmitted output is coupled to the feedback 
loops using a microwave coupler. The signal is then split 
between the tone loop and the phase sync. loop. The part of the 
transmitted output which is fed to the phase sync, loop is mixed 
with a carrier reference signal. This latter signal is derived 
by carefully multiplying up a stable reference signal which is 
say 5 MHz. Thus the multiplication factor, N, is on the order 
4,000. The result of the mixing process is filtered in an IF 
amplifier. The output of the first IF amplifier is mixed with 
a second reference at M x 5 MHz. The output is then filtered 
in a second IF amplifier. This IF filter is designed to select 
the carrier component at the output of the mixer and to reject 
those components which are associated with the tone modulation. 
The output of the IF amplifier is thus the down-converted micro- 
wave carrier component. This carrier is then divided by n in a 
frequency divider and the results compared directly against the 
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stable reference standard to create a phase error signal. 

Assuming that the transmitter output center frequency is 

s t 

higher than N x the stable reference frequency, then the 1 IF 
frequency in the phase synchronization loop is defined as 


u) 


IF1 


"J 

c 


N.o , 

ref 


(5.2) 


Where tu^^is the IF center frequency, is the K-band carrier 
output center frequency and ! o re ^ is the stable reference fre- 
quency, The second IF is similarly defined as 


- M'jj _ = (u - (N+M) to 
IF2 IFl ref c ref 


(5.3) 


The phase detector lock condition is defined by 

U). 


(0 


ref 


. 'U ~ (N+M) id 

IF2 _ c ref 

n n 


(5.4) 


Thus the output center frequency is defined in terms of the 
phase synchronizer loop parameters as follows: 


.0 = (N+M+n) m (5.5) 

c ref 

The phase stability of the output produced by the phase 
synchronizer loop is strongly dependent on the quality of the 
loop reference signals. In fact the stability of the output 
will be almost equal to the stability of the RF reference i.e. 
the stable reference carrier multiplied by N. 

The modulation of the output signal is also controlled by 
a feedback loop. Again a small sample of the output signal is 
coupled into a microwave discriminator. The discriminator de- 
modulates the frequency modulated waveform and produces a tone 
at the output whose amplitude is proportional to the modulation 
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deviation and whose frequency is equal to a tone reference fre- 
quency. This latter signal is derived directly from the stable 
reference standard by multiplying the stable reference carrier by 
k. The demodulated tone signal is compared against the 
tone reference and the difference is amplified in a loop com- 
pensation system. The loop compensation system actually em- 
ploys in-phase and quadrature demodulation of the tone to pro- 
duce in-phase and quadrature base band error signals. These 
signals are low-pass filtered and then are used to modulate the 
tone reference signal. The results are applied to the VCO con- 
trol input. 

The transmitter subsystem also provides a means for trans- 
mitting a sample digital message. This message is modulated 
directly onto the output carrier produced by the VCO. A micro- 
wave phase shifter is used for that purpose. The phase shifter 
is driven from a pseudo -random sequence generator which produces 
the sample data stream. For convenience the bit rate is chosen 
to be equal to the tone modulating frequency, 180 MHz. Digital 
logic which works at this rate is available off the shelf. Sim- 
ilarly a K-band digital phase shifter with sufficiently fast 
switching time is also readily available. Of course it would be 
more interesting to transmit at a significantly higher rate but 
components are not available at the present time for this pur- 
pose. It is possible to increase the chip rate to 500 MHz 
with a moderate effort. 

The transmitter configuration just discussed performs the 
same function for the most part as previously implemented sys- 
tems. Basic changes included in the configuration presented in 
Fig. 5.1, take advantage of new technology. The resulting 
device would be all solid state. The system of Fig. 5.1 relies 
heavily on the use of feedback techniques to assure that the 
desired output signal is being generated. This approach some- 
what reduces component performance requirements. For example, 
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the use of a microwave discriminator relaxes the specifications 
on the IMPATT VCO.In particular it is simpler to make a stable 
wide range linear microwave frequency demodulator than it is 
to make a stable microwave frequency modulator. The demodula- 
tor can be used to determine the output frequency with the aid 
of feedback. Similarly the inner phase synchronizing loop uses 
feedback so that a low power carefully generated reference sig- 
nal can be used to lock up the high powa: IMPATT oscillator 
thus transfering the stability from the low power signal to the 
high power signal. 

5 . 2 Receiver Design 

This section describes a receiver which is compatible with 
the transmitter described above and which provides incremental 
phase correction. The receiver is also compatible with the 
transmitter configuration currently employed in the ATS-F sys- 
tem. 

A block diagram of the receiver is shown in Fig. 5.3. The 

operation of the system is as follows. The received RF signal 

is converted to the first IF by mixing with a local oscillator 

reference. The local oscillator reference is generated in a 

frequency synthesizer subsystem which will be described in more 

detail below. For the moment it is noted that the first local 

oscillator frequency is coherently related to a fundamental 

* 

reference frequency denoted by •» . In particular the local 

* ■ r 

oscillator frequency is N . This down conversion is exactly 
the same as that performed in the phase synchronizer loop of 
the transmitter. Thus the first IF in the receiver is equal to 
the first IF in the transmitter's phase synchronizer loop. Ade- 
quate performance is achieved if this stage is wideband. To 
this end the first IF should be around 2 GHz. This choice is 
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similar to that exercised for the ATS-F receiver. The output 

of the wideband IF is used to drive a carrier loop and a series 

of tone demodulation circuits. The carrier loop extracts the 

carrier from the modulated signal and coherently demodulates 

to baseband. This operation is carried out in a series of 

down conversions. First the output of the wideband IF is 

mixed with a second coherent reference signal generated by the 

* 

frequency synthesizer subsystem. This signal is at M'u^ . The 
carrier component at the output of this mixing process is fil- 
tered in a narrow band IF and the results are phase detected 

* 

against a third coherent reference, at nil) . Note that these 
operations are essentially equivalent to those carried out in 
the phase synchronizer loop of the transmitter. Except that 
in the case of the transmitter the IF signal is divided by n 
whereas for the receiver it is better to multiply the phase de- 
tector reference by n. It is noted that in the transmitter we 
desire a wide pull-in range for the synchronizer loop and we 
are not specifically concerned with sensitivity of the phase 
detector since we are not dealing with a signal plus noise 
process. On the other hand, in the receiver we have a poten- 
tially fading weak signal with noise and we are interested in 
sensitivity. Furthermore, division of the fading process 
could lead to extra zero crossings which would upset the phase 
locked condition. 

The results of the phase detector process are used to 
control the VCO in the frequency synthesizer subsystem after 
suitable loop compensation as discussed in Section 3.3. 
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The carrier tracking loop provides two functions for the 

receiver. First it centers the received signal on all the 

filters in the tone demodulator system. More important, it 

★ 

provides a stable reference signal, m> , which is incrementally 
phase corrected. This reference signal is multiplied up to cre- 
ate conversion references for the carrier loop. It is also 
multiplied up in the frequency synthesizer subsystem to pro- 
vide incrementally phase corrected tone demodulation references. 

This operation is shown in the tone demodulation portion of the 
receiver in Fig. 5.3. 

Specifically, the output of the wideband IF is split eight 
ways. Each component is then bandpass filtered in an appro- 
priately tuned circuit to select a specific line in the received 
spectrum. The output of the bandpass filter is then mixed with 
a specific reference frequency which is produced in the fre- 

★ 

quency synthesizer subsystem. The resulting frequency is nio^ . 

This component is further filtered in an IF amplifier. Then 
the signal is synchroneously demodulated again using references 
generated by the frequency synthesizer subsystem to produce 
video signals which are proportional to the in-phase and qua- 
drature complex envelope components of the received tones. These 
outputs, sixteen video signals, contain the propagation data. 

A block diagram of the receiver frequency synthesizer is shown 

in Fig. 5.4. The operation of this device is straight forward. A 

VCO is controlled from the carrier loop compensation network. This 

★ 

VCO operates at the reference frequency, nuu^ . In order to be able 
to establish exceedingly narrow carrier tracking loop bandwidths 
the VCO must provide high inherent stability. It would not be 
surprising to find on further analysis that this VCO's instabilities 
place the lower limit on loop noise bandwidth. In other words it 


Quad 
H ybrid 
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Synthesizer 









is likely that the spaceborne master oscillator in the trans- 
mitter and the effect of satellite motion will be less than 
the effects of changing parameters in the receiver VCO. There 
are several methods available for realizing an extremely stable 
VCO. Further analysis should be applied to determine the most 
appropriate approach for this receiver. The VCO output is used 
to synthesize all of the required signals in the receiver. The 
manipulations required to produce these signals are relatively 
straightforward as can be seen in Fig. 5.4. The tone demodula- 
tion references are generated with the aid of relatively simple 
phase locked loops which create the required offsets. For ex- 

•k 

ample a demodulation reference at(M+k)’ii^ is generated as fol- 
lows. The VCO output which is running at the desired frequency 

is mixed against M times the reference frequency . The difference 

* 

is therefore at k x id . This is amplified in an IF amplifier 

r 

and compared against k x tu in a phase detector. The output of 

the phase detector is then used in the usual way to control the 

* 

VCO. The exact same circuitry is required to generate (M-k)-^ . 
In this case the VCO is constrained to operate below the first 
mixer reference rather than above it. In a similar way the other 
six demodulation references are generated. 

5 . 3 Compatibility with the ATS-F System 

The incremental phase correcting receiver just described 
can be implemented to be compatible with the existing ATS-F 
transmission format. The ATS-F system uses a 5 MHz master oscil- 
lator which is multiplied by a factor of 4,000 to produce a 20 
GHz carrier or by 6,000 to produce a 30 GHz carrier. A 180 MHz 
tone is generated by multiplying the master oscillator output 
by a factor of 36. In terms of the system just described com- 
patibility will be achieved if the following constraint is. 
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fullfilled . 


(0 = 4000 »o ; 

'0 

= 5 

MHz, <u 

= 20 

GHz 


c r 

r 

c 



to = 6000 id ; 

(D 

= 5 

MHz. m 

= 30 

GHz 

(5.5) 

c r 

r 

c 




Thus for operation at 20 GHz the various multiplication factors 
can be chosen as follows: 

N = 3600, M = 375, and n = 25. 

Fig. 5.5 shows a convenient synthesis technique for producing 
the three appropriate frequencies. 

It would convenient to use similar multiplication factors 
for operation at 30 GHz. In this case a suitable set is as 
follows ; 

N = 5600, M = 375, and n = 25. 

A suggested synthesis scheme for the 30 GHz operation is shown 
in Fig. 5.6. Note that this synthesizer shares a considerable 
number of components in common with the synthesizer proposed 
for 20 GHz operation. 
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SECTION 6 


DATA PROCESSING REQUIREMENTS 

This chapter is concerned with the data processing re- 
quirements for a multitone channel characterization experiment 
as outlined in Section 3. The recording and interpretation of 
received data is considered in relation to other auxiliary 
data such as rainfall rate. 

6.1 Processing of the Tone Responses 

The experiment outline given in Section 3.2 calls for the 
processing of each tone to obtain amplitude and phase charac- 
teristics. Alternatively in-phase and quadrature components 
can be extracted. This latter approach is of course closely 
related but has the advantage that phase ambiguity problems are 
completely eliminated. Figure 6.1 illustrates the data collec- 
tion procedure while Fig. 6.2 shows the way in which the data 
processing is carried out. 

After phase and amplitude distortion by the channel, the 
received tones essentially represent the transfer function H(f,t) 
sampled at frequencies [oj /2 tt] i.e. 



Thus the required properties of T (f) and T R (f,t) can be ex- 
tracted from the phase and amplitude behavior of the received 
processes {r (t) ] after translation down to individual IF chan- 
nels. We will assume here that the phase processes have 

been incrementally corrected for oscillator drift effects so 
that from Eq. (3.28); 
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Y (f , t) = (6 (f , t) + 9 (f , t) 
~ 0 (f , t) + e (f ) 


( 6.2 


where e(f) = nonlinear phase vs frequency component 

due to T (f) 
c 


jzf(f,t) = random phase introduced by scattered 

signal component T (f,t). 

R 

As shown in Fig. 6.1, the in-phase and quadrature components 
are extracted at each tone frequency, i.e. 


and 


X (t) = A (t) cos y (t) 
n n n 

Y (t) = A (t) sin y (t) 
n n n 


A n (t) “ l H fe' t ) 1 
Y n (t) =/(h' t ) 


(6.3) 


(6.4) 


As previously discussed, (X (t) 1 and £ (t) 3 are sets of 
stationary Gaussian processes; hence they are characterized 
by their means and cross covariances. This topic was discussed 
in detail in Section 2.3. The main distinction to be made here 
is the presence of additive noise generated at the receiver front 
end. The most interesting propagation conditions result in con- 
siderable attenuation of the signal, and as shown in Appendix B, 
the scattered signal may be 10 to 20 dB below the coherent or 
mean component; thus the scattered signal and receiver noise 
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component will often be comparable in power levels. Before pro- 
ceeding with this aspect we write the observed I and Q compo- 
nents in terms of their mean and random components, i.e. 


X^t) = X^(t) + x^t) + u k (t) 


Y k (t) = Y k (t) + Y k (t) + V k (t) 
where the mean components can be written. 


( 6 . 5 ) 


~ R et T C V2T, 


*)} 


U)_ 

y (t) = i It ^ 

k ml c\2rr 




( 6 . 6 ) 


and the random terms are. 


(U 

- R eK(^' *0} 


!0 

*k (t) - t )} 


( 6 . 7 ) 


^it) 


v k (t) 


= I noise in k IF channel 
(av. power = c^) 

til 

= Q noise in k IF channel 
(av. power = ct ) 


It can be shown that the processes {u^} and are a H i n< 3e- 

pendent of each other. Further individual time samples of each 
process will be approximately independent if the I and Q processes 
are sampled at tne Nyquist rate (i.e. W samples/sec where W is the 
nominal IF filter bandwidth) . 
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Note that in (6.6), the function T (f,t) differs from pre- 

c 

viously definitions because of the influence of the phase in- 
cremental correction receiver which completely removes the 
slowly-varying linear phase vs frequency component from the 
coherent signal transfer function. 

In the absence of other information pertaining to statistical 
properties of the channel, near otpimum estimates of means and 
covariances are provided by the following estimators; 

Mean X 



N 

E 

n=l 


\(nA) 


( 6 . 8 ) 


Covariance (x^t) and y^(t)) 


R xy (£ k' f /.' T) A \y t (T) 


“ N Ji (V ni >- (V nA - T > - h) 


(6.9) 


Here the sampling interval is denoted by A and the total ob- 
servation interval used in the estimate is T = NA. It is also 
assumed that t is some multiple of A. In the next section, the 
estimator performance is analyzed, and typical channel param- 
eters are used to provide concrete examples. 

6.1.1 Estimation Performance 

Consider first the mean estimate given by Eq. (6.8). The 
statistical average of the estimate is. 
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( 6 . 10 ) 


E{X } = E[X^ + x^(nA) + u^(nA)} = \ 


while the variance is, 


EU^ - \) ) = -J E{ (Z ^(nA) + u^fnA) ) } 
' N n 


- E { ( S V"*” 2 } + if 


( 6 . 11 ) 


with 


n=m 


u 


E{u (nA) u (mA) } =< 


otherwise 


Equation (6.11) can be further reduced to an expression contain- 
ing the correlation function for x, (t) ; i.e. R (t) 

k x k 


N 

Var (Xj^) = “ Z (N- 2n) R (nA) + 

N n=l *k ^ 


( 6 . 12 ) 


Note that when T=NA is much larger than the correlation time, or 

width of R (t), we can make the approximation, 
x k 


Var 


<V 


N 


Z R (nA) + a 
Ln \ 


! ] 


(6.13) 


For a Doppler bandwidth of B, R will have a width of approxi- 
mately 1/B, while the average power of x^(t) is given by 


a A R (o) 
X = *k 


Hence Z R (nA) ~ a /AB 

n *k 


(6.14) 


6-7 



Thus 


2 


(6.15) 


v * r <V = N [tS + 


Consider for example the range of parameters expected in the 
ATS-F experiment. [6.1] . 


Av. coherent power/tone = -110 dBm (no attenuation) 

Carrier/noise power density = 48.8 dB/Hz . 

For a 50 Hz IF bandwidth, i.e., A =1/50, we get. 
Carrier/noise power = 31.8 dB. 

Assuming a 12 dB attenuation and a scattered signal power 20 dB 

down from the carrier we arrive at a situation where the noise 

2 2 . 
power a , and the scattered signal power, , are both approxi- 
mately equal. Expressing the performance as a fractional standard 


deviation: 


Var(X k )| 2 


(— + 1 ) 
AB 


(6.16) 


With B ~ 1 Hz we can solve for the T = NA which gives a reason- 
able quality estimate. Thus for y = 0.1 we find 


io -2 ~ io _2 /nab 


since 


x , n -2 
. 2 = 10 


(6.17) 
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Therefore for the stated accuracy we need, 

T = 1/B 

= 1.0 sec. 

Note that the noise waveform contains an inherently larger 
number cf degrees of freedom for a given T, and as a conse- 
quence its effect is averaged out much faster than the scattered 
signal contribution. As a result, the averaging time suggested 
above represents only an approximate guide, the choice of T to 
meet a given performance being strongly dependent on the true 
scattered/coherent power ratio. 

In the second half of this section we consider the measure- 
ment accuracy attainable for an observation interval T when one 
is interested in implementing the covariance estimate given by 

a 

Eq. (6.9). The expected value of R is computed as, 

N 

E{R (t)1 = ^ E { £ [x (nA) + u (nA) ] Ty (n&- t) + v(nA-T)]} 
xy N L n ■> 

** R (t) (6.18) 

xy 

Since E{uy} = E[uv} = e{vx] = 0 

Note that the frequency subscripts have been suppressed to avoid 
complicating the notation. It has been assumed here that the 
interval T is much larger than the correlation time (hence also 
the range of t) , and that the means and are known exactly. 
In practice these quantities must be replaced (e.g. in Eq. (6.9)) 
by estimates which have non-zero error; thus the performance 
calculated below for known values, may be slightly optimistic. 

A 

The variance of R is found from; 

xy 


& 
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Var {*xy (T) } = E {~^2 ( E [*( nA) + u(nA)~|^y (nA-T)+ v(nA-T)j) 2 } 


R (t) 
xy 


(6.19) 


For Gaussian variables z^,k=l,4 we can use the property, 

e{|| z k } = E{z^z^) ^( Z 3 Z 4^ + E ^ Z 2 Z 4^ + E ^ Z 1 Z 4^ E ^ Z 2 Z 3^ 

k 

( 6 . 20 ) 

Applying this to (6.19) and using an observation time T which 
is long compared with the correlation time we get. 


var{S xy ( T )} ~ ± EE [r* ( (n-m) i)^] [ Ry ( (n-m) a)+ 6( (n-m) i)^] 


+ ~r T, T, R 2 (t) 
N 2 m n X Y 


+ — Y, T. R ( (n-m) A + tV * R ( (m-n) A + t) 
n 2 m n ) *y^ > 


BT (t) 
xy 


( 6 . 21 ) 


Assuming for example that. 


R (t) = R (t) 
x y 


R (t) = R (-t) 
xy xy 


we get. 
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( 6 . 22 ) 


Var K Y (t) } “ N {j R x (jA) + 2 V° >0 u + 

+ 5 { ?V^ T > R xy (jA+T >} 


This may be very roughly approximated by the expression; 

r* -i i R x (°) 2 4~i 

var Kv (T >/“ i 1 *i + 2 V o)o u + °uj 


(6.23) 


where 


0 .5^ 

P ' R x<°> 


and 


2 R x (0) 

S R ( j A) ra — — — B = Doppler spread 

x A3 

Following the previous definition of an RMS fractional error, 
we write. 


Var { R XV ( 

*>} 

R 2 (t) 
xy 


1 j 

_(i+p) : 

NR 2 (t) 
xy 

AB 

If (1 +p) 

R (o) 

f x ^ 

Nl AB 

\ 2 ) 
°u 


t 2 . 

x 


2 4 

+ 2R (o)ct + cr 
x u u 


] 


+ 2 


eft ♦ V(^) 

rt f n 


R ( t ) . 2 


u 


u 


(6.24) 

Consider the following parameters, which were used previously. 
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B = 1 Hz 


1/A = 50 Hz 
R ( °) 

— — = 1 (scattered power/noise power) 


R 

xy 


(t) 


_ 1 _ 

10 


R (o) 

x 


p ~ 1/10 

Y 2 ~ ^ |55 + 2 + lj 10 2 


For Y = 0.1 we thus require 

N = 5,800 
or 

T ~ 100 secs. 

Similarly when estimating R (t) and R (t), the same formula can 
be used with p set equal to unity. For example when 


R (o) ~ 10 
x 

The result is almost exactly the same. 

Note that when the scattered signal power level exceeds 
the average noise; i.e. 


R (o) 

x 


> 1 


u 


equation (6.24) is closely approximated by; 
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1 

N 


(6.25) 


2 

Y = 


R (o) 2 

(ifn) 

xy 


i.e. the error introduced by the additive noise is significant. 

2 

This is in fact true for even smaller values of R^(o)/a u due to 
the fact that all noise samples are independent while the signal 
is correlated over intervals of approximately 1/3. 

6.1.2 Summary of Tone Processing 

The in-phase and quadrature components of each tone response 
will be available in sampled form with sampling interval A. The 
quantities to be extracted for every T secs of data are; 

(a) Mean I and Q values. 


E[X k’ “ R e fT c (f k n 


E{Y, = I {T (f. ) 1 (6.26) 

k m^ c k 

(b) Time-Frequency Covariance Properties 
E^Xj- (t) x^ (t-T) } = R x (f k ,f^; t) 

E(x k (t) y / (C- t)) = R xy <W t) 

E (y k (t)y ; (t-T) } = R (fj^fj? t) (6.27) 

These correlation functions along with the mean signal 
levels given by Eq. (6.26) allow complete characterization of 
the channel. In many examples of communications channels involv- 
ing scattered signals, it has been found that the x and y components 
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have identical correlation properties and that certain symmetry 

in the cross correlation R can be exploited. However because 

xy 

of the narrow beamwidths involved in a millimeter wave link, 
such relationships can not necessarily be assumed, and before 
the design of a suitable data link can proceed it is necessary 
to experimentally verify the expected parameter values and their 
relationships. 

From the correlation processing we can obtain all required 
channel characteristics. For example the following parameters 
can be extracted. 

(a) Attenuation versus frequency. The average in-phase 
and quad, processes indicate attenuation. As discussed 
previously, variation across the band gives rise to 
waveform distortion. 

(b) Delay distortion. The components and Y n can be 
used to compute the phase shift across the band. 
Deviations from linear also distort the direct signal 
component . 

(c) The two dimensional covariance functions R , R and 

' xx xy 

R can be used to determine the Doppler and multipath 

yy 

spreads imposed by the scattering medium. (See Ref. 6.2) 

(d) Phase stability of the transmitter oscillator and 
satellite movement. The phase correction scheme allows 
these quantities to be separated from channel induced 
effects and measured if required. 

The above parameters may be correlated with available meteoro- 
logical data to verify that atmospheric models are useful in 
accurately predicting channel behavior. Then in turn, these 


6-14 


characteristics can be used to evaluate the distorting effect 
which the channel has on specific signaling waveforms and 
hence the corresponding short term error rate. 


6.2 Collection and Processing of Auxiliary Data 


The prevailing meteorological conditions obviously play 
a dominant role in determining the propagation conditions and 
resultant link performance. Consequently it is logical to 
undertake the sensing of certain meteorological paramecers 
while the channel characterization experiment is being con- 
ducted, so that postulated couplings between atmospheric con- 
ditions and link properties can be verified. Furthermore if 
these two aspects do demonstrate strong correlation, there 
remains the prospect of direct prediction of link behavior 
based on the local meteorological history. 

6.2.1 Meteorological Parameters of Interest 

We first single out some of the parameters which have a 
direct bearing on millimeter wave propagation conditions. Con- 
sider for example the following observable quantities. 

a) Rainfall rate? desired as a function of all space 
near beam intersections, but as a practical matter 
available only from a grid of spatial samples at 
ground level [6.3]. 

b) Atmospheric turbulence, in particular the atmospheric 
humidity levels. Again most readily available at 
ground level. r6.4], 

c) Refractive Index profiles (vs height) . These are of 
interest in determining beam refraction effects [6.5], 
and possible multipath propagation when low elevation 
angles are involved. T6.6]. 

In addition to the above categories it has been found that fog, 
cloud and snow have negligible effects on propagation at milli- 
meter wavelengths. 
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i 

f 

i 

f 

) 

The rainfall rate itemized in a) , has a strong influence > 

on attenuation. There has been a considerable amount of effort ; 

by others in this area [6.7, 6.8] and the basic relationships 
between rainfall rate and attenuation have been fairly well j 

j 

established. Similarly the mapping of rainfall rates over a j 

wide spatial area to determine the potential of diversity paths 
through a rainstorm has been widely discussed and is currently 
under investigation using a dense rain gauge network To. 3, 6.9], 

In light of these research efforts it does not seem profitable 
to consider an extensive rainfall gauge network implementation. 

It would still be desirable however to relate multipath and 
amplitude fluctuation data to the rainfall rate in the beam in- 
tersection volume and this can be simply satisfied using a line 
of gauges placed under the ground station beam. The characteristic 
structure of rainstorms tends towards cells of high rain inten- 
sity with dimension on the order of 1 km. Thus the spacing of 
the rain gauges should be 1 km or less r6.10].Fig. 6.3 illustrates 
the general configuration required. Rain gauges are currently 
available which will operate with high accuracy during rain rates 
of up to 10 inches/hour r6.111. The electrical output consists 
of a continuous tone, the frequency of which is related directly 
to rain rate. The carrier offset frequency is sampled at each 
gauge output and recorded on tape for later processing, including 
the conversion from frequency to rain rate. Fig. 6.4 illustrates 
the processing involved. A similar configuration of rain gauges 
has already been implemented for the ATS-V experiment r6.12" 1 . 

Atmospheric turbulence directly effects the amount of power 
scattered from the coherent beam and dictates the mean square 
amplitude and phase fluctuations of the signal. Unfortunately 
surface measurements of temperature pressure and humidity are far 
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Network 


Gauge Spacing < 1 km 
Total Number 10 - 20 


Fig. 6.3. Rain Gauge Configuration. 
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Fig. 6.4. Rain Gauge Data Acquisition. 


6-19 







from representative of atmospheric conditions at other alti- 
tudes. In general the complete profiles of temperature, 
pressure and humidity with altitude are necessary to establish 
the stability or relative turbulence of the atmosphere. These 
conclusions are largely substantiated by the results of various 
experiments where surface weather data was used. r 6.13, 6.141. 

In general the effects of high atmospheric turbulence are not 
severe, since the signal attenuation varies by only a few dB, 
and the link continues to operate under very high signal to noise 
conditions. This is in contrast to the high attenuation possible 
with rainstorms. Consequently it is not recommended that sur- 
face weather data be used in the experiment data reduction pro- 
gram, as the additional complications do not appear to be justi- 
fied. This does not mean that bandwidth and multipath param- 
eters are unimportant for propagation through the clear atmos- 
phere but simply that the range of variation is not great and 
only loosely coupled to observable meteorological variables. 

The third class of atmospheric data of interest is related to 
beam refraction effects which arise when the atmosphere becomes 
stratified into layers with different refractive index properties. 
In many instances these layers cover large horizontal areas so 
that vertical refractivity soundings provide a fairly complete 
picture of the total structure. These profiles are obtained 
from temperature pressure and humidity soundings. By utilizing 
a ray tracing program it is possible to examine the way a beam 
is distorted by the atmosphere and to investigate the possibility 
of multiple ray paths to the satellite. However it is well known 
that these refraction effects will be of minor importance in 
situations where the elevation angle exceeds 20°, such as would 
be the case for an equatorial synchronous satellite linked to 
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most points in the North American continent. r 6.5l. Conse- 
quently we will not pursue this topic further, but instead 
pass on to a brief discussion of the analysis of rainfall data. 

6.2.2 Verification of Link Parameters Using Rainfall Data 

Since rainfall causes the most severe link degradation 
and its gross properties are readily measured, we restrict fur- 
ther attention to this subject alone. The rainfall rate data 
can be utilized in two ways; 

a) Prediction of satellite to ground link properties 
such as total attenuation, coherent/scattered power 
ratio and delay spread. These predictions would be 
based on the techniques used in Appendix A and B, 

and since the true channel parameters would be avail- 
able by direct probing, correlation analysis of the 
predicted and measured data would indicate the valid- 
ity of prediction of link characteristics based on 
surface weather data alone. 

b) Interference Prediction. The intersection of antenna 
beams not associated with the same link will give rise 
to serious interference effects when a scattering 
medium such as a rainstorm is within the intersection 
volume. Predictions of the degree of interference can 
be made if the rainfall distribution is known. 

These topics which involve link parameter predictions are 
implicitly developed in Appendices A-D where analytical expres- 
sions for the various parameters may be found. 


6-21 



[ 6 . 1 ] 

[ 6 . 2 ] 

[6.3] 

[6.4] 

[6.5] 

[ 6 . 6 ] 

[6.7] 

[ 6 . 8 ] 

[6.9] 

[ 6 . 10 ] 
[ 6 . 11 ] 


SECTION 6 
REFERENCES 


Ippolito, L.J. , "The ATS- F Millimeter Wave Propagation 
Experiment," NASA X Document, X-751-71-460, October 
1971, p 3.12. 

Gallager, R. , "Characterization and Measurement of Time 
and Frequency-Spread Channels" , MIT Lincoln Lab. Tech. 
Report 352, April 1964. 

Semplak, R.A. , and H.E. Keller, "A Dense Network for 
Rapid Measurement of Rainfall Rate", BST J, Vol. 48, 

No. 6, July 1969, pp 1745-1756. 

F. Villars , and V.F. Weisskopf, "On the Scattering of 
Radio Waves by Turbulent Fluctuations of the Atmos- 
phere" , Proc. IRE, October 1955, pp 1232-1239. 

Crane, R.K. , "Propagation Phenomena Affecting Satellite 
Communication Systems Operating in the Centimeter and 
Millimeter Wavelength Bands" , Proc. IEEE, Vol. 59, No. 
2, February 1971, p 185. 

Ruthroff, C.L., "Multiple-Path Fading on Line-of-Sight 
Microwave Radio Systems as a Function of Path Length 
and Frequency", BSTJ, Vol. 50, No. 7, September 1971, 
pp 2375-2398. 

Ref. [6. 5], p 177. 

Godard, S.L., "Propagation of Centimeter and Millimeter 
Wavelengths Through Precipitation", IEEE Trans, on 
Antennas and Propagation, Vol. AP-18, No. 4, July 1970, 
pp. 530-534. 

Hogg, D.C., "Path Diversity in Propagation of Milli- 
meter Waves Through Rain", IEEE Trans. Antennas and 
Propagation, Vol. AP-15, No. 3, May 1967, p 410-415. 

Semplak, R.A. and R.H. Turrin, "Some Measurements of 
Attenuation by Rainfall at 18.5 GHz", BSTJ, Vol. 48, 
July 1969, No. 6, pp 1767-1789. 

Ref. [6.3] 


6-22 



[6.12] Ippolito, L.J., "Effects of Precipitation on 15.3 
and 31.65 GHz Earth Space Transmissions With the 
ATS-V Satellite", Proc. IEEE, Vol.59, No. 2, Febru- 
ary 1971. 

[6.13] Weibel , G.E., and H.O. Dressel, "Propagation Studies 
in Millimeter-Wave Link Systems", Proc. IEEE, Vol.55 
No. 4, April 1967, p 509. 

[6.14] Straiton, A.W. et.al., "Amplitude Variations of 
15 GHz Radio Waves Transmitted Through Clear Air 
and through Rain", Radio Science, Vol. 5, No. 3, 
March 1970, p 554. 


6-23 


SECTION 7 


CONCLUSIONS 

This section summarizes the results achieved during the 
study effort. The analysis of channel characteristics carried 
out in this program has attempted to extend and refine the 
propagation channel model. for K-band satellite-to-ground links. 
This augmentation of the channel characteristics is, in part, 
based on recently derived experimental data and, in part, on a 
more detailed definition of typical communication systems 
parameters hypothesized for K-band satellite-to— ground links. 

In a similar manner the report considered measurement techniques 
to derive channel characteristics using an experimental space- 
craft. Again, the design philosophy which is evolved in this 
report is based in part on recent experimental experience and in 
part on advances in the state-of-the-art in component technology 

7.1 Channel Characterization and Measurement Approach 

Channel characterization and measurement approach are dis- 
cussed in detail in Section 2. A relatively simple model for 
the channel is derived. This model is specifically appropriate 
for the analysis of communication and/or navigation systems oper 
ating through the propagation medium. The emphasis in the dis- 
cussion is placed on the effect of the channel characteristics 
on generalized received signals, rather than on the physical 
interactions which give rise to the channel properties. A 
measurement philosophy is derived based on the simple channel 
model. The rather general discussion of Section 2 is augmented 
by detailed analyses which are contained in Appendix A and 
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Appendix B. These latter discussions attempt to tie together 
the discussions available in the open literature so that some 
insight into the relative magnitudes of the effects to be 
observed can be derived. 

It was found, for example, that satellite- to-ground K-band 
links may not conform to the generalized channel models derived 
or assumed in the open literature. Note that there is a move 
toward narrower beamwidth antennas at K-band frequencies and 
above. It was found in this report that for narrow beamwidths 
the channel is not WSSUS as assumed in most theoretical dis- 
cussions (See for example Ref. 2.2 and 2.7). Further, exami- 
nation of the effect of narrow-beam antennas on channel char- 
acteristics uncovered a source of frequency dependence in the 
deterministic portion of the channel transfer function. 

The channel model derived in this report, and hence, the 
experimental measurement technique suggested , is valid in 
the region from 15-90 GHz. At 90 GHz, however, multiple scat- 
tering should become significant, and thus, the random scattering 
effects may become important and predominate over deterministic 
effects despite narrow beamwidths. 

7.2 Measurement Techniques 

Sections 3 and 4 considered two basic measurement techniques 
for extracting channel data. In particular, section 3 analyses 
the multi-tone technique which has been implemented for ATS-5 
and ATS-F. An incremental phase direction scheme is explored to 
augment the basic multi-tone channel probing approach. This 
technique is compatible with the receiver demodulation technique 
currently employed in the ATS hardware. The incremental phase 
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correction approach provides a little more data about the channel 
scattering characteristics. It is shown that the investigated incre 
mental phase correction scheme will provide adequate performance 
for the K-band satellite-to-ground link only if the com- 
bined effects of satellite oscillator instability and satellite 
motion could be tracked out. Chapter 3 concludes with a discussion 
of phase-locked loop techniques for incremental phase correction. 

In particular, the second and third order digital phase-locked 
loops are discussed which provide extremely narrow bandwidths . 

It is shown that such techniques will provide the necessary per- 
formance; however, detailed analysis to determine the exact level 
of system complexity and the applicability to the current satellite 
configuration have not been carried out at this time. 

The multi-tone approach analyzed and discussed in section 3 
does not provide data in a format which is directly usable by 
communication and navigation design engineers. Furthermore, 
it is shown that significant receiver complexity and processing 
is required by this technique to derive adequate data. Thus, 
in Section 4 we consider a more direct channel measurement 
approach. In particular, section 4 considered pulsed channel 
measurement techniques. It is interesting to note that although 
the pulse techniques appear to be more straightforward, more 
detailed analysis shows that these techniques suffer from essen- 
tially the same problems as were considered for the multi-tone 
approach, namely oscillator stability and satellite motion. 
Furthermore, there is a problem in the pulsed system concerning 
signal-to- noise re t: o. The problem arises because the satellite 
transmissions must consist of extremely short (sub-nanosecond) 
pulses and these are difficult to generate at high power levels 
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in the spacecraft environment. CW waveforms are generally 
considered to be more desirable in such an environment. Further- 
more the received pulses cannot be combined coherently at the 
receiver unless a highly complex receiver system is constructed. 
Thus, the results of this study reaffirm the conclusion that 
multi-tone processing, with all its drawbacks, remains optimum 
for K-band satellite-to-ground link channel measurement. 

7.3 Hardware and Data Processing Considerations 

Sections 5 and 6 consider equipment designed for channel 
measurement and data processing requirements, respectively. 

The equipment configurations described in Section 5 are 
based in part on experience obtained by NASA personnel with 
current equipment, and in part on improvements in the state-of- 
the-art of solid-state components. An all solid-state multi-tone 
channel prober transmitter is shown which relies heavily on 
feedback compensation techniques to assure the stable trans- 
mission of the prober waveform. An IMPATT oscillator is phase- 
locked in the transmitter to a suitably multiplied carrier 
reference frequency. The stability of the output is directly 
related to the stability of this reference signal. The second 
loop which includes a frequency discriminator is used to create 
the tone modulation. Again, the output of the transmitter is 
demodulated through the frequency discriminator and the resulting 
demodulating waveform is compared against the desired waveform. 

In this case the desired waveform is a simple sinusoid at the 
fundamental tone modulation frequency. In general, the feedback 
approach employed in the transmitter design somewhat reduces 
individual component performance requirements. A receiver design 
is described which is compatible with the feedback stabilized 
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transmitter. This receiver design also provides incremental 
phase correction. The receiver is furthermore compatible with 
the transmitter configuration currently employed in the 
ATS-F system provided a special frequency reference synthesizer 
subsystem is used. This latter device is also described in detail 
in Section 5. 

Section 6 deals with detailed processing requirements for the 
multi-tone channel measuring system. The discussion emphasizes 
the use of in-phase and quadrature processing rather than ampli- 
tude and phase detection techniques. An analysis is provided to 
demonstrate the advantages of the in-phase/ quadrature approach. 

Of course, the hardware changes required to convert from amplitude/ 
phase to in-phase/ quadrature are trivial. 

Using the processing techniques described in Section 6 the 
following channel characteristics can be extracted: 

a) Attenuation vs. Frequency and Delay Distortion 
These characteristics characterize the non-scattering 
distortion properties of the propagating medium. 

b) The Two-Dimensional Covariance Functions R »R and R 

aa av yy 

These correlation functions characterize the scattering 
components of the medium. 

c) Phase Stability of the Transmitter Oscillator and 
Effect of Satellite Movement. 

These characteristics constitute the third component 
in the K-band satellite-to-ground channel model. 

Section 6 also presents a discussion of the utility of various 
meteorological monitoring techniques, in particular rain gauge 
data systems are discussed in some detail. It is shown that a 
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relatively simple rain gauge configuration will provide adequate 
data to correlate observed channel characteristics with meteorological 

effects. 

7.4 Summary of Results and Recommendations 

The study reconsidered the channel characterization of K-band 
satellite-to-ground links. Several channel measurement techniques 
were considered for the purpose of collecting new channel data. 
Emphasis was placed on variations of the multi-tone probing tech- 
nique and some variations of a pulsed probing technique. It was 
found that the multi-tone technique is preferred because the state- 
of-the-art in component technology is not sufficient to support a 
pulsed techniques. Problems, in particular, were found with respect 
to available power levels and with respect to available signal band- 
widths with the pulsed techniques considered. 

The recommeded multi-tone technique is slightly modified 
relative to the approach implemented for ATS - 5 and ATS-F. Speci- 
fically, it is still uncertain whether an incremental phase correc- 
tion technique be used to track out the effect of transmitter oscilla- 
tor instability and satellite motion. If this could be done, then 
the multiple tones could be demodulated in such a way that actual 
phase differences among tones can be extracted if desired. It is 
noted, however, that direct phase measurements were not recommended. 

It was found that in-phase and quadrature demodulation of the tone 
components provided a small advantage. 

The transmitter design which is recommended reflects the latest 
technological advances. The approach will supply superior perform- 
ance in the spacecraft environment. Since the transmitter design 
is all solid-state, it should provide a significant improvement in 
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reliability relative to the current design. It is noted that the 
current design is based on technology available in the inid-1960 s. 

Some small component improvements are required to realize the 
transmitter as described in this report. In particular a small 
increase in power level above that which is available off-the-shelf 
at the present time is required. It is expected that since the 
state-of-the-art of IMPATT oscillators is advancing rapidly at 
this time that higher powers will be available shortly. 

It is difficult to evaluate the robustness of the conclusions 
derived in this report. Nonetheless it is stated that the most 
vulnerable area concerns the selection of the channel measurement 
technique. In particular, it is noted that given a significant 
advance in the state-of-the-art of digital logic components and 
RF high-power switching techniques, an alternate channel measure- 
ment technique would have been proposed. In such a case, the 
pseudo-random noise probing technique is more desirable. This 
technique does not suffer from s ignal-to- noise ratio limitations 
as does the pulse techniques of section 4. However, to implement 
the PRN approach, a bit rate on the order of 2 GHz is required. 

This certainly cannot be purchased off-the-shelf at the present 
time, nor is such performance anticipated in the near future. 
Furthermore, RF switch gear which could be used to create a 
phase shift keyed modulator is similarly unavailable at this 
time. Such switch gear would have to provide 100 picosecond 
switching transients in order to be useful. It is possible to 
realize a PRN channel prober which operates at roughly 500 
megabits per second, however such a rate is not sufficiently 
faster than proposed microwave data links to provide relevant 
performance information. 
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Several developments are recommended at the present time 
to ensure successful deployment of the equipment designed in 
this study. First it is recommended that an effort be made 
to improve the stability of the space qualified stable reference 
oscillators. For the system under consideration in this report, 
we are particularly concerned with stability over intervals on 
the order of 100 seconds or less. We are less concerned with 
very long term stability. 

Next is recommended that an effort be undertaken to develop 
space qualified IMPATT oscillators designed to operate in the 
region above 20 GHz. 

Finally, it is recommended that an effort be expended to 
develop efficient frequency synthesizer techniques for use with 
the proposed spacecraft transmitter portion of the channel prober 
system. 
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Appendix A 

RAIN SCATTERING EFFECTS FOR THE 
LOS MILLIMETER WAVE SATELLITE LINK 


A. 1 Description of the Scattering Process 
A. 1.1 General Discussion 

The discussion here will be limited to propagation through 
rain. Some of the properties of typical raindrops, which are 
assumed to be approximately spherical, are listed below. 

Diameter 1-3 mm 

Drop density 15 - 3000/m^ 

Rain rate 5 - 100 mm/hr. 

Considering a range of frequencies from 15 GHz (X = 2 cm) to 
94 GHz (X = 0.3 cm) it is apparent that the drop diameters will 
usually be less than a wavelength in diameter, while average 
spacing s 

_3 

d = (Density) meters 

will be significantly greater than a wavelength with the excep- 
tion of the highest drop densities at the lowest frequency. 

This latter conclusion is important because it allows the 
analysis to proceed in terms of independent scattering, i.e., 
the field scattered from any single drop can be studied without 
reference to other drops.* This does not preclude the possibility 
of multiple scattering of energy. However, it simply allows the 


Frequently it is possible to express the average scattered 
intensity as the sum of the intensities from individual scatterers 
but for a bistatic antenna configuration this assumption must 
be examined carefully. 
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analysis to proceed in terms of radiation fields with near field 
effects neglected. Note that the mutual spacing required for 
independent scattering is on the order of 3 times the particle 
radius (rather than a function of wavelength). 

Other media that are of interest at the higher frequencies, 
such as fog and clouds, also satisfy the criterion for indepen- 
dent scattering. 


A. 1.2 Field Scattered From a Dielectric Sphere 

The solution for the field which results when a plane wave 
is incident on a lossy dielectric sphere can be found in spher- 
ical coordinates. (See Ref. [Al]). 

The solutions are expressed in terms of the wave functions; 
,cos jzi 

*e =\ . ,P (cos 6)Z (KR) 

gn Ism 0 n n 


(A.l) 


where [P n (cos 0)} are the Legendre polynomials, and 


fj n (XR ^ J 2 KR J 


n+1/2 


( KR ) ; R < a 


Z (KR) 
n 


■< 


h n 2>(KR) H n^/2 (KR> ’ 


R > a 


= HanKel function of second Kind 
J = Bessel function of first Kind, 
a = Radius of the sphere . 

For large arguments, 


(A. 2) 


j n (kR) " KR COS (KR " (D I 1)ll) 

= _1_ [(.jjn+leJW + (j ) n+1 e-J kR ] 


(A. 3) 
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and 


h n 2)(KR) " kR (-J) n+lejKR 


(A. 4) 


The field at any point in space is the sum of two components; 


E = 


E i + E s 


E. + E 
1 t 


R > a 


R < a 


(A. 5) 


E = scattered field 
_s 

E fc = transmitted field (inside sphere). 

The incident field can be expressed in spherical coordinates as; 


E i = 1 x E 0® 


j(wt - kz) 


= V 


jwt 


£ <- J,n » <= 0 n + V 


(A. 6) 


In the far field (k 2 R » 1) with 1/R terms neglected; 


-0 X l 

n 


— : r j (KR)P (cos 8) sin 0 

sin 0 J n n 


A dP (cos 0) 

-Vn (KR) ae cos * 


(A. 7) 


. [RRj (KR)] 1 dP (cos 8) 

. n n v 

n = i — — u — t-t cos 0 

— e 8 KR d 8 

n 


1 [KRj (KR)] 1 

-i / — t” — - P (cos 0) sin 0 

0 sin 8 KR n 


(A. 8) 
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where K = *2 

A A 

, are unit vectors 
0 0 

[f(x)]‘ denotes derivative w.r.t. x. 
The scattered field is of the same form, 


E 

s 


V 


jwt 


E (-j) 


n 2n+l 
n( n+1 ) 


, s (2) , ( 2 ) . 

(a m. + jb n ) 
v n —On n — en 


(A. 9) 


n/ 2 ^ and n^ 2 ^ are computed from the expressions for m and 
-On -en (7l . un s 

n usinq h v ; in place of j , while the coefficients a and 
— en ^ n n n 

b® are computed from boundary conditions. It is apparent from 
the expression for E g along with the far field approximation for 
h^ 2 ^ ( KR ) that the scattered field is in the form of propagating 
spherical waves. 

A. 1.3 Description of Received Field 

For a plane wave Eb incident on a spherical particle, 
the resultant field observed at some distance from the scatterer 
can be expressed as (Ref. [A2] ) 
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where k = ~ 

x 0 

R = distance to observer 
S = scattering matrix 
a = particle radius 
0 = scattering angle 
n = refractive index. 

The two by two scattering matrix S gives a description of the 

scattered field at angle 0, for each polarization component 

(in the plane of scattering and orthogonal (see Fig . A. 1) ) . For 

particles large in terms of a wavelength the components of S 

are computed from a Mie series (i.e., Eq. (A.9)with boundary 

2 ira 

conditions). When — — is much less than unity, Rayleigh scat- 

K 

tering prevails and there is no scattering angle dependence. 

If a single scattering model is valid, the total received sig- 
nal can be computed from a summation over all particles. If 

e„ is the field contribution from particle l, the total field is 
l 



E\ + £ S^(n,a,0)E. 


-J KR t 

e 


(A. 11) 


( See Fig . A. 2) . 

The entries in the summation of Eq. (A. 11) can be rearranged 
into two groups, i.e., those within the first Fresnel zone and 
those outside of this zone. 
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Thus 


E_ - 


' _ - jkR 

M S (n,a, 0 )E . _ E e 

s i + E ~ ' u ser + 5 ? (n - a ’ e) -par- 

m J m K L l 




M S 
1 + I - 


m 


iKR 

m J m 


K 


+ E S 


E.e 

l 


-jkR 


l 


i < JKR <, 


(A. 12) 


(coherent term) 


(incoherent term) 


where the gross approximation has been made that all scatterers 
within the first Fresnel zone radiate coherently. Thus the 
first term above represents coherent radiation while the second 
corresponds to incoherent terms. The individual terms must of 
course be weighted by the antenna response (to plane waves) when 
computing the terminal voltage, e.g., for a linearly polarized 
antenna the voltage is ; 

-jk R . 

S.(e)E e * 

v r = g rI e c + E % (a) - — jk: — } < A - 1 3> 


where 


E 

c 


is the coherent component, 


E 

c 


M 

[1 + E 

m 


S 

jkR m ] E i 


(A. 14) 


and g ri (a) is the receiving antenna response to a plane wave 
K 

such that 


g R (°) = g r 


(A. 15) 
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Energy is removed from the incident wave by both absorption 

and scattering. These two components are accounted for in 

Eq. (A. 14) by recognizing that the refractive index n is complex. 

Consider first the nature of the coherent signal, i.e., 
contributions from within the first Fresnel zone. The summation 
over m in Ea. (A. 16) represents a complex constant which modifies 

* —4 

both the phase and magnitude of the original incident field E i « 
This results in: 

(a) An effective attenuation due to the intervening 
scattering medium (loss of energy from absorption 
and scattering). 

(b) A phase shift. 

The discrete distribution of scatterers which contributes to the 
coherent signal can thus be replaced with a uniform lossy medium 
having an equivalent refractive index; 


n 

e 


n + jn . 
er ei 


(A. 16) 


where n accounts for the losses introduced (see Fig. A. 3) • 
ei 

At first sight it might appear that the phase and amplitude 
of the coherent component fluctuate with the variation in number 
of scattering particles present. This is in fact true, but the 
effect is not very significant as is shown in Section B.2.4. This 
topic is postponed temporarily because conclusive arguments re- 
quire estimates of rain drop densities and these are contained 
in Appendix B. For most purposes we can safely assume the coher- 
ent signal equivalent refractive index to be time invariant. 

Outside of the Fresnel zone the scattered contributions 
to the received field tend to be incoherent, i.e., the param- 
eter R modulo X is uniformly distributed over X. That is to 
K* 
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Fig. A. 3 Equivalent Propagation Medium for Rain 



say, 

R = NX + v (A. 17) 

V 

where v is uniformly distributed on the interval [0,\], and X is 
the wavelength. Because the particles are in constant motion, 
the relative phase introduced by each is constantly changing and 
variable over 2tt radians. Consequently the superposition of 
these incoherent components gives rise to a signal with deep 
fading characteristics and average power determined by the sum 
of the individual scattered power levels. The fading rate is 
determined by the velocity of the particles and the geometry in 
question, while the increased path length from transmitter to 
receiver via a scatterer gives rise to the multipath phenomenon. 

To determine the relative importance of the incoherent and 
coherent fields one must examine the specific system geometry in 
question. It is easily seen that the. coherent contribution is 
derived from an integration over the first Fresnel zone, while 
the incoherent component is found by integrating over the remain- 
ing volume within the antenna beam. For millimeter wavelength 
antennae the Fresnel zone can be comparable to the 3 dB beamwidth. 
The current trends indicate future emphasis on higher frequency 
bands and larger diameter antenna structures, so that the effect 
being discussed will tend to become more pronounced. A more de- 
tailed analysis of the statistical properties to be encountered 
is deferred until the next section which does not use the rather 
simplistic approach of dividing the scattering region into co- 
herent and incoherent zones. 
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A. 1.4 Statistical Properties of the Rain Scattering Transfer 
Function 


The lowpass complex field response to a sinusoid of frequency 

^ j2TTft 

f is given in one form by Eq.(A.ll). (The carrier component e 
has been suppressed in this formulation.) i.e.. 


.2TTf 


E o = E. + E S (n,a,cp) E. 
R 1 l * 1 


~T 


,.2TTf „ 

o— > R i 


(A. 18) 


Including the antenna gain pattern g (ci.) and neglecting depolariza- 
tion effects, the complex received voltage can be expressed in 
terms of the transmitter terminal voltage as 


_. 2TTfR^ 

V R “ G R V T l 1 + T , W S t (n ' a 'V ) (A ' 19) 


where the constant G includes all required proportionality con- 

R 

stants. The subscript l denotes the l' th scattering particle 
with a. equal to the angle subtended by the scatterer with the 

'v 

transmitter-receiver path. (i.e., cp = 180- Ct) . Note that the 
quantity S will implicitly be a function of frequency because 
of the diameter/wavelength change in particles as a function of 
frequency. This aspect will really only be of great significance 
when the particle diameters are a significant portion of a wave- 
length. In this section we suppress the implied frequency 
dependence. 


The equivalent transfer function can be written in an alter- 
native integral form; 


-J- 


2TTf 


H ( f ) = oil + P g (a) n (r , a) 

X\ w X\ 


<^T R > 


dr 


} 


(A. 20) 
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where 


H = transfer function 

r = position vector 

R = Distance receiver to scatterer 

AR = Distance from transmitter to receiver minus 

the line of sight path length 

n = impulsive distribution function 

G now includes the linear phase term corresponding to line of 
R 

sight propagation. As a first approximation we assume n(r,ct) 
is a function of space only, i.e., the particle scattering pat- 
tern is broad compared with the antenna beamwidth. The distri- 
bution function is random in nature since it corresponds to 
the instantaneous distribution of scattering particles. It has 
mean and covariance properties expressible as 

E{n (r) } = n (A. 21) 


and 

E{(n(? 1 ) - n) (n(? 2 ) - n) } = o Q 6 (^ - r 2 ) (A. 22) 


Time variation of V and n can be included but this is not neces- 

X\ 

sary to illustrate the principles being discussed. In addition 
we note that n(r) and g(a) are approximately real, so that the 
in-phase and quadrature components of H(f) are found to be 


X(f) 


Re{H(f) 3 


G r( 1+ I g R {al 


n(r) cos (2TTfAR/c) 

2EE R 


(A. 23) 


with 


A- 13 


e{x ( f ) ) 



An (r) = n (r) - n 

the I and Q components x(f,t) and y(f,t) defined in Section 2.3.1 
can be written for a given instant of time as, 

x(f ,t) 
y(f,t) 

The complete time and frequency covariance functions for x and y 
are difficult to evaluate because of the need for time- frequency 

—t 

correlation properties of An(r) and these are generally not known. 
However, the frequency correlation functions can be established 


,-\ cos (2nfAR/c) 
in(r) — dr 


s in ( 2TTfAR/c) - 
2'fTfR/c r 
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Let 


x(f 1 ) = x { f ^ <t) 

y(f 2 ) = y(f 2 >t) 

Then 


vw = 

E{x (fj^) 

X(f 2 )} 

(A. 25) 

‘VW = 

E{x(f x ) 

y(f 2 )) 

(A. 26) 

VW = 

E{y (f 1 ) 

y(f 2 )) 

(A. 27) 


These functions can be evaluated from (A. 24) using (A. 22) 


» _ cos 2rrf Ar/c cos 2TTf Ar/c _ 

R X (f l' f 2 } = G rJ g R (Cl) ° 0 [ 2TTf 1 R/c ^ 2TTf 2 R/c - 1 ^ 

(A. 28) 


R (f. , f_) = g_ f (a) a [■ 
xy 1 2 R J o 


sin 217^ Ar/c cos 2rrf Ar/ c 


2TTf R/c 


-ID 


2TTf 2 R/ C 


] dr 
(A. 29) 


_ sin 2TTf AR/c sin 2nf Ar/c 

VW = G rJ 9 R (a) a o^" — — ] Idr 


2TTf r/c 


2TTf 2 R/c 


(A. 30) 


Clearly these functions involve arguments such as (f^ + an( 3 80 
are non- stationary in frequency. However the sum frequency com- 
ponents will be rapidly oscillating with respect to Ar and may 
not be important. To further illustrate the behavior of the func- 
tions involved, consider a one-dimensional version of the problem, 
i.e., using the Fresnel approximation let 
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(A. 31) 


2 

AR « — 

2Z 

(See Fig. A. 4 for notation) Note also that variations in R are 
relatively slow compared with other terms under the integrals. 
The covariance functions therefore involve integrals approxi- 
mately of the form: 


V £ l' f 2> ■ 

A ll 

/ x \ 
2 

o 

2 2 

COS (TTf ^X /c) COS (TT^X /c) dx 

(A. 32) 

In particular we find 

that 



R x (f l' f) 

A 1 I 

1 x \ 
g R Z 

o 

cos 2 (tt fn 2 / c ) dx (A . 3 3 ) 

R (f,f) = 

xy 

a 2 J 

/ x \ 
g R z o 

2 3 

sin (TTf-jX /c) cos (rrf^x /c) dx 

R y (f,f) = 

A 3J 

g R Z 

o 

2 2 

sin (TTf-j^x /c) dx 


9 2 , 

The functions sin kx^ and cos kx are sketched in Fig. A. 5 . 

Clearly if g is narrow compared with the first Fresnel zone we 
R 

have 

R (f,f) ? R (f,f) 
x y 

and 

R (f,f) ? 0 (A. 34) 

xy 

This result is important in the discussion at the end of 
Section 2. 
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Fig. A. 4. One-dimensional Scattering Model 


v- 
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Fig. A. 5. Integral Weighting Functions 
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Appendix B 

PREDICTION OF CHANNEL CHARACTERISTICS 
B.l Introduction 

The transmission medium between satellite and receiver 
will be subject to marked variations with changing weather 
conditions. Rather than distinguishing between different 
climatic conditions and their associated probability of oc- 
currence, we shall simply choose examples which might intui- 
tively be construed as worst case situations. Since the com- 
munication error rate will be strongly influenced by these 
abnormal low probability events anyway, they will be the most 
interesting as far as channel characteristics are concerned, 
i.e., the basic system design will include a suitable margin 
for fading effects. Below this fading level the error rate 
will be very high, while above it the error rate will be per- 
haps several orders of magnitude smaller. 

The two effects to be considered are rain scattering and 
atmospheric turbulence. The former will be the real culprit 
in terms of worst case system performance but the latter is 
included because it will be the dominant scattering mechanism 
most of the time. Where appropriate some attention will also 
be given to ionospheric effects and scattering from cloud 
particles . 

B. 2 Scattering From Rain 

As discussed in Appendix A the response to a transmitted 
sinusoid can be decomposed into the sum of a time varying 



coherent component and an incoherent component. Alternatively 
we can think of the received signal as the sum of a time in- 
variant specular component plus a random contribution. 

The following assumptions and approximations will be used 
in the analysis. 

(a) A single scattering model is sufficiently 
accurate . 

(b) Both the incident and first order scattered 
component fields will suffer from attenuation 
due to absorption and further scattering. 

If the intensity of the incident wave is 1^, and the dis- 
tance propagated through the scattering medium is Z (see Fig. 
B.l), then from Ref. [B-l], 


I (Z) = l.e _YZ 
SPEC x 


(B.l) 


where 


SPEC 


Y 


intensity of "specular" component 
(mean coherent field) 

extinction constant 


Note that 


Y = Y ab + Y s 


(B.2) 


where v accounts for energy losses due to absorption, while 
ab 

y indicates scattering losses. It can be computed by averag- 
s 

ing the scattering cross section a (8) over a sphere enclosing 
the scatterer. a(9), which describes the power scattered in 
direction 8 per unit solid angle for a unit volume of scatter- 
ers, is itself derived from the scattering matrix S(n,a,8) 
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(Ref. [B-l]) . An averaging over the drop size distribution 
is also required. Thus for a total rain layer thickness of T 
we have an emerging field intensity, 


"SPEC 


= I.e 
1 


- yT 


Neglecting further free space attenuation due to spherical 
wave propagation, this will be the strength of the specular 
signal at the receiving aperture. The received power from 
the- unscattered component is therefore; 


P - A I . e 

SPEC e 1 


-yT 



A 

e 


equivalent 
antenna area 


(B. 3 ) 


Consider now the scattered component of the field. At 
any point in the scattering medium, the intensity of the field 
incident on a particle is roughly, 

- yz 

I! = I.e Y 
1 i 

The power scattered from a particle then undergoes further 
attenuation as it crosses the medium to the receiver, due to 
both absorption and secondary scattering. Assuming none of 
these losses are regained through multiple scattering, the 
wave from the particle reaches the receiver with an intensity 
which is reduced b” a factor; 


-Y (T-z) 
e 


Thus for a unit volume of scatterers d^ meters from the re- 
ceiver, the flux density at the antenna aperture is. 


r. 

{ 

L 



4rrd 


a ( 0 ) 

2 

r 




- Y (T-z) 
e T 


(B.4) 


Hence the total average scattered power is obtained by inte- 
grating over the antenna power gain pattern G^_, i.e., 

G (0) = 9 r <e)G or with g r (0) = 1 


Then, 



Power density Antenna 

equiv. area 


(B.5) 


Note that this integral form could not be used if strong mul- 
tiple scattering was in evidence. Strictly speaking the inte- 
gration volume only covers the region outside the first few 
Fresnel zones, but generally we shall over-estimate the scat- 
tered power levels by using the complete beam intersection with 
the rain volume. 
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B.2.1 Scattered/Specular Power Ratio 


Let 

(i) 

(ii) 


us assume that we wish to find. 

The specular component power level p S p EC 
The ratio of scattered to unscattered power 

_ — p /p 

" SCAT SPEC' 


Consider (ii) 


first using Eqs, (B.3) and (B.5). 




P /P 
SCAT' 


SPEC 


■ J 

vol 


o9 r (e) 

T 

4rrd 

r 


du 


(B.6) 


Here it has been assumed that a is approximately independent 
of 0. Also 

d - (h + ~ - z) (B.7) 

r 0 2 

for narrow beamwidths. The integration volume is determined 
primarily by. the layer dimensions and the function g^CG). 

As a first approximation let the integration volume be cylind- 
rical with diameter determined by the beam thickness at the 
center of the layer. 

Thus 



du 


Tn(e h ) 
r 0 


(B.8) 
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is the 3 dB beamwidth. 


where e 

r 


But 


Hence 


For example consider the ATS-F experiment with the following 
parameters , 

(0.23° BW) 


Z = 46 x 10 m 6 /m 3 (46 x 10 4 mm^/m 3 ) 

This choice of Z corresponds to a 6"/hr. rainfall rate (see 
[B-3]) when the average drop size distribution is assumed to 
be that used by Laws and Parsons. The accompanying attenua- 
tion [B-3] at 20 GHz is 10 dB/Km and about 20 dB/Km at 30 GHz. 
Then the ratio p computed using Eg. (B. 10) is. 


-3 

e = 4 x 10 rad. 
r 

3 

T = 4 x 10 meters 


“ 


oTe 

3 

16 


(B.9) 


a = 


h ' k i 2z 


„ rr 


(Ref. [B- 2 ] ) 


T) = 


2 

e rr TZ 
r 


16X 


(B. 10) 
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(X in meters) 


1.84 tt 5 10‘ 15 
n = 4 

X 

Thus for f = 30 GHz, 

n = 5.8 x io" 5 
= -42 dB 

This ratio is of course smaller for lower rainfall rates 
(smaller Z) or lower frequency. 

The total attenuation effect on the specular signal as 
itemized in (i) , can readily be obtained from theoretical and 
experimental evidence relating attenuation to the parameter 
Z or rainfall rate [B-3, B-4]. Note that rainfall rates which 
give rise to more than 10 dB total attenuation through the 
atmosphere can be expected to occur less than 0.1% of the time 
at frequencies of 30 GHz. The example computed involves a 
4 Km thickness for the rain layer and a total of 40 dB attenua- 
tion, thus indicating that we are indeed dealing with an ex- 
treme case. 

We conclude this subsection by noting two assumptions that 
could possibly be violated. 

(a) When a Rayleigh scattering cross section is not ap- 
propriate (i.e., > .1) the equivalent Z may be significantly 

larger for a given drop size distribution and rainfall rate 
due to strong forward scattering. This would be important at 
higher frequencies such as 90 GHz. 


B-8 


(b) Multiple scattering may occur at the high frequency 
bands which would change the scattered power level slightly. 

(c) Experimental results for line of sight terrestrial 
links often involve larger beamwidths and larger paths which 
could conceivably give rise to a significant scattered signal. 

We now consider the signal dispersion effects which result 
from multipath scatter propagation. 

B.2.2 Multipath Spread Due to Rain Scattering 

Consider now the dispersion in arrival time of a short 
transmitted pulse due to the variety of different scattering 
paths constituent in the total scattered signal. When single 
scattering predominates the delay introduced by each scatterer 
can be computed from its spatial coordinates. Here we capital- 
ize on the symmetry of the antenna beams and use r to denote 
distance from the transmitter-receiver center line, and h to 
denote height above the earth. 

The delay introduced into an individual scattering path 
relative to line of sight can be expressed as, 

?(r,h) = (D 1 + D 2 )/c - t q (B. 11) 


where 



distance to transmitter 
distance to receiver 
line of sight transit time 
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Thus 


_ .,2 2 . 1 / 2 ^ , 

§ = [D 1 + (h + r ) ]/c - t 0 

2 

~ £ 

2 hc 


(B. 12) 


Since 


,, 2 2 1/2 « , A r 

(h + r ) - * + 2h 


The scattering elements have a broad scattering pattern so 
that the maximum delay experienced will be primarily deter- 
mined by the receiving antenna beamwidth rather than the ele- 
ments themselves. 

The following calculation is for the value of 5 computed 
at the 3 dB antenna beamwidth, and will correspond roughly to 
the 1/2 power value of the delay power spectrum, which spec- 
ifies average power as a function of delay. 

Let 


„ 2 r 
9 1/2 “ n 

= .23 


(3 dB beamwidth at 20 GHz) 


and 


h = 10 Km. 
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Then 


he 


1/2 


' 1/2 


8c 


-9 

= 0.061 x 10 sec 


(B. 13) 


Because the scattering medium is limited to the lower portion 
of the atmosphere, we would not expect h to exceed 10 Km, so 
that the above figure represents a bound on the expected mul- 
tipath, at least for a homogeneous medium-single scatterer 
model. Substantial increases in multipath spread would occur 
for a propagation path at low elevation angles (hence allowing 
a larger value of h) , or for multiple scattering situations. 


B.2.3 Doppler Spread Due to Rain Scattering 

The constant and random movement of the scatterers gives 
rise to a spectral broadening of any transmitted signal. 
Expressions for the Doppler shift in frequency due to motion 
of a single particle give a rough indication of the degree 
of spreading involved. Consider a scatterer at position (r,h) 
moving with velocity v. Using small angle approximations we 
can easily derive the following equations for Doppler shift 
in terms of velocity components either transverse or parallel 
to the transmitter receiver center line (see for example [B-6 ]) . 
For the transverse component v^. 


f 

r 


rv 
r 

x o h 


(B. 14) 
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But for the parallel component v , 


f h = 


2r 


. 2 
2r v. 




Using 0 = — we find that 
h 


v 9 
£ r - ^ 


V.0 

f ~-h- 
h 2 Xr 


«***•-. 


(B. 15) 


Since the height of the scatterer is not a parameter (within 
the range of these approximations) , we compute worst case 
Doppler shifts using only the 3 dB antenna beamwidth. Further- 
more we compute the Doppler shift per unit velocity (1 meter/ 
sec) since even typical velocities are difficult to predict. 
Thus 


-3 

0^2 = 3.8 x 10 rad 

X Q = 1.5 cm 

f,/v - 1.26 x 10 
h r 

-4 

f, /v - 4.8 x 10 
h h 


(.23°) 

(20 GHz) 
Hz/meter/ sec 

Hz/meter/ sec 


Vertical velocities of 10 meters/sec or less might be typical, 
while the RMS horizontal velocity (v ) due to turbulence will 
be on the order of several meters/sec. Based on the above 
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calculations then, we would predict a Doppler spread of 
about 1 Hz. Because there is a collection of contributing 
scatterers each with different velocity and different spatial 
position, we get a continuum of Doppler shifts at any instant. 

B.2.4 Coherent Propagation Through Rain 

Here we refer to all energy derived from a spatial region 
encompassing the first few Fresnel zones, and include both 
scattered and unscattered signals. It is assumed that the in- 
dividual scattered signals have identical phase, so that phase 
and amplitude fluctuations of the total signal in time and fre- 
quency, are derived from variable instantaneous number and size 
distributions . 

When single scattering is appropriate the equivalent med- 
ium refractive index can be represented as in Section A. 1.3, 


n 

eq 


n 


R 



To proceed with the evaluation of n^ as a function of fre- 
quency one needs to assume forms for drop size and number dis- 
tributions which adequately describe the range of distributions 
to be encountered. The actual computation is quite involved 
and beyond the scope of this study so we will merely rely on 
the results of other workers in this area. In particular con- 
sider the work of Crane fB-7i who evaluated n using measured 

o 

(convective) rainstrom drop size distributions. At 18 C in 

the region of 20 - 30 GHz he found that the refractivity 
6 

(n - 1) 10 was approximately linear, 

R 
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Rain rate 8" /hr . 


N = (n - 1) 10 6 
R 

( 7.0 - 0.3 (f - 20) 

1.2 - 0.05 (f - 20) Rain rate 2"/hr. 

with f in GHz. The effective phase shift experienced at the 
receiver is computed from, 

f(n - 1) x 

/ = — cycles 

c 

where x = thickness of rain layer. For example consider a 
5 Km thickness. Then, 

cycles 
8" /hr. 

2" /hr. 

again with f in GHz. Several points are evident from these 
expressions . 

(a) At any frequency, fluctuations in instantaneous rain 
rate from 2"/hr. to 4" /hr. give rise to phase variations of 
several cycles. 

(b) The linear phase shift with frequency gives rise 
to a time varying group delay. The change in delay in chang- 
ing from 2" /hr. to 8" /hr. rain rate is approximately 0.1 nsec. 


^(f) = 0.016 fN 

(o.016f (7.0 - 0.3 (f - 20)) 


0.016f {1.2 - 0.05 (f - 20)) 


(c) The quadratic phase vs. frequency behavior gives 
rise to distortion effects, i.e., a short pulse transmitted 
through the medium is spread in time (see Ref. [B-7]). 

Similarly one can predict the variation of attenuation with 
frequency and rain intensity using the imaginary component of 
the refractive index n . With linear approximations to Crane's 
data (Fig. B. 3 Ref. [B-7]) the following rough guidelines can 
be used; 

Attenuation at 20 GHz 

( 25 dB/Km at 8" /hr. 

7 dB/Km at 2" /hr. 

Attenuation gradient at 20 GHz 


As in the case of quadratic phase versus frequency, the vari- 
able attenuation across the band will cause pulse distortion. 

Generally speaking the instantaneous fluctuations in n 
will not be as extreme as the 2 "/hr. to 8" /hr. examples shown 
above. Since both n D and n can be shown to bear an approxi- 
mately linear relationship to the number of contributing rain- 
drops, we can gauge the RMS variation by drop number variations 
Assume for example that the number of drops is a 


\2 dB/Km/ GHz at 8" /hr. 


lO . 4 dB/Km/ GHz at 2 ,, Ar. 
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Poisson variable. Then, if M is the average number of drops 
which gives rise to a particular average rainfall rate, the 
RMS spread in n is proportional to «/M. However one must be 
careful to obtain the correct proportionality constant. The 
instantaneous phase, assuming a linear relationship, is, 

= am m = instantaneous number of 

drops in volume 


Then 


<^Vs ■ 2>2]1/2 


= a[(m - m) 2 ] 1/2 



where 



Eliminating a we obtain the expression 


(A*0 


RMS 



The average number of raindrops in the volume of interest will 
be quite large. For example a rain rate of 2"/hr. corres- 
ponds to a raindrop density of approximately; 

-4 

p - 10 drops/cc 
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At an altitude of 10 Km the \/4 Fresnel zone has a diameter 
of about 15 meters so that for a 5 Km thick layer of rain, 
the total volume is, 

3 

V - 5 x 10 x it x 49 cu.m. 


Thus 


M = p V 

—4 5 6 

- 10 x 7.5 x 10 x 10 


- 10 


8 


Hence 


‘^RMS * 2 X 1°" 

This is a small fraction of a cycle. We can readily conclude 
that any observable fluctuations in amplitude and phase with 
time for the coherent component, will be due to inhomogeneity 
in the medium structure rather than density variations as 
predicted using a Poisson distributed scatterer model. 

B. 3 Propagation Through the Atmosphere 

The atmosphere gives rise to two distinct propagation 
phenomena. These are: 

(a) Scattering of energy from the refractive index 
inhomogeneities which result from minute pressure, 
temperature and humidity fluctuations, 

(b) Distortion of the transmitted signal because of 
the weak frequency dependence of the mean refrac- 
tive index. 
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It should be noted that the atmosphere can be modelled as a 
layer with an exponential decrease in refractive index as a 
function of height. The refractivity value at ground level 
is on the order of 300 N units, i.e., the refractive index is, 

n = 1 + 300 x 10" 

Superimposed on this trend we have random fluctuations which 
give rise to scattering. The RMS variation is typically 
several N units. 

The next few subsections contain estimates of various 
parameters of interest, most of which have been obtained directly 
from the literature on this subject. 

B.3.1 Amplitude Scintillation 

The amplitude fluctuations of a CW plane wave source 
after propagation through a turbulent atmosphere, have been 
computed by several investigators, using the theory presented 
in [b- 9}. For example Lane [B-10] predicts the scintillation 
fading for a path through the atmosphere at the zenith angle 
using two different models. The magnitude of the received 
carrier can be thought of in terms of a fixed mean component 
Aq and a random component 6 a. The parameter usually computed 
from Tatarski's equations is, 

x 2 = [log A-j/Aq] 2 

« (6a) 2 /A q 2 for x 2 « 1 


where 



Generally attention is centered on the fading level corres 
ponding to an amplitude reduction of which corresponds 

to the RMS value of X for weak scintillation. 

Using the formulae given in Tatarski and two different 
models for the structure of the stmosphere. Lane [b-IOj has 
computed typical fading levels. At 35 GHz we find 

L = fading depth relative to signal without 
turbulence 

± 0.9 dB Model 1 (eq. (4) [C-10] ) 

2.2 dB Model 2 (eq. (5) [C-10] ) 

Experimental results at 10 GHz reported by Crane Cb— 11J indicate 
a general confirmation of these figures. 

B.3.2 Multipath Spread Due to Turbulence 

It was demonstrated in section B.2.2 that the multipath 
spread for rain can be computed directly from geometrical con- 
siderations, using the 3 dB antenna beamwidth to gauge the 
width of the delay power spectrum. For atmospheric turbulence 
the situation is somewhat different since the scattered energy 
received is confined to a cone much narrower than the beamwidth. 
Hence we would expect a multipath spread considerably less than 
the value computed in B.2.2. 
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APPENDIX C 

DELAY POWER SPECTRUM ESTIMATES FOR INTERFERING BEAMS 


C. i Delay Power Spectrum Calculations 

For a scattering cross section n ( A ) , the average received 
power for any ellipsoidal range shell can be expressed as; 


Q(«?) d? 


* p t I 

V(?) 


U3 G 
t r 

. 2.2 

d l d 2 


(4-rr) 


dv 


(C.l) 


where V( ? ) 

P t 

n 

h 




volume corresponding to a delay shall 
transmitter power 

scattering cross section/unit volume 
depolarization and absorption loss 
antenna gain functions 

distance to scatterer from receiver and transmitter, 
respectively 

scatter angle. 


Figure C-l shows some of the parameters involved. The antenna 
beam patterns can be expressed as the product of their maximum 
gains, G Qt and G Qr , with a renormalized pattern function. 
Assuming axial symmetry around the beam center line the antenna 
beam patterns can be expressed as a function of the angle sub- 
tended with the antenna center line ( Q-) , i.e., 


g^. ( o.) - GQ t gt ( °) (C.2) 

(similarly for G r , the receiving antenna). For example when g t 
is Gaussian shaped ; 
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, talk. 


g t (a) 



(C.3) 


with a 3 dB beamwidth given by 

e fc = CT t /0. 6 . (C.4) 

Note that in Fig. B1 the receiving and transmitting azimuth and 
elevation pointing angles are denoted by and 

C. 2 Simplified Integration Procedure 

Equation (C-l) i-s difficult to evaluate unless some simpli- 
fications are made. In particular we choose the following 
approximation (refer to Fig.c-2.) For the volume V( ? ) of interest 
it will be assumed that the integral Eq. (C-l) is simply the 
integral of a constant, i.e., the receiver beamwidth is suf- 
ficiently narrow so that changes in quantities such as n , dj, 

d„, and G. are negligible since the receiver beamwidth is on the 

^ ^ o 

order of 0.5 . Hence, 


o(n d? 


A o I 

v(c) 


nG^G^d 

2 2 
d X d 2 


dv 



(C. 5) 


Furthermore the receiving antenna beamwidth is sufficiently nar- 
row so that the volume V(f) in question can be approximated by 
the intersection of parallel planar surfaces at F and f + dF,' 
with the receiving \r. ;enna beam. Initially we shall consider the 
distribution of power as a function of distance along the receiv- 
ing antenna center line, i.e., let 
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Rec 


3dB Beam contour 


Fig. C.2 Delay Shell Intersection 


Trans 
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d i = D i + r 


then we can compute a function 


P(r ) dr = J ° r dV 

' d l d 2 V-(r) 


(C.6) 


much more easily than Q(p). V' (r) denotes the new integration 
volume shown in Fig.c.3. But 


J G r dv 
V' (r) 


2/2 


dr 


rf G.e 0 do djzi 


.2 


I *® 5 ,2 2 

dr 2 TT G t J p e a l°t d p 
0 


dr TT G t (d 1 rr t ) 


(C.7) 


Returning to Eq„.(c.6) we therefore find; 


P X 

P<r) = G ° tG ° r "' 0 - 6d i £r> 






r)e 


air) 


2 2 3 

d £ d 2 10 


'I 


(C.8) 


with r d l and d 2 in KmS 

X c(r) in meters 

i!r = angle separation from transmitting beam center line. 
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Note that it is possible to obtain a rough estimate of 
Q(P) from P(r) using a derived relationship between r and P, 
but its validity will be restricted to situations where the P 
plane almost intersects the beam at a right angle. 

The angle \ji can be expressed in terms of r and £ (see 
Fig. C.4) . 


r sin ( TT - ^ ) 

D 2 - r cos (n - 8 ) 


tan ,,r 


(C. 9) 


tlr 


tan 


f r sin P 
ID 2 + r cos 


Furthermore the antenna gains G Qt and G Qr can easily be expressed 
in terms of their beamwidths when the beamwidth is on the order 
of a few degrees or less. 

G(* ) = G 0 e“" 2/ ® 2 (C-10) 

For a radiation efficiency e; 

r = e G ( 0 ) 

0 jj G('>') dA/4^R 2 
R 

But 


[[ G(vli) dA * J*J G(|) pdp djz* 

R 



0 
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Path a : L a s a, + a 2 

Path b : L b = b, + b 2 



Fig. C.5 Differential Attenuation 
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For an efficiency of 0.6 this yields, 



0.6 

2 /„ 

a /4 


= 6.66/ g 2 



where e is the half power beamwidth and 


a = 0,6 e 

Returning to Eq. (C.8) with P^. = 1, we find that. 


(C.ll) 


(C. 12) 


P(r ) 


r.36^ 

64 ^ - 


6. 66 ^2-| | n(r)e 


-th 


d 




= .25 x 10 


_, lr 2 /2a| , \ 

-4 \ 2 e r rr(r)h^) 


€ 2 d 2 
r 2 
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C. 3 Polarization and Absorption Loss 

In general the polarization decoupling loss will be a com- 
plicated function of the antenna pointing angles. However the 
differential effect over the common volume will be quite small 
and here will be neglected completely. The power levels com- 
puted in this program will be based on minimum coupling loss, 
but corrections can easily be made to taKe into account the de- 
polarization effect. 
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The field incident at any point in the scattering volume 
will be attenuated relative to that which would exist in the 
absence of scatterers due to both absorption and scattering. 
Both effects have an exponential behavior with increasing path 
length for a uniform scattering medium, and will be lumped 
together, (see Fig.C.5) 

Let 


A = attenuation coefficient dB/Km 

then the intensity of the incident field I in terms of the free 
space intensity 1^ can be written, 

I = I 0 e“ aL (C. 14) 

L = total path length through scattering medium 
2 30 2 

a = (-— - q — )A attenuation constant. 

When the scattering medium is uniform and encloses the beam 
intersection completely; 


L = 


d l + 


d 2 “ 


D„ - D, 


= r + 


d 2 - 


(C. 15) 


Thus 



- . 23A( r + d 2 - D 2 ) 
e 


(C. 16 
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C.4 Scattering Cross Sections 


Assuming a single scattering model is valid, the scattering 
cross section o for the two most important scattering mechanisms 
can be expressed as; 

( i ) Rayleigh (rain) scattering 



with | K | ~ 1 


(ii) Atmospheric Turbulence 


o 


128TT 6 

X 4 


(0.033 


2/ 4tt 
c n ( ~T 


Sin 


i) 

2 


• 11 / 3 , 


(C. 18) 


where 0 = the scattering angle (change of ray direction at the 
scatterer ) . 

These expressions can be used in Eq. (C. 13) and the profile of power 
versus r is computed directly. With an approximately linear 
relationship between delay and r the resultant function P(r) can 
be interpreted in terms of power versus delay. Alternately 
when linearity cannot be assumed, more involved transformations 
can be performed. 


C. 5 Description of the Program SCATER 

This program computes path loss (includes antenna gains), the 
power profile P(r) from Eq. (C.13) , and transforms r into the equi- 
valent differential delay t relative to the intersection point 
of the two beams. The required input variables are listed below. 
All are in Kms or degrees unless noted. 
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(a) Frequency (f Q in GHz), Receiver beamwidth (BW), 
and a parameter ISIG; 


ISIG = 


r 


0 

1 


Rayleigh cross section 

Atmospheric turbulence 
(eq. (C. 18) ) 


(eq. (C. 17) ) 
cross section 


Note that in the program the following parameters are fixed. 


Z = 



K = 1.0 

c 2 = 10- 14 
n 

(b) Spacing (D) , Bearing angle to Nth— Sth line of trans- 
mitting site (BEAR), beamwidth of transmitting beam 
(EPS) 

(c) Scan coordinates. Radius (R) , azimuth (ALPH) , height 
of intersection (ZO). 

On completion the program reverts to a point where it accepts 
a new set of scan coordinates (i.e. , (c)). Input of zero inter- 

section height causes return to monitor. 

The total' path loss in dB. is computed from the sum of the 

samples of P(r). Corrections can be made to this value for absolute 

attenuation (for the beam intersection path) and variations in the 

2 

cross section parameters c and Z. 
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C . 6 Delay Power Spectrum and Path Loss Calculations 

The orthogonal baseline configuration previously considered 
was examined in terms of its delay power profile and path loss. 
These quantities were computed for different intersection points 
in the proposed scan volume with the following fixed parameters. 


f = 7.8 GHz 
BW = 0.5° 

D = 90 and 180 Km 
e = 1.45° 

Z0 (intersection ht.) =7.0 Km 

The resulting path loss values as a function of scan radius 
R, are shown in Figs . C. 6 and C. 7 . The former is with the atmospheric 
turbulence scattering cross section, Eq. (C.18), while Fig.c.7 
applies to a Rayleigh scattering medium (Eq. (C.17)). 

In the absence of differential attenuation effects (ATT — 0.), 
it was found that the power profiles P(r) (Eq. (C. 13) )were approximately 
Gaussian in form. i.e., P(r) is determined almost entirely by the 
transmitting beam radiation pattern. Thus we can approximately 
represent P(r) as; 

P ( r ) ~ f (Y r . V V a t ) 

i.e., a constant times the antenna beam pattern, where f( ) is 
a function of the antenna pointing angles only. 

The following points are also of interest concerning the 
behavior of path loss. 

(a) Along the great circle scan, path loss varies only 

with the change of common volume length for Rayleigh 
scattering and therefore follows an approximate cosecant 
law with respect to the receiving antenna elevation 
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Orthogonal Scan 
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Fig. C .6 . Path Loss vs Scan Coordinate R (Atmospheric Scattering) 


Loss, dB 
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Fig. C.7 Path Loss vs Scan Coordinate R (Rain Scattering) 



(b) 


angle X^.. For atmospheric turbulence there is a 
strong scattering angle dependence so that in the forward 
scatter mode losses are very much reduced. In the 
orthogonal scan configuration the scattering angle 
changes only slightly. (see Fig. C.6) Furthermore, 
there is no change of volume length L across the scan 
which accounts for the constant loss characteristic with 
a Rayleigh cross section. 

Choice of antenna spacing D. 

There are two ways of approaching this choice. One 
is for fixed transmitting antenna beamwidth e fc , and 
the other is for an illuminated volume at the receiver 
of fixed vertical extent. i.e., the beam thickness of 
the transmitting antenna at the receiver is given by; 


Ah ~ D£^. 


(see Fig. C.8) 


This determines the vertical extent of the scan volume. 
Consider variation of path gain with fixed assuming 
that the receiver beam is pointed in a vertical direc- 
tion. 

-2 

Gain a Common Vol. x D 

a Ah . D~ 2 
-1 

a D 


Hence the 3 dB improvement in Figs. C.6 and C.7 when the 
distance is changed from 180 Km to 90 Km. For fixed Ah; 
- 2 . 


Gain a 
a 
a 
a 


Ah x D x Antenna Gain 

-2 -2 

D x e. 

-2 2 
D x D 


const . 

Therefore the path loss can be maintained at a constant value as 
long as the beamwidth e^_ can be reduced sufficiently. 
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Figure C.8 Common Volume Dimensions 


C-18 


C. 7 Conclusions for the Proposed Bistatic Interference Experiment 


In designing an interference experiment there are two 
possible philosophies. 

a. Restrict the E.S. antenna beam elevation to greater 
than say 30° (see Fig.C.9) and measure interference through 
main lobe coupling or through coupling of the E.S. sidelobes 
with the TS beam. This represents a direct simulation of 
satellite link interference. The ES radiation pattern can be 
modeled (see Fig.C.10)by the sum of an omnidirectional pattern 
and a pencil beam. 

b. Recognize that sidelobe interference contributions 
while not as strong as main lobe effects, will occur more fre- 
quently. Hence use the E.S. main lobe over the complete range 
of elevation angles 0 - 90°, as a probing tool. This results 
in higher signal levels than would be available using the side- 
lobes for the same measurement as in a. (This approach has 
been used by Crane.*) 

The previous analysis and work by others, notably Crane, 
tends to support the following points. 

(i) The principal interference mechanisms will be rain, 
atmospheric turbulence and ducting. 


Predictions of Transhorizon Field Strengths Using Modeling 
Techniques. Lincoln Lab. TN. 
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Fig. C.9 Antenna Beam Configuration 
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Fig. C.10 ES Beam Pattern 


(ii) There is a 60 - 70 db difference between the inter- 
ference signal levels for rain and atmospheric turbulence 
for scattering through the main lobe (assuming E.S. eleva- 
tion angles greater than 30°). 

( iii ) Power contributed through the sidelobes will be important 
and will frequently exceed the CCIR recommended levels. 

All three mechanisms in (i) will be competitive, for low 
scattering angle, forward scatter situations. (i.e., if 
they do contribute significantly it will be from a preferred 
direction and elevation angle). 

(iv) At times ducting through the sidelobes can exceed 
main lobe rain scattered signal levels. 

Thus it can be seen that the real discrimination problem will be 
one of two kinds: Distinction between rain inside the main lobe 

and ducting through the sidelobes ( Fig. c. 12 ),or alternatively 
between rain, atmospheric turbulence and ducting, all through 
the sidelobes (see Fig. C. 11) .Note that (iii) indicates the prime 
advantage of a bistatic measurement over monostatic radar . In 
particular, ducting into the sidelobes will be a major inter- 
ference contributor, but this mechanism is not readily observable 
with a monostatic configuration. On the other hand, without an 
absolute delay reference the actual location of the dominant 
scattering medium will be unknown, and Doppler discrimination 
techniques which generally require such knowledge, will be in- 
effective. 

The exceptions, which fortuitously match our discrimination 
requirements are; 

(l) Ducting as opposed to rain scattering, will exhibit 
low Doppler shifts and spreads. 


3? 
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(2) Along and transverse to the great circle path, forward 
scatter from atmospheric inhomogeneities will result in 
negligible Doppler spreads (similarly for ducting) while 
rain is predominantly vertical in motion, and will contribute 
signals with substantial Doppler spreading. (Note that 
for sidelobe interference, the common volume is the inter- 
section of a TS beam with an omni beam. ) 

Conclusions . The orthogonal TS configuration originally proposed 
should be reduced to a single ES-TS pair. It would be desirable 
to have an E.S. beam steerable down to 0° elevation. In its 
normal operating mode (elevation greater than 30 ) it would 
be used to map rain intensity distributions, while in a forward 
scatter mode (elevation ~ 0°) ifc would give higher signal levels 
for investigation of sidelobe interference. 
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INTERFERENCE EXPERIMENT CONFIGURATION DESIGN 


D.l Outline of the Material Presented 

In view of the fact that the reader may not desire to 
follow the detailed analysis and design procedure presented, 
this Appendix has been organized in a slightly unorthodox 
way. Section D.2 contains an outline of the various param- 
eters involved for each mechanism (Section D.2.1) followed 
directly by a discussion concerning choice of configuration 
parameters (Section D.2.2). A detailed description of the 
program and computed data is then available in Section D.3. 

D. 2 General Objectives 

In choosing a suitable geometrical configuration, the 
basic objectives of the experiment must be kept in mind. 

(i) The accumulation of valid experimental data for 
the interference coupling between a simulated communi- 
cations satellite link and terrestrial microwave relay 
systems. (Provides justification for use of bi-static 
system over monostatic radar ) . 

(ii) Verification of the extent to which meteorological 
history for a given region can be used to predict inter- 
ference levels by extrapolating the results of direct 
coupling measurements. 

To achieve the second of these goals, it is necessary 
to make a correct assessment of the mechanism causing the 
interference and to obtain as much information as possible 
about its characteristics for correlation with observed meteoro 
logical parameters. Thus the system must be capable of carry- 
ing out the following functions. 

(a) Identification of the scattering mechanism with 
some degree of certainty. 
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(b) Measurement of parameters which will be important 

in characterizing the interference effects, e.g., average 
power levels, Doppler spreading and multipath spreading 
in the case of wideband communication systems. 

(c) Measurement of spatial structure of the scattering 
medium or equivalent parameters which can be related 
directly to meteorological data. 

In summary then, the system must perform an identification and 
subsequent mapping of a spatial region which corresponds to 
realistic beam intersection geometries. (See Fig. D.l.) 



■ TS2 

Fig. D.l Intersection Geometry. 



The spatial zone of interest will be defined arbitrarily in 
terms of an elevation angle greater than 30°. Very low ele- 
vation angles may be particularly useful in the investigation 
of specific mechanisms, however (see Appendix C ) to facilitate 
the mapping of a large spatial region without excessive dupli- 
cation of equipment, it appears essential to use the ES antenna 
in a scanning mode. This would provide a high degree of 
flexibility and also allow the extraction of large amounts of 
data from relatively rare but important events, such as high 
intensity rainstorms. A fixed beam experiment is restricted 
in the sense that only an isolated spatial region is sampled, 
with a resultant lack of confidence in the predicted 
long term statistics. 

It is interesting to note in passing that a scattering 
cross section mapping experiment as outlined above would also 
provide useful information for characterization of the satellite 
to ground line of sight link; i.e., given a complete descrip- 
tion of the spatial structure of the medium and estimates of 
defining parameters, the influence of the scattering medium 
for a line of sight configuration could be found. With such 
an approach long term statistics could be built up from the 
experimental data. 

D.2.1 Important Interference Mechanisms 

The most important scattering mechanisms that will be 
encountered and which must be separated, can be grouped as 
follows ; 

(a) Precipitation in the form of columnar or con- 
vective rain showers, snow, hail shafts. 

(b) Scattering from homogeneous turbulent atmospheric 
inhomogeneities. 

(c) Ducting effects and partial reflection mechanisms 
involving large scale horizontal statif ication of the 
atmospheric refractive index. 

(d) Reflections from aircraft and birds. 
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Table D.l MECHANISM ATTRIBUTES 
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Table D.llists some of the attributes of interest for each 
mechanism group. Ideally one would like to select an arrange- 
ment of terrestial station antennas, TS1 , TS2 ... etc., in a 
way that when all were intersecting the ES beam at a single 
point in space, individual TS beams were particularly sensi- 
tive to different attributes, (see Fig.D.l) This would allow 
identification of the existing mechanism for the mutual 
intersection volume in question. 

D.2.2 Configuration Parameters 

If attention is restricted to a scanning ES beam concept 
such as illustrated in Fig.D.l, the remaining system parameters 
to be selected are as follows. 

(a) Number of transmitting stations TS1 , TS2... etc. 

(b) Location of each in terms of distance D and angle 8. 

(c) Beamwidth of terrestrial stations (e). 

(d) Frequencies of operation (f Q ) 

(e) Probing signal bandwidth (W) 

Note that the receiving antenna (ES) beamwidth is assumed to 
be representative of a real satellite terminal (e.g. , 0.5 at 
8 GHz). Similar considerations should dictate the choice of 
transmitting site beamwidths. 

A computer program described in Section D.3 was used to 
give a preliminary evaluation of proposed configurations and 
the insight gained from sample calculations will now be pre- 
sented . 

(a) Number of transmitting sites . This choice is 
dominated by the desire for simplicity. We should 
choose the minimum number required to give an unambiguous 
identification of the scattering medium and also to 
allow the ES to scan the full volume proposed (i.e., 
the scan cone in Fig. D.l). 
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(b) Variation with D. Our analysis has demonstrated 
that over the scan volume of interest the spatial re- 
solution and vertical to horizontal Doppler component 
discrimination capabilities are only weakly affected 
by choice of D. The primary consequence of increasing 
D for a fixed beamwidth e, is the reduced power level 
at the common volume. 

A further consideration is the use of movable TS 
beams to allow variation in the intersection altitude 
(and hence scanning over a larger vertical extent.) If 
D is too large, shadowing effects due to earth curvature 
place a bound on the lower limit of the scan volume. 

(c) Beamwidth of terrestrial stations . The beamwidth 
should be representative of those used in microwave 
links. A wider beam gives a more extensive scan capability 
with lower demands on pointing accuracy and fewer adjust- 
ments of the TS beams to cover a given volume. (Note 

that 1.45° beamwidth corresponds to a beam thickness of 
2.3 Km at D=90 Km). The power levels expected for a given 
spatial resolution cell will be reduced as the inverse 
square of the beamwidth because of reduced antenna gain. 

(d) Frequencies of operation . The different scattering 
mechanisms exhibit wavelength dependences which could' be 
profitably used in their identification (using a multiple 
wavelength system). However the added complexity might 

be undesirable for an initial system design. This applies 
only to discrimination capabilities; accumulation of data 
at various frequencies would be useful but preferably 
over different observation periods. 

(e) Probing signal bandwidth . The analysis in Section d. 3 
demonstrates that for the scan volume proposed, the spatial 
resolution is determined almost entirely by W the bandwidth. 
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(i.e., it is almost independent of D the TS-ES spacing). 

From the example given in Table D-2. it is seen that 10 MHz 
bandwidth allows 190 meter resolution (D = 180 Km) in the 
extreme forward scatter mode, and 16 meter resolution in the 
extreme backscatter mode. As discussed, the received power 
levels are directly related to the volume in a spatial 
resolution cell, so that reduced resolution might be desirable 
if backscatter measurements are important. Alternatively a 
variable bandwidth scheme could be implemented, with resolution 
tailored to the scan position. 

A Proposed Configuration 

Based on the previous comments it appears that an ortho- 
gonal system of the form shown in Fig. D. 2 should be implemented. 
The parameters D (spacing), e (beamwidth) and W (bandwidth) 
cannot be determined without a more detailed study of signal 
levels. (See Appendix c) . 

Despite the fact that backscatter measurements will be 
quite difficult to implement, they will still provide data 
which is unique in terms of building up statistical knowledge 
of interference. Simply in terms of scanning the cone o^ 
interest (see Fig. D.2 ) , third antenna located at point B 
would be desirable, but it must be recognized that the time 
involved in a complete volume scan would be unrealistic any- 
way; hence the best strategy would seem to involve a broad 
scan in the forward scatter mode, with occasional limited 
sampling of backscatter conditions. 
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B.3.1 Description of the Configuration Analysis 

Consider an array of terrestrial transmitting antennas 
which have beams intersecting a narrow satellite earth station 
beam (ES) (See Fig.D.3). The terrestrial beams are denoted 
by TSl to TS4 for the system under consideration. The coordi- 
nate system is centered on ES while all azimuthal angles are 

measured with reference to an arbitrary line such as the Nth- 

Sth line. Each terrestrial station TS^ is located according to 

its radius D and angle 0 . (See Fig. D.4 for details of the 
n n 

currently configured RIPP experiment.) 

The five beams are assumed to intersect at a point deter- 
mined by parameters (R, a, h) in a cylindrical coordinate sys- 

tem. (See Fig. D.3.) A common intersection point for all 
beams is somewhat artificial in terms of a real system network, 
but is necessary for mechanism identification and to avoid over 
complication of the analysis. SIGNATRON has a computer program 
(SCAN) which gives a quick estimate of the delay and Doppler 
parameters for each coupling between the ES beam and terrestrial 

beam.* The following assumptions are made: 

(i) The ES beam is of neglibible width compared with 
terrestrial station beams . (Pencil beam). 

(ii) The volume common to both beams is completely and 
uniformly filled with elemental scatterers, i.e. 

The scattering cross section is not a function of 
.space. 


*Note that interference between terrestrial stations such as 
a T1 - T2 beam intersection can be analyzed using conven- 
tional tropospheric scattering prediction methods with the 
appropriate parameters. This will be considered if time 
permits . 
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Fig. D.3 Antenna Pointing Coordinate System 



Fig. D.4 Current RIPP Configuration 
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2. Norfolk 

22.7 

332.7 

3 . Ft . Lee 

86.0 

99.9 

4. Eastville 

50.8 

237.9 
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(iii) 


The gross parameters computed are found without 
considering the weightings due to scattering 
cross section, inverse path length etc. Only 
the intersection of the 3 dB antenna patterns 
is of interest. 

D- 3 . 2 Parameters Computed 

Figure D. 5 shows a detailed perpective of the common 

volume. The following parameters were computed for the up- 
per and lower intersection points. 



Fig. D. 5 Common Volume Geometry 
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Frequency shift in Hz per unit velocity 

(meter/sec) for x,y and z components of 
scatterer velocity. 


% Path delay, (psec.) 

L Volume length in Kms. 

t Diameter of ES beam at common volume mid point 


In addition the differential delay spread across the volume 
was computed from the upper and lower delays. i.e., 


* - s u - (Da> 

One should note that for a fixed temporal resolution the 
corresponding spatial resolution can be found approximately 

from, 


AL ~ (?) I (D ' 2) 

This assumes a linear relationship between delay and distance 
along the common volume cylinder, which is only approximately 
true. Furthermore, the slope of the ellipsoidal delay surface 
relative to the ES beam ultimately determines the characteristics 
of a spatial resolution cell. 

In utilizing the Doppler coefficients f x » f y» f z > two P oints 
are of concern. First the (x.y.z) coordinate system is rela- 
tive to the chord joining the ES and terrestrial station . 

Thus, for a network of terrestrial antennas the appropriate 
transformations of x and y components must be made. The 
second point is that at each location in the common volume 
there will generally be a spectrum of velocities for each 
direction. Hence there will be a corresponding spectrum of 
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Doppler shifts. The coefficients above should therefore only 
be interpreted in terms of Doppler sensitivity at a given point. 

The program operates in the following manner: 

(a) Read parameters R, a, z Q , f Q . e, 9, D 
r,a, z Q = intersection coordinates 
f = frequency in GHz 

e =. terrestrial station beamwidth 
D, 0 = coordinates of terrestrial station 

(b) Compute intersection coordinates in an (x,y,z) frame 
x = rCos ( a- 0) 

y = rSin(a-B) (D.3) 

2 = z o 


(c) Compute path lengths to center point (see Fig. D.6) 


Djl = (r 2 + z 2 )* 5 

D 2 = [(D - x ') 2 + y 2 + z 2 ]^ 


(D.4) 


(d) Find coordinates of upper and lower intersection 
points . 


1 D 1 2 +D 2 2 - d2 
e > = Cos [ ^^~ 2 — ] 


(D.5) 


(scattering angle at common volume center) 


= Tan ^(z _/r) (D.6) 

o o 

d u = D 2 Sin (e/2/Sin( P-e/2) (D.7) 

d^ = D 2 Sin (e/2)/3in(TT-p-e/2) 

r u .t = (D 1 - d u,t ) Cos <D - 8) 
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D.6 Intersection Geometry 
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Hence 


x 


u , l 


y 


u , l 


z 


u 


l 


u , l 

Cos ( a-3 ) 


(D. 

9) 

D l± 

d u,l )C ° S 

Cos (a- 8) 



u , l 

Sin ( a- 9) 


(D. 

10) 

D l± 

d u,i )Cos 

Sin ( a- 6) 



D l± 

a u,<.) Sin 

1 ^o 

(D. 

ID 


(e) Calculate the delay and Doppler components for 
upper and lower intersection points using the 
formulae above for and D 2 . i.e., 


lu 


D i < x u- v 2 u> 


Then for both upper and lower points; 


5 = 


D 1 +D 2~ D 


f - _. ° ( x _ (D-x) 
X " "c D 2 


) 


f ° (A- + -L) y • f = 

y c ^ D i d 2 ;Y ’ z 


The delay spread t is indicated by, 


etc. 



(D. 12 ) 
(D. 13) 


7 = - h (D - 14) 

In addition the total volume length and diameter are computed 
from 


L = d + d, (D. 15) 

l = D 1 x (BW) 

where BW is the beamwidth of the ES beam. 
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D.3.3 Evaluation of the Current RIPP Configuration 


A preliminary examination of possible configurations has 
indicated that an orthogonal system of antennas would simplify 
matters and yet still be sufficiently general. For this 
reason the scans shown in Fig.D.7 were evaluated using the 
current site configuration as a tentative example. The inter- 
section point was moved along the dotted lines corresponding 
to scan A and scan B at heights 10 Km and 3 Km (Freq. 3.67 GHz, 
B Widths TS1 , TS2 , TS3 , TS4 1.8°, ES, 0.62°. 

Table D. 2 contains the results. The Doppler coefficients 
at the upper and lower edges of the common volume were found 
to i>© within 10% of Gsch othor Blrnost without sxc©ption« 

Hence the table contains only the upper common volume Doppler 
coefficients . 

From a detailed study of the table one can draw the 
following broad conclusions. 

(a) The volume resolution AL specified in Eq. (D.2) will 
be of the same order of magnitude as the beam resolution l 
( see Eq. (D. 15) ) when the bandwidth W is about 5 MHz. For 
example, in a backscatter mode (see Fig. D.8a ) , this bandwidth 
gives AL ~ 40 meters (ES 30° elevation angle) while in the 
forward scatter mode (see Fig . D.8b), AL ~ 130 meters. An ex- 
cessively small value of AL compared with -t does not serve any 
useful purpose; on the contrary it is detremental in the sense 
that the measurement signal to noise is proportional to AL x L 
The volume resolution shows little dependence on spacing D for 
a backscatter mode and only minor variation with D in the for— 
wajcd scatter mode. In the backscatter mode the system performs 
almost as a monostatic radar would, so that bandwidth and 
resolution AL can be traded directly. 


D-17 




Figure D.7 Scan Patterns 
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Great circle path 

a. Backscatter Mode 




c. Vertical Mode 

Fig. D. 8 Great Circle Paths 
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Intersection ht 






























































































Intersection ht: lCkm 
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TABLE d- 2 (Continued) 






























































TABLE D-2 (Continued) 
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Intersection ht: 10 Km 































































'TABLE'D-2 (Continued) 

























































































































TABLE D-2 (Continued) 





























































(b) The total delay spread t becomes excessive for back- 
scatter modes of operation. This is only important however if 
the spatial resolution AL is correspondingly reduced to a value 
much less than l, which causes a severe loss in signal strength 
A large value of t is not too difficult to handle if a tapped 
delay line processor with sliding taps is implemented. Further 
more, it may be true that only a portion of the delay range 
corresponds to an interesting spatial region. 


(c) The ability to distinguish between vertically and 
horizontally induced Doppler effects is highly desirable, due 
to the predominance of z velocity components in rain, and x 
and y components in aircraft and atmospheric turbulence. The 
ratio of f ^/f is maximized by choosing to operate in a back- 
scatter mode, while it is minimized by operating in a forward 
scattering mode. Ratios of about 8:1 are possible in extreme 
cases. In the vertical mode (Fig. D-8c) horizontal and 
vertical velocity components cannot be distinguished, although 
a transverse velocity v^ will give zero contribution. (This 
is approximately true for all Great Circle Paths). 


(d) For the orthogonal path geometry shown in Fig. D.9, 
all three Doppler coefficients are found to be of the same 
general magnitude. There may be some advantage however in 
using both the great circle and orthogonal geometries in 
obtaining sensitivity to both of the orthogonal horizontal 
velocity components. 
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D.3.4 Evaluation of an Orthogonal Configuration 


As a test case the sample configuration shown in Fig. D.10 
was run. Two different values of D show the sensitivity in- 
volved in selection of this parameter. 

Parameters of interest are, 

D = 180, 90 Km 

e = 1.45° Transmitting beam widths 
BW =0.5° ES Beamwidth 
h = 7.0 Km Intersection height 

r was varied over a range corresponding to a minimum ES 
evaluation angle of 30°. The results are shown in Table D-3. 
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Intersection ht: 7 . 0 Km 
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TABLE D-3 (Continued) 




































































